Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,137 T→A 7.5% E9V (GAA→GTA) ‡ sdaC ← serine:H+ symport permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,1370TA 7.5% 514.7 / 23.7 160E9V (GAA→GTA) ‡sdaCserine:H+ symport permease
Reads supporting (aligned to +/- strand):  ref base T (65/83);  new base A (5/7);  total (70/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 168 total aligned reads displayed.

GCAGCGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTT  >  NC_002947/4080988‑4081286
                                                                                                                                                     |                                                                                                                                                     
gcagcGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCgttgtt                                                                                                                                                       >  2:2718262/1‑150 (MQ=255)
 cagcGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTc                                                                                                                                                      >  2:2086316/1‑150 (MQ=255)
 cagcGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTc                                                                                                                                                      <  2:3596945/150‑1 (MQ=255)
  agcGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCa                                                                                                                                                     >  2:2900607/1‑150 (MQ=255)
      gCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGc                                                                                                                                                 >  2:3695657/1‑150 (MQ=255)
        cAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGAGTCTGCAAAGCGTTGTTCAACGCtt                                                                                                                                               >  1:1456219/1‑150 (MQ=255)
         aGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGAACAACGCttt                                                                                                                                              <  2:154388/150‑2 (MQ=255)
         aGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGAACAACGCttt                                                                                                                                              <  2:4939015/150‑2 (MQ=255)
         aGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGAACAACGCttt                                                                                                                                              <  1:3262495/150‑2 (MQ=255)
          gAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTgg                                                                                                                                             <  2:2797586/150‑1 (MQ=255)
          gAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGAACAACGCtttg                                                                                                                                             <  1:2554637/150‑3 (MQ=255)
           caGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTggg                                                                                                                                            <  1:787922/149‑1 (MQ=255)
           aaGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTggg                                                                                                                                            <  1:765723/150‑1 (MQ=255)
            aGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCCAAGCGTTGAACAACGCtttgca                                                                                                                                           <  1:1735785/150‑5 (MQ=255)
             gccgccAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGc                                                                                                                                          <  1:4740480/150‑1 (MQ=255)
              ccgccAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGcc                                                                                                                                         >  1:3409796/1‑150 (MQ=255)
               cgccAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCt                                                                                                                                        <  2:1392164/150‑1 (MQ=255)
                gccAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTg                                                                                                                                       >  1:4514399/1‑150 (MQ=255)
                gccAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTg                                                                                                                                       >  1:4092409/1‑150 (MQ=255)
                  cAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCt                                                                                                                                     >  2:1354023/1‑150 (MQ=255)
                   aaGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCtc                                                                                                                                    <  1:394135/150‑1 (MQ=255)
                      gCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCtcatt                                                                                                                                 >  1:3424730/1‑150 (MQ=255)
                          gCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTcattcatt                                                                                                                             <  1:4230709/150‑1 (MQ=255)
                           cGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCattcattc                                                                                                                            <  1:1094274/150‑1 (MQ=255)
                            gTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTcc                                                                                                                           >  2:4405335/1‑150 (MQ=255)
                            gTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTcc                                                                                                                           >  1:1638083/1‑150 (MQ=255)
                             ttGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCg                                                                                                                          <  2:842906/150‑1 (MQ=255)
                              tGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCgg                                                                                                                         <  2:5024536/150‑1 (MQ=255)
                              tGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCgg                                                                                                                         <  1:5019940/150‑1 (MQ=255)
                               gATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCggg                                                                                                                        <  1:4206794/150‑1 (MQ=255)
                                aTCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGt                                                                                                                       >  2:862651/1‑150 (MQ=255)
                                 tCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTg                                                                                                                      >  1:2639003/1‑150 (MQ=255)
                                 tCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTg                                                                                                                      >  1:5035556/1‑150 (MQ=255)
                                 tCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGACGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTg                                                                                                                      >  2:761598/1‑150 (MQ=255)
                                  cGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGa                                                                                                                     <  1:3832403/150‑1 (MQ=255)
                                   ggCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGaa                                                                                                                    <  1:2387434/150‑1 (MQ=255)
                                    gCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAg                                                                                                                   <  1:862651/150‑1 (MQ=255)
                                      aGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCt                                                                                                                 >  1:3225911/1‑150 (MQ=255)
                                       ggAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTc                                                                                                                <  2:4277722/150‑1 (MQ=255)
                                        gAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTcc                                                                                                               >  2:203969/1‑150 (MQ=255)
                                            aGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTc                                                                                                           <  1:4236331/150‑1 (MQ=255)
                                               gtCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCa                                                                                                        <  2:1034089/150‑1 (MQ=255)
                                                ttcctgcgcCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCaa                                                                                                       >  1:3035124/6‑150 (MQ=255)
                                                tCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCaa                                                                                                       <  1:4045651/150‑1 (MQ=255)
                                                tCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCaa                                                                                                       >  1:2238922/1‑150 (MQ=255)
                                                  ccGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCGCAGGTGCCGAGGGCGGCGGGGGTCGTTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGa                                                                                                     >  2:4463435/1‑150 (MQ=255)
                                                      cgcCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGc                                                                                                 <  1:268033/150‑1 (MQ=255)
                                                        ccGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCaa                                                                                               <  2:1842139/150‑1 (MQ=255)
                                                         cGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGGCGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAAc                                                                                              <  1:460035/150‑1 (MQ=255)
                                                                tgtgCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTga                                                                                       >  2:739945/1‑150 (MQ=255)
                                                                tgtgCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTga                                                                                       <  2:1783235/150‑1 (MQ=255)
                                                                   gCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACACGCAAGACTGCAACAGGCTGAGaa                                                                                    >  1:2843649/1‑150 (MQ=255)
                                                                      gAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATgg                                                                                 <  1:5019010/150‑1 (MQ=255)
                                                                      gAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATgg                                                                                 <  1:2989344/150‑1 (MQ=255)
                                                                            gCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAAt                                                                           <  1:557422/150‑1 (MQ=255)
                                                                             cccAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATa                                                                          <  2:3642578/150‑1 (MQ=255)
                                                                               cAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGa                                                                        >  1:489577/1‑150 (MQ=255)
                                                                               cAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTAATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGa                                                                        >  1:3468019/1‑150 (MQ=255)
                                                                                      cAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCg                                                                 >  2:3429317/1‑150 (MQ=255)
                                                                                        ggtggtGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGcc                                                               >  2:150697/1‑150 (MQ=255)
                                                                                         gtggtgTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCaa                                                                                                       >  1:2574527/1‑109 (MQ=255)
                                                                                         gtggtgTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCaa                                                                                                       <  2:2574527/109‑1 (MQ=255)
                                                                                            gtgtCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTTCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAAc                                                           <  1:4317162/150‑1 (MQ=255)
                                                                                                  tGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAAcc                                                     <  1:2193424/150‑1 (MQ=255)
                                                                                                    aCGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGACCAACTCAACCgg                                                   >  2:1130203/1‑150 (MQ=255)
                                                                                                          ccAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAaga                                             <  1:4836322/150‑1 (MQ=255)
                                                                                                          ccAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAaga                                             <  2:4593097/150‑1 (MQ=255)
                                                                                                          ccAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAaga                                             <  1:4413723/150‑1 (MQ=255)
                                                                                                             ggTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAggg                                          <  2:3594813/150‑1 (MQ=255)
                                                                                                             ggTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAggg                                          <  2:1873231/150‑1 (MQ=255)
                                                                                                             ggTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAggg                                          <  1:1928832/150‑1 (MQ=255)
                                                                                                               tGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGgc                                        <  2:379321/150‑1 (MQ=255)
                                                                                                                  cGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGcgcg                                     >  1:2084234/1‑150 (MQ=255)
                                                                                                                  cGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGcgcg                                     >  1:1693971/1‑150 (MQ=255)
                                                                                                                         ggcggGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGc                              <  2:4449752/150‑1 (MQ=255)
                                                                                                                          gcggGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCa                             <  1:4651462/150‑1 (MQ=255)
                                                                                                                              gggTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATCGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGa                         >  2:4483626/1‑150 (MQ=255)
                                                                                                                              gggTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGa                         >  1:3815097/1‑150 (MQ=255)
                                                                                                                               ggTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGaa                        <  2:4514399/150‑1 (MQ=255)
                                                                                                                                gTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGaaa                       >  2:5007570/1‑150 (MQ=255)
                                                                                                                                    aTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGtt                   <  1:1340364/150‑1 (MQ=255)
                                                                                                                                    aTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGtt                   <  2:3419259/150‑1 (MQ=255)
                                                                                                                                    aTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGtt                   <  1:361581/150‑1 (MQ=255)
                                                                                                                                    aTTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGtt                   <  2:3004739/150‑1 (MQ=255)
                                                                                                                                     ttCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTg                  >  1:2504863/1‑150 (MQ=255)
                                                                                                                                      tCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGc                 <  2:4191898/150‑1 (MQ=255)
                                                                                                                                      tCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATCGTGCACGGAAAGGTTGc                 >  1:3516117/1‑150 (MQ=255)
                                                                                                                                       cTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCa                <  2:1507520/150‑1 (MQ=255)
                                                                                                                                       cTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCa                <  1:3679542/150‑1 (MQ=255)
                                                                                                                                          ccaaGCGTTGAACAACGCTTTGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACtt             >  2:1436287/3‑150 (MQ=255)
                                                                                                                                          ccaaGCGTTGAACAACGCTTAGGGACTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGCTTGCACtt             >  2:4594772/3‑150 (MQ=255)
                                                                                                                                          cAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACtt             <  2:4092409/150‑1 (MQ=255)
                                                                                                                                           aaaGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTg            >  2:516034/1‑150 (MQ=255)
                                                                                                                                            aaGCGTTGAACAACGCTTTGGGACTGCACATTAATTCAGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACTTTTCGCACAACTCAACAGGAAAAGATGGGCGCGATTGTGCACGGACAGGATGCACTTGc           >  2:3023024/1‑150 (MQ=255)
                                                                                                                                             aGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCg          <  1:4405335/150‑1 (MQ=255)
                                                                                                                                              gCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGa         <  2:1033896/150‑1 (MQ=255)
                                                                                                                                                   gttCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGccc    <  2:4149944/150‑1 (MQ=255)
                                                                                                                                                    ttCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCt   >  1:510670/1‑150 (MQ=255)
                                                                                                                                                    ttCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCt   >  1:1868647/1‑150 (MQ=255)
                                                                                                                                                     tCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCtt  >  2:448088/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCGGCCAGAAGCCGCCAAGGCCGGCGTTGATCGGCAGGAACAGTGTCCCGGCGCCGATGGCTGTGCCGAACAGGCCCAGCATCCAGGTGGTGTCGTGACGCGCCCAGGTGCCGAGGGCGGCGGGGGTCGATTCTGCAAAGCGTTGTTCAACGCTTGGGGCCTGCTCATTCATTCCGGGTGAAGCTCCACTCGCAAGACTGCAACAGGCTGAGAATGGCGAAATAGACGTTTCGCCCAACTCAACCGGAAAAGAGGGGCGCGATTGTGCACGGAAAGGTTGCACTTGCGAAGCCCTT  >  NC_002947/4080988‑4081286

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