Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,422,998 C→G 5.2% N51K (AAC→AAG)  ulaE → L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,9980CG 5.2% 415.2 / 10.5 134N51K (AAC→AAG) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (72/55);  new base G (5/2);  total (77/57)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 135 total aligned reads displayed.

TTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTA  >  NC_000913/4422849‑4423146
                                                                                                                                                      |                                                                                                                                                    
ttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAAc                                                                                                                                                      >  1:1012551/1‑150 (MQ=255)
  gtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                                    >  1:1722768/1‑150 (MQ=255)
  gtCCAAACAAATCCCGCTTGGCATATATGATTAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGGCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                                    <  2:695894/150‑1 (MQ=255)
     caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                                 >  1:1372656/1‑150 (MQ=255)
     caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTAGATATTGTCGAAATGTCGGTAGAGG‑AAACTGACGATCGCCTGTCGCGCCTCTACTGGAg                                                                                                                                                 >  1:1051841/1‑150 (MQ=255)
      aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGc                                                                                                                                                >  2:1666671/1‑150 (MQ=255)
       aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCTAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               >  1:658253/1‑150 (MQ=255)
       aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               <  1:368264/150‑1 (MQ=255)
       aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               >  1:43735/1‑150 (MQ=255)
          aaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCcgcg                                                                                                                                            <  2:1058347/150‑1 (MQ=255)
            aTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGGCTCCAGTTGAGGCGCGac                                                                                                                                          <  1:905777/150‑2 (MQ=255)
             tCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGGCTCCAGTTGAGGCGCGaca                                                                                                                                         <  1:444849/150‑3 (MQ=255)
               cctcgtGGCATCTCTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAAGGTCGGTAGATC‑ATACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGagcag                                                                                                                                       <  1:1576156/145‑1 (MQ=255)
                   ttGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                                                                                                                                   <  2:1017003/150‑1 (MQ=255)
                      gCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTgg                                                                                                                                >  2:1481295/1‑150 (MQ=255)
                       cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgc                                                                                                                               <  1:993562/150‑1 (MQ=255)
                       cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgc                                                                                                                               >  2:1542890/1‑150 (MQ=255)
                         tCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgcgc                                                                                                                             >  2:1595546/1‑150 (MQ=255)
                            aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCggg                                                                                                                          >  1:1512597/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑ACACTGACGATCGCCTGTCGCTCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTc                                                                                                                        >  2:158767/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCACCTGGCGCCGGGAGCAGCGGCTGGCGCTGGTc                                                                                                                        >  2:164183/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTc                                                                                                                        >  1:619571/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATGTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTc                                                                                                                        >  2:100943/1‑150 (MQ=255)
                                aaaaGCGCTCCCCGACGGGGAGTGCTGGCTGGAACGCCTGCAAATTGCAAAAACGTTAGGCTTCTAATTTGTCGAAATGTCTGTATATG‑AAAATTACGATAGCCTGTCGAGCCTCAACTGGAGCCGCGACCATAGTCTGGCGCTTGTCaa                                                                                                                      >  2:611029/1‑150 (MQ=255)
                                 ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAAt                                                                                                                     >  2:1628248/3‑150 (MQ=255)
                                 ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGAAG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCGGGTCAAt                                                                                                                     >  2:1270902/3‑150 (MQ=255)
                                   aGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc                                                                                                                   <  2:1494394/150‑1 (MQ=255)
                                   aGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCGTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTTTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc                                                                                                                   <  2:1395684/150‑1 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGa                                                                                                                 >  1:1299726/1‑150 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                >  2:1250512/1‑150 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGAACGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                <  2:6467/150‑1 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGACTGCAACTGGCAAAAACGTTAGGCATCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGAGCTGGTCAATGCGAt                                                                                                                >  1:1121640/1‑150 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACCCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                >  2:708955/1‑150 (MQ=255)
                                       cTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAtt                                                                                                               >  1:1265707/1‑150 (MQ=255)
                                         ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgt                                                                                                             >  2:317225/1‑150 (MQ=255)
                                         ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgt                                                                                                             >  1:2049/1‑150 (MQ=255)
                                           ccgccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCATAAACGTTAGGCTGCGATTTTGTCGAATTGTCGGTAGATG‑ATACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgttg                                                                                                           <  1:1247782/150‑1 (MQ=255)
                                            cgccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCATTGCGATTGTTGa                                                                                                          >  1:1287336/1‑150 (MQ=255)
                                                ggggAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAAcc                                                                                                      >  1:1426376/1‑150 (MQ=255)
                                                     gTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGgcgt                                                                                                 <  2:155173/150‑1 (MQ=255)
                                                        ctggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCTATTGTTGAAACCGgcgtgcg                                                                                              >  2:121867/1‑150 (MQ=255)
                                                        ctggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGgcgtgcg                                                                                              >  2:496623/1‑150 (MQ=255)
                                                                     cTGCAACTGGCAAAATCGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGGCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtg                                                                                 <  1:714510/150‑1 (MQ=255)
                                                                     cTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtg                                                                                 <  2:916364/150‑1 (MQ=255)
                                                                       gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                               >  1:529320/1‑150 (MQ=255)
                                                                       gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                               <  2:1012551/150‑1 (MQ=255)
                                                                        cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                              <  2:1379624/150‑1 (MQ=255)
                                                                        cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                              <  1:317225/150‑1 (MQ=255)
                                                                         aaCTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGcc                                                                             <  1:1416014/150‑1 (MQ=255)
                                                                             ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                         >  2:314375/1‑150 (MQ=255)
                                                                             ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                         >  1:1055927/1‑150 (MQ=255)
                                                                               cAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGAGAGCAGCGTCAGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTg                                                                       >  2:203323/1‑150 (MQ=255)
                                                                                aaaaaCGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGc                                                                      >  2:818607/1‑150 (MQ=255)
                                                                                      gTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcg                                                                <  2:619571/150‑1 (MQ=255)
                                                                                        tAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                              >  2:229288/1‑150 (MQ=255)
                                                                                        tAGGCTGCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                              <  1:881141/150‑1 (MQ=255)
                                                                                          ggCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                             <  1:1413578/149‑1 (MQ=255)
                                                                                          ggCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                             >  2:1413578/1‑149 (MQ=255)
                                                                                            cTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGttt                                                          >  1:1044867/1‑150 (MQ=255)
                                                                                              tCGACTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGAGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTcc                                                        <  1:895252/150‑1 (MQ=255)
                                                                                                 aTTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGTTCAATGCGATTGTTGAACCCGGCGTGCGCGAGACGTCCATGTGCCTTTCTGCGAATCGTCGTATccccc                                                     >  2:403618/1‑148 (MQ=255)
                                                                                                 aTTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGc                                                     <  1:227385/150‑1 (MQ=255)
                                                                                                     tgacGAACTGTCGGTCGATG‑AAACTGACGATCGACTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAACCCGGCGTGCGCGTGCCGTCCATGAGCCTTTCTGCTCATCGTCGTTTCCCGCTggg                                                 <  1:101637/147‑1 (MQ=255)
                                                                                                     tGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTggg                                                 <  2:2049/150‑1 (MQ=255)
                                                                                                      gTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGc                                                >  1:1136465/1‑150 (MQ=255)
                                                                                                       tCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCa                                               >  1:381469/1‑150 (MQ=255)
                                                                                                       tCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCa                                               >  1:549402/1‑150 (MQ=255)
                                                                                                        cGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCCGCGTCTGGCGCGGGTCACTTCGAGTGTTGAAACCGGCGGGAGCGTGCCGGCCATGTGCCGTTCTGCTCCTCGTCGTGTCCCGCTGGGCAg                                              >  2:286148/1‑150 (MQ=255)
                                                                                                         gAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGGTGGGctgg                                             >  1:1160707/1‑147 (MQ=255)
                                                                                                         gAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGt                                             >  2:1270894/1‑150 (MQ=255)
                                                                                                             tGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAg                                         <  2:1549088/150‑1 (MQ=255)
                                                                                                              gTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGGGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGa                                        >  1:1019805/1‑150 (MQ=255)
                                                                                                                cGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATg                                      >  2:1132413/1‑150 (MQ=255)
                                                                                                                cGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATg                                      >  1:1723978/1‑150 (MQ=255)
                                                                                                                 ggTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCTCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGa                                     >  2:1452290/1‑150 (MQ=255)
                                                                                                                  gaaGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGAc                                    <  2:934554/148‑1 (MQ=255)
                                                                                                                        g‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTg                              <  2:51928/150‑1 (MQ=255)
                                                                                                                             cTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCggg                          <  1:1481295/150‑1 (MQ=255)
                                                                                                                             cTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCggg                          >  1:423881/1‑150 (MQ=255)
                                                                                                                              tGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGgc                         <  1:252559/150‑1 (MQ=255)
                                                                                                                               gACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGgcg                        >  2:406147/1‑150 (MQ=255)
                                                                                                                                 cGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCa                      >  2:1596069/1‑150 (MQ=255)
                                                                                                                                  gATCGCCTGTCGCGCCGCAACTGGAGCCGCGAGCTGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATTTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAg                     <  2:1051841/150‑1 (MQ=255)
                                                                                                                                   aTCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAgg                    <  2:381469/150‑1 (MQ=255)
                                                                                                                                    tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGATGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg                   <  1:979712/150‑1 (MQ=255)
                                                                                                                                      gCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCGGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGc                 >  1:1452610/1‑150 (MQ=255)
                                                                                                                                         tGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTgg              <  2:1055927/150‑1 (MQ=255)
                                                                                                                                         tGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTgg              <  2:43735/150‑1 (MQ=255)
                                                                                                                                         gctCGCGGCTCCAGTTGAGGCGCGATCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGTCGTGCGCGTGCAGTCCATGTGCCTTTCTGCTCATAGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCTGGCGCAGGGGCTgg              >  2:241555/3‑150 (MQ=255)
                                                                                                                                          ctCGCGGCTCCAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCGTCGTCGTTTCCCGCTGGGCAGTGAAGATGACCCGGTGCGGGCGCAGGGGCgggg             >  1:170657/2‑149 (MQ=255)
                                                                                                                                          ctCGCGGCTCCAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACAGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTGCCGCTGGGCAGGGAAGATGACGAGGTGCGGGCGCAGGGGCTggc             >  1:875919/2‑149 (MQ=255)
                                                                                                                                           tCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGgag            <  1:502928/150‑1 (MQ=255)
                                                                                                                                                 cTCAACTGGAGCCGCGAGCAGCGTCTGGCGGTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTGATCGTCGTTTCCCGCAGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCAGGAGATTATg      >  1:778245/1‑150 (MQ=255)
                                                                                                                                                 cTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATg      <  2:1722768/150‑1 (MQ=255)
                                                                                                                                                  gcAACTGGAGCCGCGAGCAGCGGCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGc     <  1:1562563/149‑1 (MQ=255)
                                                                                                                                                   cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGTCGTGCGCGTGCCGTCCATGTGCCTTTATGCTCATCGTCGTTTACCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCg    >  2:398134/1‑150 (MQ=255)
                                                                                                                                                   cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATAGTCGTTTACCGCTGGGCAGTGAAGATGAAGAGGTGCGTGCGCAGGGGCTGGAGATTATGCg    >  2:391938/1‑150 (MQ=255)
                                                                                                                                                    aaCTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGt   <  1:1250512/150‑1 (MQ=255)
                                                                                                                                                    aaCGGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGGGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGGCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGt   <  1:708955/150‑1 (MQ=255)
                                                                                                                                                     aCTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTa  <  2:484221/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTA  >  NC_000913/4422849‑4423146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑