Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,423,004 C→G 5.5% S53R (AGC→AGG)  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,0040CG 5.5% 384.3 / 8.5 123S53R (AGC→AGG) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (62/54);  new base G (5/2);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 128 total aligned reads displayed.

AAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTAAAGC  >  NC_000913/4422855‑4423150
                                                                                                                                                      |                                                                                                                                                  
aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGc                                                                                                                                                    >  2:1666671/1‑150 (MQ=255)
 aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCTAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                                   >  1:658253/1‑150 (MQ=255)
 aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                                   >  2:534861/1‑150 (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)
         cctcgtGGCATCTCTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAAGGTCGGTAGATC‑ATACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGagcag                                                                                                                                           <  1:1576156/145‑1 (MQ=255)
             ttGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                                                                                                                                       <  2:1017003/150‑1 (MQ=255)
                gCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTgg                                                                                                                                    >  2:1481295/1‑150 (MQ=255)
                gCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCGGGCTGGAACGCCGGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTgg                                                                                                                                    <  1:1063054/150‑1 (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)
                     tATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGGCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTg                                                                                                                               <  2:1001479/150‑1 (MQ=255)
                      aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCggg                                                                                                                              >  1:1512597/1‑150 (MQ=255)
                      aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTgg                                                                                                                              >  2:1672748/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)
                         aaaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCa                                                                                                                           >  2:33782/1‑150 (MQ=255)
                          aaaaGCGCTCCCCGACGGGGAGTGCTGGCTGGAACGCCTGCAAATTGCAAAAACGTTAGGCTTCTAATTTGTCGAAATGTCTGTATATG‑AAAATTACGATAGCCTGTCGAGCCTCAACTGGAGCCGCGACCATAGTCTGGCGCTTGTCaa                                                                                                                          >  2:611029/1‑150 (MQ=255)
                           ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGAAG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCGGGTCAAt                                                                                                                         >  2:1270902/3‑150 (MQ=255)
                            caGCACTCCCCGGCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATg                                                                                                                        >  2:400093/2‑150 (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)
                                  tCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AACCTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttg                                                                                                                  >  1:1439460/1‑150 (MQ=255)
                                   ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAttgt                                                                                                                 >  2:317225/1‑150 (MQ=255)
                                                  ctggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCTATTGTTGAAACCGgcgtgcg                                                                                                  >  2:121867/1‑150 (MQ=255)
                                                    ggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATTGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGgcgtgcgcg                                                                                                <  1:681335/150‑1 (MQ=255)
                                                             gCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCa                                                                                       <  1:1024617/150‑1 (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)
                                                                  cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                                  <  1:317225/150‑1 (MQ=255)
                                                                  cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                                  <  2:1379624/150‑1 (MQ=255)
                                                                   aaCTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGcc                                                                                 <  1:1416014/150‑1 (MQ=255)
                                                                      tGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACCGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCttt                                                                              >  2:1596154/1‑150 (MQ=255)
                                                                       ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCttta                                                                             >  2:1502295/1‑149 (MQ=255)
                                                                       ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                             >  2:314375/1‑150 (MQ=255)
                                                                         cAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTg                                                                           >  1:484221/1‑150 (MQ=255)
                                                                         cAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGAGAGCAGCGTCAGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTg                                                                           >  2:203323/1‑150 (MQ=255)
                                                                          aaaaaCGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGc                                                                          >  2:818607/1‑150 (MQ=255)
                                                                                  tAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                                  >  2:229288/1‑150 (MQ=255)
                                                                                  tAGGCTGCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                                  <  1:881141/150‑1 (MQ=255)
                                                                                   aGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAACCCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                                 >  2:475811/1‑150 (MQ=255)
                                                                                    ggCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                                 >  2:1413578/1‑149 (MQ=255)
                                                                                    ggCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                                 <  1:1413578/149‑1 (MQ=255)
                                                                                      cTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGttt                                                              >  1:1044867/1‑150 (MQ=255)
                                                                                        tCGATTTTGTCGATATGTGGGTAGAGG‑AATCTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGTGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCTTCGTCGTTTcc                                                            <  1:901309/150‑1 (MQ=255)
                                                                                        tCGACTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGAGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTcc                                                            <  1:895252/150‑1 (MQ=255)
                                                                                          gATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCg                                                          <  1:187565/150‑1 (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)
                                                                                                  cGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCCGCGTCTGGCGCGGGTCACTTCGAGTGTTGAAACCGGCGGGAGCGTGCCGGCCATGTGCCGTTCTGCTCCTCGTCGTGTCCCGCTGGGCAg                                                  >  2:286148/1‑150 (MQ=255)
                                                                                                   gAAATGTCGGTAGATGAAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGCCGTTTCCCGCTGGGCAg                                                  <  2:1299726/150‑1 (MQ=255)
                                                                                                   gAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGGTGGGctgg                                                 >  1:1160707/1‑147 (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:1389453/1‑150 (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)
                                                                                                                              tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg                       <  1:1270902/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:43735/150‑1 (MQ=255)
                                                                                                                                   tGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTgg                  <  2:1055927/150‑1 (MQ=255)
                                                                                                                                   gctCGCGGCTCCAGTTGAGGCGCGATCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGTCGTGCGCGTGCAGTCCATGTGCCTTTCTGCTCATAGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCTGGCGCAGGGGCTgg                  >  2:241555/3‑150 (MQ=255)
                                                                                                                                    gTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCAGCTCATCGTCGTTTCCCGCTGGGCAGTGAACATGACCCGGTGCGGGCGCAGGGGCTggt                 >  1:793647/1‑149 (MQ=255)
                                                                                                                                    ctCGCGGCTCGAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGga                 >  1:610870/2‑150 (MQ=255)
                                                                                                                                    ctCGCGGCTCCAGTTGAGGCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCGTCGTCGTTTCCCGCTGGGCAGTGAAGATGACCCGGTGCGGGCGCAGGGGCgggg                 >  1:170657/2‑149 (MQ=255)
                                                                                                                                     tCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGgag                <  1:502928/150‑1 (MQ=255)
                                                                                                                                           cTCAACTGGAGCCGCGAGCAGCGTCTGGCGGTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTGATCGTCGTTTCCCGCAGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCAGGAGATTATg          >  1:778245/1‑150 (MQ=255)
                                                                                                                                            gcAACTGGAGCCGCGAGCAGCGGCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGc         <  1:1562563/149‑1 (MQ=255)
                                                                                                                                             cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGTCGTGCGCGTGCCGTCCATGTGCCTTTATGCTCATCGTCGTTTACCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCg        >  2:398134/1‑150 (MQ=255)
                                                                                                                                             cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCg        <  1:314375/150‑1 (MQ=255)
                                                                                                                                             cAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATAGTCGTTTACCGCTGGGCAGTGAAGATGAAGAGGTGCGTGCGCAGGGGCTGGAGATTATGCg        >  2:391938/1‑150 (MQ=255)
                                                                                                                                              aaCTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGt       <  1:1250512/150‑1 (MQ=255)
                                                                                                                                               aCTGGAGCCGCGCGCAGCGTCGGGCGCTGGTCAATGCGATTGTTGAACCCTGCGTGCGCGTGCCCTCCATGTGCCTTTCTGCTCAGCGTAGTTTCCCGCTGGGCAGTGAAGAGGACGCGGTGCGGGCGCTGGGGCTGGAGATGATGCgga      >  2:671679/1‑148 (MQ=255)
                                                                                                                                               aCTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTa      <  2:484221/150‑1 (MQ=255)
                                                                                                                                                   gAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTAAAGc  <  2:946064/150‑1 (MQ=255)
                                                                                                                                                   gAGCCGCGAGCAGCGGCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTAAAGc  <  1:1542890/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTATGCGTAAAGC  >  NC_000913/4422855‑4423150

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: ‑