Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,403 A→G 5.2% intergenic (‑241/+143) sdaC ← / ← amaC serine:H+ symport permease/D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,4030AG 5.2% 709.3 / 3.6 239intergenic (‑241/+143)sdaC/amaCserine:H+ symport permease/D‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (100/113);  new base G (11/2);  total (123/116)
Fisher's exact test for biased strand distribution p-value = 9.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 243 total aligned reads displayed.

TGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTC  >  NC_002947/4081254‑4081552
                                                                                                                                                     |                                                                                                                                                      
tgtgCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGA                                                                                                                                                        <  1:1197012‑M1/150‑1 (MQ=255)
     aCGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGCTAggg                                                                                                                                                   >  2:2348777/1‑150 (MQ=255)
      cGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGc                                                                                                                                                  <  1:8242868/150‑1 (MQ=255)
         aaaGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGc                                                                                                                                               <  1:6921303/150‑1 (MQ=255)
                 cACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACACGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGc                                                                                                                                       <  1:3938838/150‑1 (MQ=255)
                   cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGTGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                                     <  1:1424003/150‑1 (MQ=255)
                   cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccc                                                                                                                                     <  1:7996863/150‑1 (MQ=255)
                   cTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGccc                                                                                                                                     <  2:7362301/150‑1 (MQ=255)
                    ttGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccct                                                                                                                                    <  2:8455644/150‑2 (MQ=255)
                       cGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTTTGCGGccctatc                                                                                                                                 <  1:8766187/150‑5 (MQ=255)
                         aaGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGGGGCGCGGGTGGGGCGCCAAGCGGCCCCAGcca                                                                                                                               >  1:8080262/1‑148 (MQ=255)
                          aGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAt                                                                                                                              <  1:3814357/150‑1 (MQ=255)
                              cTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGCGGCGCGATGGGGGCGCAGAGCCGCGCCAGCAAATTcc                                                                                                                          >  1:2080596/1‑149 (MQ=255)
                              cTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCa                                                                                                                          <  2:7478606/150‑1 (MQ=255)
                               ttttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCCGTAGGGCCGCCAAGCCGCCCCACCAATATCAg                                                                                                                         >  2:8840474/1‑150 (MQ=255)
                               ttttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCCCGGTAGGGCCGCAAAGCGGCCCCCCCAAATTCAg                                                                                                                         >  2:1573330/1‑150 (MQ=255)
                                tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                                        <  2:7060173/150‑1 (MQ=255)
                                tttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGc                                                                                                                        <  2:7276812/150‑1 (MQ=255)
                                 ttCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCg                                                                                                                       <  2:5413877/150‑1 (MQ=255)
                                    cGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCa                                                                                                                    >  1:4187216/1‑150 (MQ=255)
                                      aggagCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCaaa                                                                                                                  <  1:2602820/150‑1 (MQ=255)
                                       gcagCGGTTGCACTAGTCTTTACAAGAGCAAGGCATTGCGGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTATGGGGCGGCCTGTGGGAGCGGCCAATCGTCGCGAAAGGGCCCCAAAGCGGCCCCAGCAAATTCAGCGGCAAAg                                                                                                                 >  2:4406129/3‑150 (MQ=255)
                                         agCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCt                                                                                                               <  2:2845928/150‑1 (MQ=255)
                                          gCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTg                                                                                                              <  1:6301440/150‑1 (MQ=255)
                                           cGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGc                                                                                                            <  2:3226899/150‑1 (MQ=255)
                                             gTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCaa                                                                                                          >  1:8890186/1‑150 (MQ=255)
                                                 cACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGGCGCAAAGCGGCCCCAGCAATTTAAGCGGCACAGCTGCCATCGt                                                                                                      >  2:5479500/1‑150 (MQ=39)
                                                  aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                                     <  1:4152675/150‑1 (MQ=255)
                                                  aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                                     <  2:4187216/150‑1 (MQ=255)
                                                  aCTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTa                                                                                                     <  2:7607988/150‑1 (MQ=255)
                                                   cTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCCCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCAGCAAACGTAg                                                                                                    >  2:4370574/1‑150 (MQ=255)
                                                    ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTTCGTCGCGATAGGGCCCCAAAGCGGCCCCACCAATTACACCCGCAAAGCGGCAAACGTAgg                                                                                                   >  2:583078/1‑150 (MQ=255)
                                                    ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAgg                                                                                                   <  1:450392/150‑1 (MQ=255)
                                                    ttGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCGCAAATTCAGCGGCAAAGCGGCAAAAGTAgg                                                                                                   >  1:7594764/1‑150 (MQ=38)
                                                         tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAAATTCAGCGGCAAAGCTGCAAACGTAGGGgcgc                                                                                              >  2:8972103/1‑150 (MQ=255)
                                                         tttACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTCGTCGCGATATGGCCGCAAAGCGGCCCCACCAAATTAAGCGGCAAAACTGCAATAGGATGGgcgc                                                                                              >  2:3320720/1‑150 (MQ=255)
                                                            aCAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTcc                                                                                           <  2:5555087/150‑1 (MQ=255)
                                                              aaGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTc                                                                                                                           >  1:3639703/1‑117 (MQ=255)
                                                              aaGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTc                                                                                                                           <  2:3639703/117‑1 (MQ=255)
                                                               aGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCa                                                                                        <  2:4864030/150‑1 (MQ=255)
                                                                gCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCCGCGCGCTAGGGCCGCCAACCGGCCCCCCCCAATTAAGCGGCCAAGCTGCGCACGTAGGGGCGcgccgccc                                                                                       >  2:8479780/1‑148 (MQ=255)
                                                                  ggAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCAccc                                                                                     <  1:3536783/150‑1 (MQ=255)
                                                                   gAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTCGTCGCGATAGGGCCGCAAAACGGCCCCACCAAATTAAGCGGCAAAACTGCAATAGTAGGGGCGCGCCGCAccca                                                                                    >  2:3442109/1‑149 (MQ=38)
                                                                   gAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGTCGGGATTGGGCCGCAAAGCGGCCCCCGCAATTACAGCCGCAAAGCTGCAATCGAAGGGGCGCTCCGCAcccg                                                                                    >  1:5285633/1‑149 (MQ=38)
                                                                      ggCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCAAGCGTCCCCATAATGCCCCAAAACGCCCCCAGCAAATTAAGCGGCAAAACAGAAAAAGTACGGGCGCACCCCCCCCattt                                                                                 >  2:1180419/1‑150 (MQ=255)
                                                                         gtaactcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTcgc                                                                              >  1:5027794‑M1/54‑150 (MQ=255)
                                                                         aTTTCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgc                                                                              >  2:6626189‑M2/1‑146 (MQ=255)
                                                                         aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCTTGCGTGGCGATAGGGCCGCAAATCGGCCCCAGCAATTTAAGCGGCAAATCTGCAATAGTTGGGGGGCTCCGCACCCATTtagc                                                                              <  2:3772292‑M2/150‑5 (MQ=255)
                                                                         aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgc                                                                               <  2:6082919‑M1/150‑74 (MQ=255)
                                                                         aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCa                                                                                                                          <  1:4389114/107‑1 (MQ=255)
                                                                         aTTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgc                                                                              >  1:724769‑M2/1‑146 (MQ=255)
                                                                           tGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcga                                                                            >  2:2640343‑M2/1‑144 (MQ=255)
                                                                            gCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTCtcgcgac                                                                           >  2:2253155‑M2/1‑143 (MQ=255)
                                                                            actcgccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGAc                                                                           >  1:2706324‑M1/51‑150 (MQ=255)
                                                                              tGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGGCGTGGGGGGCGGGAGGGCCGCAAACCGGGCCCCCCAAATTCACCGGCAAACCAGCAAAAGGATGGGCGCGCCGCCCCCGTTTCGcgccgc                                                                         >  1:2883303/1‑150 (MQ=255)
                                                                              tGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgc                                                                         >  1:1828228‑M2/1‑141 (MQ=255)
                                                                               ggTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCCATCGTGGGGGCGCTCCGCACCCATTtcgcgtcgca                                                                        <  2:8213483‑M2/150‑11 (MQ=255)
                                                                               ggTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCGCAAACCGGCCCCACCAATTTGAGCGGCAAAGCAGCAATCGGAGGGGCGCTCCGCACCCGTTtcgcgacgcc                                                                        <  2:2274024‑M2/150‑11 (MQ=255)
                                                                                gtgtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaa                                                                       >  1:2883869‑M2/1‑139 (MQ=255)
                                                                                 ccaacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAg                                                                      >  1:469682‑M1/46‑150 (MQ=255)
                                                                                  gtCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaagg                                                                     >  1:2740188‑M2/1‑137 (MQ=255)
                                                                                  caacgggccaactgcggcccctcctgtggggcgcggttggacttTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAgg                                                                     >  1:8975962‑M1/45‑150 (MQ=255)
                                                                                   tCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaaggc                                                                    >  2:418931‑M2/1‑136 (MQ=255)
                                                                                    cGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaaggcc                                                                   >  1:2909127‑M2/1‑135 (MQ=255)
                                                                                     gATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCGTTtcgcgacgcaaggccg                                                                  <  2:8135124‑M2/150‑17 (MQ=255)
                                                                                      aTCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtc                                                                                >  2:2641899‑M2/1‑133 (MQ=255)
                                                                                      aTCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCCCAAAGCGGCCCCAGCAAATTAAGCGGCAAACCTGCAAAAGTAGGGGCGCTCCGCCCCCGTTtagcgacgcaaggccgc                                                                 <  2:4462068‑M2/150‑18 (MQ=255)
                                                                                        cGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTTGGTCGCCATAAGGGCCCAAAACCGCCCCACCAAATTAAGCCGCAAAGCCGCCAACCGGAGGGCGCGCCCCACCCCGTTCGCGACCCAAcgccgcga                                                               >  2:77063/1‑148 (MQ=255)
                                                                                          gTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGCTAGGGCCGCAAAGCGGCCCCCCCAAATTCAGCGGCCAAGCTGCAATCGTTGGGGCGCTACGCCCCCATTtcgcgacgcaaggccgctccc                                                             <  2:5074453‑M2/150‑22 (MQ=255)
                                                                                               aTGTTGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                                        <  2:5323724/150‑1 (MQ=255)
                                                                                               aTGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACaaa                                                        <  2:6809019/150‑1 (MQ=255)
                                                                                                          tGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCGTGCCGCGCCCTAGGGCCCCAAACCGGCCCCAGCCAATTAAGCCGCCAAACGGCAATCGGAGGGGCGCGCCGCCCCCGTTTAGCGACCCAAGGCCGCTCCCACAAACTACAAccccc                                             >  2:1938835/1‑147 (MQ=255)
                                                                                                              cgcCTGTGGGGCGGCCTGTGGGAGCGGCCATGCGTCGCGATAGGGCCGCAAAGCCGCCCCAGCAAATTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACCAAGTCCAACCGcccccc                                         >  2:6086633/1‑150 (MQ=37)
                                                                                                                         cggcCGGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGAAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTCCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGc                              <  2:6263639/150‑1 (MQ=255)
                                                                                                                                tgGGAGCGGCCTTGCGTCGCTAaatgggggcggcgcgcccctccgtttgcagctttgccgcttaaattggtggggccgttttggggccctatcgcgacgcacggccgctcccacaggccgccccacaggcgcccaggccgctacattcca                        <  2:4749269‑M1/150‑129 (MQ=255)
                                                                                                                                TGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAACGCTGCAATCGTTGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGtt                       >  2:6466626‑M1/1‑150 (MQ=255)
                                                                                                                                 gggCGCGGCGTTGCGTGGCGAactgggtgcggagcgcccctacgtttgccgctttgcggctgaatttgcgggggccgctttgcggccatatcgcgccgcacggccgctcccacaggccgccccacaggcgcccaggccgctacattccac                       <  2:4740136‑M1/150‑130 (MQ=255)
                                                                                                                                 gggCGCGGCCTTGCGTCGCTAaatgggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctaccgcgaggcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccac                       <  2:8797806‑M1/150‑130 (MQ=255)
                                                                                                                                 gggCGCGGCCTTGCGTCGCGAaatgggtgcggagcgcccctacgattgcagctttgccgcttaaattgctggggccgctttgcggccctaccgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccac                       <  2:2186871‑M1/150‑130 (MQ=255)
                                                                                                                                   GAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCACTCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACCGGAGGGGCCGCAGTTGGc                    <  2:5758591‑M1/150‑1 (MQ=255)
                                                                                                                                    aGCGGCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccga                    >  2:2334961‑M1/1‑18 (MQ=255)
                                                                                                                                    aGCGGCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccga                    >  1:6032526‑M1/1‑18 (MQ=255)
                                                                                                                                        gCCTTGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcga                >  2:5336156‑M1/1‑14 (MQ=255)
                                                                                                                                          ctggCGTCGCGATAGGGCCGCAACGCGGCCCCCGCAATTTCCGCGGCACAGCTGCAATCGTCGGGGCGCTCCGCACCCTTTTCGCGCCGCACGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg             <  2:7039745/147‑1 (MQ=38)
                                                                                                                                          CTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTgg             <  1:5484406‑M1/150‑1 (MQ=255)
                                                                                                                                           ttGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacac             >  1:3328000‑M1/1‑11 (MQ=255)
                                                                                                                                           ttGCGTCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacac             >  1:3394206‑M1/1‑11 (MQ=255)
                                                                                                                                           TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAACGCTGCAATCGTAGGGGCGCTCCACACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTGGGCCCGTTGGc            <  1:6115994‑M1/150‑1 (MQ=255)
                                                                                                                                           TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc            <  1:8468469‑M1/150‑1 (MQ=255)
                                                                                                                                           TTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGc            <  1:3865304‑M1/150‑1 (MQ=255)
                                                                                                                                               ggcgcgaactgggtgcggcgcgcccctacgttggcagctttgcggctgaatttgctggggcggctttgcggccctatcgcgacgcacggccgctcccacaggccgccccacaggcgcccatgccgctacattccaccgatcgacaccagc         <  2:611104‑M1/148‑144 (MQ=255)
                                                                                                                                               gTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTtcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacaccagc        <  2:8216180‑M2/150‑75 (MQ=255)
                                                                                                                                                tCGCGAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacaccagca        >  2:2640339‑M1/1‑6 (MQ=255)
                                                                                                                                                 CGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTa      <  1:726413‑M1/150‑1 (MQ=255)
                                                                                                                                                  gcgAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgccccggccgctacattccaccgatcgaccacagcaaa      >  1:8145897‑M1/1‑4 (MQ=255)
                                                                                                                                                  gcgAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacaccagcaat      >  1:724258‑M1/1‑4 (MQ=255)
                                                                                                                                                   cgAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacaccagcaatg     >  1:5842068‑M1/1‑3 (MQ=255)
                                                                                                                                                   CGATCGGGCCGCAACGCGCCCCCAGCAATTTCAGCGGCATAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACt    >  2:3049895‑M1/1‑150 (MQ=255)
                                                                                                                                                   CGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACt    <  2:1362406‑M1/150‑1 (MQ=255)
                                                                                                                                                    gAaaagggtgcggagcgcccctacgattgcagctttgccgctgaaattgctggggccgctttgcggccctatcgcgacgcaaggccgctcccacaggccgccccacaggcgcccaggccgctacattccaccgatcgacaccagcaatgc    >  2:8890343‑M1/1‑2 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGTGCACGGAAAGGTTGCACTTGCGAAGCCCTTTTCCGAGGAGCGGTTGCACTTGTCTTTACAAGCGGAAGGCATTGCTGGTGTCGATCGGTGGAATGTAGCGGCCTGGGCGCCTGTGGGGCGGCCTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCAGCAATTTCAGCGGCAAAGCTGCAATCGTAGGGGCGCTCCGCACCCTTTTCGCGACGCAAGGCCGCTCCCACAAAGTCCAACCGCGCCCCACAGGAGGGGCCGCAGTTGGCCCGTTGGCGAGTTACTTC  >  NC_002947/4081254‑4081552

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