Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,423,000 G→T 5.4% W52L (TGG→TTG)  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,423,0000GT 5.4% 364.2 / 9.5 131W52L (TGG→TTG) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (71/53);  new base T (5/2);  total (76/55)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 133 total aligned reads displayed.

GTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTA  >  NC_000913/4422851‑4423146
                                                                                                                                                      |                                                                                                                                                  
gtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                                    >  1:1722768/1‑150 (MQ=255)
gtCCAAACAAATCCCGCTTGGCATATATGATTAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGGCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                                    <  2:695894/150‑1 (MQ=255)
   caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                                 >  2:1562563/1‑150 (MQ=255)
   caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                                 >  2:227385/1‑150 (MQ=255)
   caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                                 >  1:1372656/1‑150 (MQ=255)
   caaacaaaTCCCGCTTGGCATCTATGAAAAAGATCTCCCCGCCACGGAGTCCTGGCAGGAACGCCTCCAACTGGCAAAAACGTAAGGCTTCGAAATTGTCGAAATGTAGGAAGATGAAAATGACGATCGCCTGACGCGCCTCAACTGGAg                                                                                                                                                 >  1:1225988/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)
     aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                               >  2:534861/1‑150 (MQ=255)
        aaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCcgcg                                                                                                                                            <  2:1058347/150‑1 (MQ=255)
        aaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCcgcg                                                                                                                                            <  1:1501262/150‑1 (MQ=255)
          aTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGGCTCCAGTTGAGGCGCGac                                                                                                                                          <  1:905777/150‑2 (MQ=255)
                 ttGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                                                                                                                                   <  2:1017003/150‑1 (MQ=255)
                   ggCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTg                                                                                                                                 >  2:252559/1‑150 (MQ=255)
                    gCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCGGGCTGGAACGCCGGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTgg                                                                                                                                <  1:1063054/150‑1 (MQ=255)
                     cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgc                                                                                                                               <  1:993562/150‑1 (MQ=255)
                     cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGgc                                                                                                                               >  2:1542890/1‑150 (MQ=255)
                          aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCggg                                                                                                                          >  1:1512597/1‑150 (MQ=255)
                          aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTgg                                                                                                                          >  2:1672748/1‑150 (MQ=255)
                            gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGACACTGACGATCGCCTGTCGCTCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTc                                                                                                                        >  2:158767/1‑150 (MQ=255)
                            gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTc                                                                                                                        >  1:619571/1‑150 (MQ=255)
                            gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATGTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTc                                                                                                                        >  2:100943/1‑150 (MQ=255)
                             aaaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCa                                                                                                                       >  2:33782/1‑150 (MQ=255)
                              aaaaGCGCTCCCCGACGGGGAGTGCTGGCTGGAACGCCTGCAAATTGCAAAAACGTTAGGCTTCTAATTTGTCGAAATGTCTGTATATGAAAATTACGATAGCCTGTCGAGCCTCAACTGGAGCCGCGACCATAGTCTGGCGCTTGTCaa                                                                                                                      >  2:611029/1‑150 (MQ=255)
                               ccaGCACTCCCCGGCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTAGGTAGATGAAACTGACGATCGCCTGACGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAAt                                                                                                                     >  2:1071202/3‑150 (MQ=255)
                               ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAAt                                                                                                                     >  2:1628248/3‑150 (MQ=255)
                               ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGAAGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCGGGTCAAt                                                                                                                     >  2:1270902/3‑150 (MQ=255)
                                caGCACTCCCCGGCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATg                                                                                                                    >  2:400093/2‑150 (MQ=255)
                                 aGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc                                                                                                                   <  2:1494394/150‑1 (MQ=255)
                                    gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                >  2:1250512/1‑150 (MQ=255)
                                    gcTCCCCGCCGGGGAGTGCTGGCTGGAACGACTGCAACTGGCAAAAACGTTAGGCATCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGAGCTGGTCAATGCGAt                                                                                                                >  1:1121640/1‑150 (MQ=255)
                                    gcTCCCCGCCGGGGAGTGCTGGCTGGAACCCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                                >  2:708955/1‑150 (MQ=255)
                                     cTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAtt                                                                                                               >  1:1265707/1‑150 (MQ=255)
                                      tCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAACCTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttg                                                                                                              >  1:1439460/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)
                                                                 gCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCa                                                                                   <  1:1024617/150‑1 (MQ=255)
                                                                   cTGCAACTGGCAAAATCGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGGCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtg                                                                                 <  1:714510/150‑1 (MQ=255)
                                                                   cTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtg                                                                                 <  2:916364/150‑1 (MQ=255)
                                                                     gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                               <  2:1012551/150‑1 (MQ=255)
                                                                     gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                               >  1:529320/1‑150 (MQ=255)
                                                                      cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                              <  2:1379624/150‑1 (MQ=255)
                                                                      cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                              <  1:317225/150‑1 (MQ=255)
                                                                          tGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACCGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCttt                                                                          >  2:1596154/1‑150 (MQ=255)
                                                                           ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCttta                                                                         >  2:1502295/1‑149 (MQ=255)
                                                                           ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                         >  2:314375/1‑150 (MQ=255)
                                                                           ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                         >  1:1055927/1‑150 (MQ=255)
                                                                             cAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTg                                                                       >  1:484221/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)
                                                                                       aGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAACCCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                             >  2:475811/1‑150 (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)
                                                                                              gATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCg                                                      <  1:187565/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)
                                                                                                      cGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCCGCGTCTGGCGCGGGTCACTTCGAGTGTTGAAACCGGCGGGAGCGTGCCGGCCATGTGCCGTTCTGCTCCTCGTCGTGTCCCGCTGGGCAg                                              >  2:286148/1‑150 (MQ=255)
                                                                                                       gAAATGTCGGTAGATGAAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGCCGTTTCCCGCTGGGCAg                                              <  2:1299726/150‑1 (MQ=255)
                                                                                                       gAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGGTGGGctgg                                             >  1:1160707/1‑147 (MQ=255)
                                                                                                            gTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGGGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGa                                        >  1:1019805/1‑150 (MQ=255)
                                                                                                               ggTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCTCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGa                                     >  2:1452290/1‑150 (MQ=255)
                                                                                                                gaaGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGAc                                    <  2:934554/148‑1 (MQ=255)
                                                                                                                           cTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCggg                          >  1:1389453/1‑150 (MQ=255)
                                                                                                                           cTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCggg                          >  1:423881/1‑150 (MQ=255)
                                                                                                                             gACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGgcg                        >  2:406147/1‑150 (MQ=255)
                                                                                                                               cGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCa                      >  2:1596069/1‑150 (MQ=255)
                                                                                                                                 aTCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAgg                    <  2:381469/150‑1 (MQ=255)
                                                                                                                                  tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg                   <  1:1270902/150‑1 (MQ=255)
                                                                                                                                  tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGATGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg                   <  1:979712/150‑1 (MQ=255)
                                                                                                                                  tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGGCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg                   <  2:1265707/150‑1 (MQ=255)
                                                                                                                                   cGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAgggg                  <  2:1439460/150‑1 (MQ=255)
                                                                                                                                    gCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGc                 <  1:415288/150‑1 (MQ=255)
                                                                                                                                    gCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCGGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGc                 >  1:1452610/1‑150 (MQ=255)
                                                                                                                                       tGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTgg              <  2:1055927/150‑1 (MQ=255)
                                                                                                                                       gctCGCGGCTCCAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCGATGTGCCTGTCTGCTCATCGTCGTTTCCCGCTGGGGAGTGAAGATGCCGCGGTGCGGGCGCAGGGGCTgg              >  1:1602066/3‑150 (MQ=255)
                                                                                                                                        ctCGCGGCTCGAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGga             >  1:610870/2‑150 (MQ=255)
                                                                                                                                        ctCGCGGCTCCAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCGTCGTCGTTTCCCGCTGGGCAGTGAAGATGACCCGGTGCGGGCGCAGGGGCgggg             >  1:170657/2‑149 (MQ=255)
                                                                                                                                        ctCGCGGCTCCAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACAGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTGCCGCTGGGCAGGGAAGATGACGAGGTGCGGGCGCAGGGGCTggc             >  1:875919/2‑149 (MQ=255)
                                                                                                                                         tCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGgag            <  1:502928/150‑1 (MQ=255)
                                                                                                                                             gcCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTTGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTCCACCGCGTCGGTTCCCGCTGGGCAGTGAAGATGCCGCGGGGCGGGGGCCGGGGCTGGAGAtga        >  2:870968/1‑148 (MQ=255)
                                                                                                                                               cTCAACTGGAGCCGCGAGCAGCGTCTGGCGGTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTGATCGTCGTTTCCCGCAGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCAGGAGATTATg      >  1:778245/1‑150 (MQ=255)
                                                                                                                                               cTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATg      <  2:1722768/150‑1 (MQ=255)
                                                                                                                                                gcAACTGGAGCCGCGAGCAGCGGCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGc     <  1:1562563/149‑1 (MQ=255)
                                                                                                                                                 cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCg    <  1:314375/150‑1 (MQ=255)
                                                                                                                                                 cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATAGTCGTTTACCGCTGGGCAGTGAAGATGAAGAGGTGCGTGCGCAGGGGCTGGAGATTATGCg    >  2:391938/1‑150 (MQ=255)
                                                                                                                                                  aaCGGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGGGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGGCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGt   <  1:708955/150‑1 (MQ=255)
                                                                                                                                                   aCTGGAGCCGCGCGCAGCGTCGGGCGCTGGTCAATGCGATTGTTGAACCCTGCGTGCGCGTGCCCTCCATGTGCCTTTCTGCTCAGCGTAGTTTCCCGCTGGGCAGTGAAGAGGACGCGGTGCGGGCGCTGGGGCTGGAGATGATGCgga  >  2:671679/1‑148 (MQ=255)
                                                                                                                                                   aCTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTa  <  2:484221/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTA  >  NC_000913/4422851‑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: