Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 284,514 A→C intergenic (‑16236/+61990) unknown ← / ← unknown SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑Arg‑ACG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1284,5140AC98.7% 250.0 / NA 79intergenic (‑16236/+61990)unknown/unknownSSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑Arg‑ACG
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (30/48);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGAATACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCAAGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGTGGCGGTTCTGGAGGTGGCGGTGGTGGCGCTTCTGGTGGTAACGATTCTAGTGGTGCTGAGACTGAGAATTCTG  >  1/284267‑284760
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
tatGAATACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTc                                                                                                                                                                                                                                                     <  1:286504/251‑1 (MQ=255)
     aTACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGG‑GAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               <  1:813436/251‑1 (MQ=255)
      tACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               <  1:306027/251‑1 (MQ=255)
      tACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               <  1:54656/251‑1 (MQ=255)
      tACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               <  1:355145/251‑1 (MQ=255)
        ccTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCAAACCATTGACGAGGCGCAAAACATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGGCCGGGGGCAATGACATCCAAAGCGGATAATGTGATATTTTAGA‑CCGCCGGGGCACGTCTTAcccc                                                                                                                                                                                                                                              >  1:629857/1‑249 (MQ=255)
              tGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAAACATTACTGTAGGCAAGGAGTATCAGTTAATCGCCAATAACGGGCCGGTGTCAATGACAACCAAGGCGGGTAATGGGGTATTTTAGA‑CCGCCGGGGCACGTCTTACCCTATTTaaa                                                                                                                                                                                                                                       >  1:25782/1‑250 (MQ=255)
               gCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCAAAACATTGACGAGGCGCAAAACATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAAAGGGCCGGGGGCAATGACATCCAAAGCGGGTAATGTGATATTTTAGA‑CCGCCGGGGCACGTCTTACCCTATTTAAGa                                                                                                                                                                                                                                      >  1:55209/1‑251 (MQ=255)
                cTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     <  1:240091/251‑1 (MQ=255)
                cTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     <  1:604256/251‑1 (MQ=255)
                cTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     <  1:596517/251‑1 (MQ=255)
                cTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAAACATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGGGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     >  1:656621/1‑251 (MQ=255)
                cTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGACAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     <  1:834518/251‑1 (MQ=255)
                  cAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     <  1:619987/249‑1 (MQ=255)
                         cGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAatat                                                                                                                                                                                                                            <  1:464649/251‑1 (MQ=255)
                          gTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               <  1:756546/231‑1 (MQ=255)
                           tGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGGGAGGAAACCCAAACCATTGACAAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAGCGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTc                                                                                                                                                                                                                                                     <  1:471447/224‑1 (MQ=255)
                           tGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTc                                                                                                                                                                                                                          <  1:814863/251‑1 (MQ=255)
                             ccACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCAAACCATTGACGAGGCGCAAAACATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGGCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGGCTTACCCTATTTTAGaa                                                                                                                                                                                                                                     >  1:59216/1‑238 (MQ=255)
                              cACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAcc                                                                                                                                                                                                                                                <  1:855305/226‑1 (MQ=255)
                              cACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAAc                                                                                                                                                                                                                       <  1:661780/251‑1 (MQ=255)
                                          cGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAGAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTggaagggaag                                                                                                                                                                                                            <  1:874321/250‑1 (MQ=255)
                                          cGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTggaagggaag                                                                                                                                                                                                            <  1:2363/250‑1 (MQ=255)
                                          cGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTggaagggaag                                                                                                                                                                                                            <  1:477543/250‑1 (MQ=255)
                                             gAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCCAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCCATAACGGTCCGGGGTCAATGACATCCCAGGCGGATAATGTGATATTTTAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAAAGGGAATATTCAAc                                                                                                                                                                                                                       >  1:844927/1‑236 (MQ=255)
                                             gAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCAAACCATTGACGAGGCGCAAATCATTACTGTAGGGAAGGATTATCAGGTAATCGCCAATAAAGGGCCGGGGTCAAAGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGGGGACGTCTTACCCCATTTAAGAAAGGGAATATTTAAATGGAAGGGAAGGCgg                                                                                                                                                                                                        >  1:484255/1‑251 (MQ=255)
                                                   caAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTTACGAGGCGCAAATCATTACTGTAGGCAAGGAATATCAGTTAATCGCCAATAAAGGTCCGGGGTCAATGACATCCAAAGCGGATAATGTGATATTTTAGACCCCCCGGGGC‑AGTCTTACCCTATTTAAGAACGGGGAAATTTCACTGGAAGGGAAGTCGgtgggtt                                                                                                                                                                                                  >  1:866153/1‑247 (MQ=255)
                                                   caAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTc                                                                                                                                                                                                                                                     <  1:814443/200‑1 (MQ=255)
                                                        aaaCCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTaa                                                                                                                                                                                              <  1:632322/250‑1 (MQ=255)
                                                        aaaCCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTaa                                                                                                                                                                                              <  1:264906/250‑1 (MQ=255)
                                                        aaaCCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTaa                                                                                                                                                                                              <  1:522404/250‑1 (MQ=255)
                                                          aCCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATgg                                                                                                                                                                                           <  1:754688/251‑1 (MQ=255)
                                                          aCCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATgg                                                                                                                                                                                           <  1:217607/251‑1 (MQ=255)
                                                              tCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAg                                                                                                                                                                                       <  1:460613/251‑1 (MQ=255)
                                                                       aacggAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCAAACCATTGACGAGGGGCAAAACATTACTGTAGGCAAGGATTATCAGTTAATCGCCAAAAAAGGTCCCGTGTCAATGACATCCAAGGCGGGTAATGTGATATTTTAGA‑CCGCCGGTGGACGGCTTACCCCATTTAAAAAAGGGAATATTTAAATGGAAGGGGAGTCGGGGGGTATTAATGGTAAGCTCATTGa                                                                                                                                                                               >  1:298933/5‑250 (MQ=255)
                                                                             aaGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATgg                                                                                                                                                                                           <  1:613014/232‑1 (MQ=255)
                                                                             aaGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATgg                                                                                                                                                                                           <  1:71305/232‑1 (MQ=255)
                                                                             aaGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCAAAACATTTACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAAAGGGCCGGGGTCAATGACATCCAAGGCGGGTAATGTGATATTTGAGA‑CCGCCCGGGCACGTCTTACCCTATTTAAGAAAGGGAATATTTAAATGGAAGGGAAGGCGGGGTGTATTTATGGTAAGCTCATTGACTTGaaa                                                                                                                                                                        >  1:691165/1‑250 (MQ=255)
                                                                               gTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAAACATTGACGAGGCGCAAAATATTACTTTAGGCAAGGACTATCAGTTAATCGCCAATAACGGTCCGGGGTCAATGAAATCCAAGGCGGATAAAGTGATATTTGAGA‑CCGCCGGGGCACGTCTTAACCTATTTAAGaa                                                                                                                                                                                                                                     >  1:194143/1‑188 (MQ=255)
                                                                               gTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAAACCCAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAAAGGTCCGGGGTCAATGACATCCAAGGCGGGTAATGTGATATTTGAGA‑CCGCCGGTGCACGTc                                                                                                                                                                                                                                                     >  1:57234/1‑172 (MQ=255)
                                                                                ttAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATc                                                                                                                                                                       <  1:895733/249‑1 (MQ=255)
                                                                                    tgtCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGt                                                                                                                                                                  <  1:742697/250‑1 (MQ=255)
                                                                                                  gTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTAttt                                                                                                                                                                                                                                          <  1:460596/164‑1 (MQ=255)
                                                                                                   ttATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGTTATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGt                                                                                                                                                                  <  1:763608/235‑1 (MQ=255)
                                                                                                     ataAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACtt                                                                                                                                                                            <  1:274959/223‑1 (MQ=255)
                                                                                                          gaagTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               <  1:461292/151‑1 (MQ=255)
                                                                                                             gTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTAt                                                                                                                                                                                                 >  1:26778/1‑194 (MQ=255)
                                                                                                               ggCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAg                                                                                                                                         >  1:827034/1‑248 (MQ=255)
                                                                                                                    gagaGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTAACAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTc                                                                                                                                  <  1:13600/250‑1 (MQ=255)
                                                                                                                            aCCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTc                                                                                                                                                                                                                          >  1:581553/1‑154 (MQ=255)
                                                                                                                                 aaCCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTaa                                                                                                                                                                                              <  1:588377/177‑1 (MQ=255)
                                                                                                                                   ccATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGggc                                                                                                                                                             >  1:443535/1‑208 (MQ=255)
                                                                                                                                   ccATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGAt                                                                                                                  >  1:408864/1‑251 (MQ=255)
                                                                                                                                              gcgcAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGa                                                                                                       >  1:786825/1‑251 (MQ=255)
                                                                                                                                                 cAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTAtt                                                                                                                                                                                                <  1:373689/159‑1 (MQ=255)
                                                                                                                                                    aTCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAgg                                                                                                 <  1:29691/251‑1 (MQ=255)
                                                                                                                                                     tCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAGGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGa                                                                                                                >  1:612125/1‑235 (MQ=255)
                                                                                                                                                          aCTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAAc                                                                                                                                                                                                                                    <  1:159875/114‑1 (MQ=255)
                                                                                                                                                                ggCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATc                                                                                                                                                                       >  1:810786/1‑169 (MQ=255)
                                                                                                                                                                 gCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCgg                                                                                     <  1:144845/250‑1 (MQ=255)
                                                                                                                                                                      gATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGaa                                                                                                                                                                                                                                     >  1:200788/1‑101 (MQ=255)
                                                                                                                                                                      gATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCt                                                                                                                        >  1:300576/1‑210 (MQ=255)
                                                                                                                                                                        ttATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTAccc                                                                                                                                                                                                                                               >  1:303252/1‑89 (MQ=255)
                                                                                                                                                                            cAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGc                                                                               <  1:72960/245‑1 (MQ=255)
                                                                                                                                                                                      ccAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTc                                                                                                                                                                                    >  1:673053/1‑134 (MQ=255)
                                                                                                                                                                                      ccAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGGGCTGGTGGCGGGTCttgg                                                               >  1:31284/1‑248 (MQ=255)
                                                                                                                                                                                         aTAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAg                                                                                                                                                                                                                                       <  1:795303/80‑1 (MQ=255)
                                                                                                                                                                                          tAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGGCGGTGGCGGTCCTggggggg                                                           >  1:292268/1‑246 (MQ=255)
                                                                                                                                                                                               gTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTAtt                                                                                                                                                                                                <  1:462417/113‑1 (MQ=255)
                                                                                                                                                                                               gTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGTGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGTGGGGGTTCTgggggggggggg                                                      >  1:831770/1‑241 (MQ=255)
                                                                                                                                                                                                          aTGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGGGCTGGTGGCGGGTCTGGGGGGGGCggggggggggcctt                                           >  1:360290/1‑239 (MQ=255)
                                                                                                                                                                                                                      gCGGATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGTGGCGGTTCTGGAGGTGGCGGTGGTGGCGCTTCTGGTGGTAACg                                <  1:433220/250‑1 (MQ=255)
                                                                                                                                                                                                                         gATAATGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGGGCTGGGGGCGGGTTTGGGGGGGGCGGGGGGGGGGCTTCTGGTGGGAACGagtt                            >  1:783711/1‑248 (MQ=255)
                                                                                                                                                                                                                             aTGTGATATTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTc                                                                                                                                  <  1:511077/145‑1 (MQ=255)
                                                                                                                                                                                                                                  atatTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGGGCTGGGGGCGGGTCTGGGGGGGGGGGGGGGGGGGCTTCtggtggggaaggttttagtggtgg                   >  1:635457/1‑233 (MQ=255)
                                                                                                                                                                                                                                    atTTGAGA‑CCGCCGGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGa                                                                                                                                                                          <  1:693962/98‑1 (MQ=255)
                                                                                                                                                                                                                                             ccgccgGTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGTGGCGGTTCTGGAGGTGGCGGTGGTGGCGCTTCTGGTGGTAACGATTCTAGTGGTGctgagactgag         <  1:670984/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   gTGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGGGGCGGGTCTGGGGGGGGCGGGGGGGGGGCTTCTGGGGGGAACGGTTCTAGGGGGGCTGGGACTGAGaaat     >  1:532842/1‑247 (MQ=255)
                                                                                                                                                                                                                                                    tGCACGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGTGGCGGTTCTGGAGGTGGCGGTGGTGGCGCTTCTGGTGGTAACGATTCTAGTGGTGCTGAGACTGAGAATTCTg  <  1:719087/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
TATGAATACCTTAGTGCTCAATGACCGTGCCACTAATGTGGGCGGGAATCACAAAGAAACCGTCAAAGGCGAACAGAAAGTTAGTGTCGGCAAAGCCCGTTATAAAGAAGTGGCGGGAGAGGAAACCCAAACCATTGACGAGGCGCAAATCATTACTGTAGGCAAGGATTATCAGTTAATCGCCAATAACGGTCCGGTGTCAATGACATCCAAGGCGGATAATGTGATATTTGAGA‑CCGCCGGTGCAAGTCTTACCCTATTTAAGAACGGGAATATTCAACTGGAAGGGAAGTCGGTGTGTATTAATGGTAAGCTCATTGACTTGAATCCGGGTGGGGCGGCGACATCGAGTGCTGAAGCTGATTCTGCGCCGGCTTCGGATGACGGTGGCGACGGAGGTGCTCGTGGCGGTGGTGCTGGTGGCGGTTCTGGAGGTGGCGGTGGTGGCGCTTCTGGTGGTAACGATTCTAGTGGTGCTGAGACTGAGAATTCTG  >  1/284267‑284760

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑