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→AAGulaEL‑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.

TTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTA  >  NC_000913/4422849‑4423146
                                                                                                                                                      |                                                                                                                                                    
ttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAAc                                                                                                                                                      >  1:1012551/1‑150 (MQ=255)
  gtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                                    >  1:1722768/1‑150 (MQ=255)
  gtCCAAACAAATCCCGCTTGGCATATATGATTAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGGCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                                    <  2:695894/150‑1 (MQ=255)
     caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                                 >  1:1372656/1‑150 (MQ=255)
     caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTAGATATTGTCGAAATGTCGGTAGAGGAAACTGACGATCGCCTGTCGCGCCTCTACTGGAg                                                                                                                                                 >  1:1051841/1‑150 (MQ=255)
      aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGc                                                                                                                                                >  2:1666671/1‑150 (MQ=255)
       aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCTAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               >  1:658253/1‑150 (MQ=255)
       aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               <  1:368264/150‑1 (MQ=255)
       aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               >  1:43735/1‑150 (MQ=255)
          aaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCcgcg                                                                                                                                            <  2:1058347/150‑1 (MQ=255)
            aTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGGCTCCAGTTGAGGCGCGac                                                                                                                                          <  1:905777/150‑2 (MQ=255)
             tCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGGCTCCAGTTGAGGCGCGaca                                                                                                                                         <  1:444849/150‑3 (MQ=255)
               cctcgtGGCATCTCTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAAGGTCGGTAGATCATACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGagcag                                                                                                                                       <  1:1576156/145‑1 (MQ=255)
                   ttGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                                                                                                                                   <  2:1017003/150‑1 (MQ=255)
                      gCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTgg                                                                                                                                >  2:1481295/1‑150 (MQ=255)
                       cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgc                                                                                                                               <  1:993562/150‑1 (MQ=255)
                       cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgc                                                                                                                               >  2:1542890/1‑150 (MQ=255)
                         tCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgcgc                                                                                                                             >  2:1595546/1‑150 (MQ=255)
                            aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCggg                                                                                                                          >  1:1512597/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGACACTGACGATCGCCTGTCGCTCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTc                                                                                                                        >  2:158767/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCACCTGGCGCCGGGAGCAGCGGCTGGCGCTGGTc                                                                                                                        >  2:164183/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTc                                                                                                                        >  1:619571/1‑150 (MQ=255)
                              gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATGTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTc                                                                                                                        >  2:100943/1‑150 (MQ=255)
                                aaaaGCGCTCCCCGACGGGGAGTGCTGGCTGGAACGCCTGCAAATTGCAAAAACGTTAGGCTTCTAATTTGTCGAAATGTCTGTATATGAAAATTACGATAGCCTGTCGAGCCTCAACTGGAGCCGCGACCATAGTCTGGCGCTTGTCaa                                                                                                                      >  2:611029/1‑150 (MQ=255)
                                 ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAAt                                                                                                                     >  2:1628248/3‑150 (MQ=255)
                                 ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGAAGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCGGGTCAAt                                                                                                                     >  2:1270902/3‑150 (MQ=255)
                                   aGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc                                                                                                                   <  2:1494394/150‑1 (MQ=255)
                                   aGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCGTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTTTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc                                                                                                                   <  2:1395684/150‑1 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGa                                                                                                                 >  1:1299726/1‑150 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                >  2:1250512/1‑150 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGAACGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                <  2:6467/150‑1 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACGACTGCAACTGGCAAAAACGTTAGGCATCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGAGCTGGTCAATGCGAt                                                                                                                >  1:1121640/1‑150 (MQ=255)
                                      gcTCCCCGCCGGGGAGTGCTGGCTGGAACCCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                >  2:708955/1‑150 (MQ=255)
                                       cTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAtt                                                                                                               >  1:1265707/1‑150 (MQ=255)
                                         ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgt                                                                                                             >  2:317225/1‑150 (MQ=255)
                                         ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgt                                                                                                             >  1:2049/1‑150 (MQ=255)
                                           ccgccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCATAAACGTTAGGCTGCGATTTTGTCGAATTGTCGGTAGATGATACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgttg                                                                                                           <  1:1247782/150‑1 (MQ=255)
                                            cgccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCATTGCGATTGTTGa                                                                                                          >  1:1287336/1‑150 (MQ=255)
                                                ggggAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAAcc                                                                                                      >  1:1426376/1‑150 (MQ=255)
                                                     gTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGgcgt                                                                                                 <  2:155173/150‑1 (MQ=255)
                                                        ctggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCTATTGTTGAAACCGgcgtgcg                                                                                              >  2:121867/1‑150 (MQ=255)
                                                        ctggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGgcgtgcg                                                                                              >  2:496623/1‑150 (MQ=255)
                                                                     cTGCAACTGGCAAAATCGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGGCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtg                                                                                 <  1:714510/150‑1 (MQ=255)
                                                                     cTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtg                                                                                 <  2:916364/150‑1 (MQ=255)
                                                                       gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                               >  1:529320/1‑150 (MQ=255)
                                                                       gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                               <  2:1012551/150‑1 (MQ=255)
                                                                        cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                              <  2:1379624/150‑1 (MQ=255)
                                                                        cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                              <  1:317225/150‑1 (MQ=255)
                                                                         aaCTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGcc                                                                             <  1:1416014/150‑1 (MQ=255)
                                                                             ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                         >  2:314375/1‑150 (MQ=255)
                                                                             ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                         >  1:1055927/1‑150 (MQ=255)
                                                                               cAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGAGAGCAGCGTCAGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTg                                                                       >  2:203323/1‑150 (MQ=255)
                                                                                aaaaaCGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGc                                                                      >  2:818607/1‑150 (MQ=255)
                                                                                      gTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcg                                                                <  2:619571/150‑1 (MQ=255)
                                                                                        tAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                              >  2:229288/1‑150 (MQ=255)
                                                                                        tAGGCTGCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                              <  1:881141/150‑1 (MQ=255)
                                                                                          ggCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                             <  1:1413578/149‑1 (MQ=255)
                                                                                          ggCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                             >  2:1413578/1‑149 (MQ=255)
                                                                                            cTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGttt                                                          >  1:1044867/1‑150 (MQ=255)
                                                                                              tCGACTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGAGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTcc                                                        <  1:895252/150‑1 (MQ=255)
                                                                                                 aTTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGTTCAATGCGATTGTTGAACCCGGCGTGCGCGAGACGTCCATGTGCCTTTCTGCGAATCGTCGTATccccc                                                     >  2:403618/1‑148 (MQ=255)
                                                                                                 aTTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGc                                                     <  1:227385/150‑1 (MQ=255)
                                                                                                     tgacGAACTGTCGGTCGATGAAACTGACGATCGACTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAACCCGGCGTGCGCGTGCCGTCCATGAGCCTTTCTGCTCATCGTCGTTTCCCGCTggg                                                 <  1:101637/147‑1 (MQ=255)
                                                                                                     tGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTggg                                                 <  2:2049/150‑1 (MQ=255)
                                                                                                      gTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGc                                                >  1:1136465/1‑150 (MQ=255)
                                                                                                       tCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCa                                               >  1:381469/1‑150 (MQ=255)
                                                                                                       tCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCa                                               >  1:549402/1‑150 (MQ=255)
                                                                                                        cGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCCGCGTCTGGCGCGGGTCACTTCGAGTGTTGAAACCGGCGGGAGCGTGCCGGCCATGTGCCGTTCTGCTCCTCGTCGTGTCCCGCTGGGCAg                                              >  2:286148/1‑150 (MQ=255)
                                                                                                         gAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGGTGGGctgg                                             >  1:1160707/1‑147 (MQ=255)
                                                                                                         gAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGt                                             >  2:1270894/1‑150 (MQ=255)
                                                                                                             tGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAg                                         <  2:1549088/150‑1 (MQ=255)
                                                                                                              gTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGGGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGa                                        >  1:1019805/1‑150 (MQ=255)
                                                                                                                cGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATg                                      >  2:1132413/1‑150 (MQ=255)
                                                                                                                cGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATg                                      >  1:1723978/1‑150 (MQ=255)
                                                                                                                 ggTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCTCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGa                                     >  2:1452290/1‑150 (MQ=255)
                                                                                                                  gaaGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGAc                                    <  2:934554/148‑1 (MQ=255)
                                                                                                                        gAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTg                              <  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)
                                                                                                                                                      |                                                                                                                                                    
TTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTA  >  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: