Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,422,992 C→G 5.5% R49R (CGC→CGG)  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,9920CG 5.5% 407.0 / 11.5 128R49R (CGC→CGG) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (70/51);  new base G (5/2);  total (75/53)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 130 total aligned reads displayed.

GATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTA  >  NC_000913/4422843‑4423140
                                                                                                                                                      |                                                                                                                                                    
gATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTACTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCgcgc                                                                                                                                                      >  2:415288/1‑150 (MQ=255)
  tatGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCt                                                                                                                                                    <  2:144835/150‑1 (MQ=255)
  tatGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCt                                                                                                                                                    >  1:1379624/1‑150 (MQ=255)
    tgttttCAAAAAAAATCCGGTTTGGCATCTATGAAAAAGCGCTCCCCGCCGTGGAGTGCTGGCTGGAACGCCTGCAACTGGCAATAACGTTAGGATTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCa                                                                                                                                                  <  1:327550/150‑1 (MQ=255)
     gttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCaa                                                                                                                                                 >  1:1494394/1‑150 (MQ=255)
      ttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAAc                                                                                                                                                >  1:1012551/1‑150 (MQ=255)
       tgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACt                                                                                                                                               <  2:1039207/150‑1 (MQ=255)
        gtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                              >  1:1722768/1‑150 (MQ=255)
        gtCCAAACAAATCCCGCTTGGCATATATGATTAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGGCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTg                                                                                                                                              <  2:695894/150‑1 (MQ=255)
           caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                           >  2:1562563/1‑150 (MQ=255)
           caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAg                                                                                                                                           >  2:227385/1‑150 (MQ=255)
           caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTAGATATTGTCGAAATGTCGGTAGAGG‑AAACTGACGATCGCCTGTCGCGCCTCTACTGGAg                                                                                                                                           >  1:1051841/1‑150 (MQ=255)
           caaacaaaTCCCGCTTGGCATCTATGAAAAAGATCTCCCCGCCACGGAGTCCTGGCAGGAACGCCTCCAACTGGCAAAAACGTAAGGCTTCGAAATTGTCGAAATGTAGGAAGATG‑AAAATGACGATCGCCTGACGCGCCTCAACTGGAg                                                                                                                                           >  1:1225988/1‑150 (MQ=255)
            aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGc                                                                                                                                          >  2:1666671/1‑150 (MQ=255)
             aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCTAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                         >  1:658253/1‑150 (MQ=255)
             aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                         <  1:368264/150‑1 (MQ=255)
             aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                         >  1:43735/1‑150 (MQ=255)
             aacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGcc                                                                                                                                         >  2:534861/1‑150 (MQ=255)
                aaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCcgcg                                                                                                                                      <  1:1501262/150‑1 (MQ=255)
                aaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCcgcg                                                                                                                                      <  2:1058347/150‑1 (MQ=255)
                   tCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGGCTCCAGTTGAGGCGCGaca                                                                                                                                   <  1:444849/150‑3 (MQ=255)
                     cctcgtGGCATCTCTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAGTTTGTCGAAAGGTCGGTAGATC‑ATACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGagcag                                                                                                                                 <  1:1576156/145‑1 (MQ=255)
                         ttGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTc                                                                                                                             <  2:1017003/150‑1 (MQ=255)
                           ggCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTg                                                                                                                           >  2:252559/1‑150 (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)
                                    gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑ACACTGACGATCGCCTGTCGCTCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTc                                                                                                                  >  2:158767/1‑150 (MQ=255)
                                    gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCACCTGGCGCCGGGAGCAGCGGCTGGCGCTGGTc                                                                                                                  >  2:164183/1‑150 (MQ=255)
                                     aaaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCa                                                                                                                 >  2:33782/1‑150 (MQ=255)
                                      aaaaGCGCTCCCCGACGGGGAGTGCTGGCTGGAACGCCTGCAAATTGCAAAAACGTTAGGCTTCTAATTTGTCGAAATGTCTGTATATG‑AAAATTACGATAGCCTGTCGAGCCTCAACTGGAGCCGCGACCATAGTCTGGCGCTTGTCaa                                                                                                                >  2:611029/1‑150 (MQ=255)
                                       ccaGCACTCCCCGGCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTAGGTAGATG‑AAACTGACGATCGCCTGACGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAAt                                                                                                               >  2:1071202/3‑150 (MQ=255)
                                       ccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAAt                                                                                                               >  2:1628248/3‑150 (MQ=255)
                                        caGCACTCCCCGGCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATg                                                                                                              >  2:400093/2‑150 (MQ=255)
                                         aGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCGTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTTTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGc                                                                                                             <  2:1395684/150‑1 (MQ=255)
                                            gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGa                                                                                                           >  1:1299726/1‑150 (MQ=255)
                                            gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                          >  2:1250512/1‑150 (MQ=255)
                                            gcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGAACGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGAt                                                                                                          <  2:6467/150‑1 (MQ=255)
                                            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                                                                                                       >  1:2049/1‑150 (MQ=255)
                                                  cgccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCATTGCGATTGTTGa                                                                                                    >  1:1287336/1‑150 (MQ=255)
                                                      ggggAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAAcc                                                                                                >  1:1426376/1‑150 (MQ=255)
                                                              ctggctggAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGgcgtgcg                                                                                        >  2:496623/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                                                                         <  2:1012551/150‑1 (MQ=255)
                                                                             gCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCAtgtg                                                                         >  1:529320/1‑150 (MQ=255)
                                                                              cAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGc                                                                        <  2:1379624/150‑1 (MQ=255)
                                                                                   ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGCGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCttta                                                                   >  2:1502295/1‑149 (MQ=255)
                                                                                   ggCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTc                                                                   >  1:1055927/1‑150 (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)
                                                                                              tAGGCTGCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtc                                                        <  1:881141/150‑1 (MQ=255)
                                                                                                ggCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                       >  2:1413578/1‑149 (MQ=255)
                                                                                                ggCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCAtcgtcg                                                       <  1:1413578/149‑1 (MQ=255)
                                                                                                    tCGATTTTGTCGATATGTGGGTAGAGG‑AATCTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGTGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCTTCGTCGTTTcc                                                  <  1:901309/150‑1 (MQ=255)
                                                                                                    tCGACTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGAGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTcc                                                  <  1:895252/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)
                                                                                                             tCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCa                                         >  1:549402/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)
                                                                                                               gAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGt                                       >  2:1270894/1‑150 (MQ=255)
                                                                                                                   tGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAg                                   <  2:1549088/150‑1 (MQ=255)
                                                                                                                    gTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGGGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGa                                  >  1:1019805/1‑150 (MQ=255)
                                                                                                                      cGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATg                                >  1:1723978/1‑150 (MQ=255)
                                                                                                                      cGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATg                                >  2:1132413/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)
                                                                                                                                    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)
                                                                                                                                          tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg             <  1:1270902/150‑1 (MQ=255)
                                                                                                                                          tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGATGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg             <  1:979712/150‑1 (MQ=255)
                                                                                                                                          tCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGGCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAggg             <  2:1265707/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)
                                                                                                                                               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)
                                                                                                                                                     gcCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTTGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTCCACCGCGTCGGTTCCCGCTGGGCAGTGAAGATGCCGCGGGGCGGGGGCCGGGGCTGGAGAtga  >  2:870968/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATG‑AAACTGACGATCGCCTGTCGCGCCTCAACTGGAGCCGCGAGCAGCGTCTGGCGCTGGTCAATGCGATTGTTGAAACCGGCGTGCGCGTGCCGTCCATGTGCCTTTCTGCTCATCGTCGTTTCCCGCTGGGCAGTGAAGATGACGCGGTGCGGGCGCAGGGGCTGGAGATTA  >  NC_000913/4422843‑4423140

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