Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,201 C→A 5.9% A130E (GCA→GAA)  yfeD → DUF1323 family putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2010CA 5.9% 1756.5 / 12.6 534A130E (GCA→GAA) yfeDDUF1323 family putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (246/239);  new base A (2/47);  total (248/286)
Fisher's exact test for biased strand distribution p-value = 1.09e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 541 total aligned reads displayed.

ATACGCTCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATC  >  NC_000913/2519052‑2519349
                                                                                                                                                     |                                                                                                                                                    
aTACGCTCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGc                                                                                                                                                      <  1:230330/150‑1 (MQ=255)
    gCTCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCataa                                                                                                                                                  <  2:5479493/150‑1 (MQ=255)
     cTCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCataaa                                                                                                                                                 <  2:5016282/150‑1 (MQ=255)
      tCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCataaat                                                                                                                                                <  2:10913740/150‑1 (MQ=255)
      tCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCataaat                                                                                                                                                <  1:2133557/150‑1 (MQ=255)
        cTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCataaataa                                                                                                                                              <  2:7415177/150‑1 (MQ=255)
         tGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCAtaaataaa                                                                                                                                             <  1:4029597/150‑1 (MQ=255)
           gAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGCCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATaaaaa                                                                                                                                           >  2:58730/1‑150 (MQ=255)
            aGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATaaaaaa                                                                                                                                          <  1:4862164/150‑1 (MQ=255)
             gCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAAc                                                                                                                                         <  1:9766464/150‑1 (MQ=255)
                cgcgCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGc                                                                                                                                      >  1:621611/1‑150 (MQ=255)
                 gcgcCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCa                                                                                                                                     >  1:10205313/1‑150 (MQ=255)
                 gcgcCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCa                                                                                                                                     >  2:9685685/1‑150 (MQ=255)
                 gcgcCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCa                                                                                                                                     >  1:6809287/1‑150 (MQ=255)
                 gcgcCCCGGCGTATCGACAAATCATCAACGCGAGAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTGTGTCACGGGAAGGAATTCGCACTTTCCTCGCCCGTCTCGATCTTGACGAATCAGCATAAATAAAAAACGGCa                                                                                                                                     >  2:1216498/1‑150 (MQ=255)
                   gcCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAgg                                                                                                                                   >  1:2517078/1‑150 (MQ=255)
                   gcCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAgg                                                                                                                                   >  1:6646676/1‑150 (MQ=255)
                             aTCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtctc                                                                                                                         >  2:5641938/1‑150 (MQ=255)
                              tCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTcc                                                                                                                        >  1:7330203/1‑150 (MQ=255)
                              tCGACAAATCAACAACGCGATATATAATATGACAGACATTGAGCAGTAAAAAGTCGCGAAATTTTTGTCACGTGAAGTAATTCGCAATTTCCTCGCCCGTCCCGATATTGCCGAATCAGAATAAATAAAAAACGGCAGGATATTATCTcc                                                                                                                        <  2:2904516/150‑1 (MQ=255)
                               ctaCAAATCATCGACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCAGGTGAAGGATTTCGCAATTTCCTCGCCCGTCTCGTTATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCt                                                                                                                       <  2:1888704/148‑1 (MQ=255)
                                gACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                      >  2:8900450/1‑150 (MQ=255)
                                 aCAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGc                                                                                                                     <  2:5973143/150‑1 (MQ=255)
                                   aaaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                                                   >  2:6220217/1‑150 (MQ=255)
                                            aCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttt                                                                                                          >  1:7087750/1‑150 (MQ=255)
                                            aCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATAttt                                                                                                          >  1:1884456/1‑150 (MQ=255)
                                              gcgATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                                                        >  1:10148836/1‑150 (MQ=255)
                                              gcgATACATAATCTGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTAACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                                                        <  2:2438880/150‑1 (MQ=255)
                                                 atagataATATGACAGACATTGAGCAGGAGAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                                                     >  2:9261273/1‑150 (MQ=255)
                                                    gataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAATTAAAAAACGGCAGGATTTTTTCTCCTGGGGTTTATCttttttgacgg                                                                                                  >  2:4330512/1‑144 (MQ=255)
                                                    gataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                                                  >  2:3458723/1‑150 (MQ=255)
                                                       aataTGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTTTGTTTTTAAACGCTaa                                                                                               >  1:3945289/1‑150 (MQ=255)
                                                          atGACAGACATTGCGCAGTAAAATGTCGCGAATTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCTATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                                            <  1:784228/150‑1 (MQ=255)
                                                          atGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                                            >  2:5722863/1‑150 (MQ=255)
                                                           tGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTGTATCTTTTTACACGCTAATTAc                                                                                           >  1:1061891/1‑150 (MQ=255)
                                                            tacagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGATGTATTTCGCAATTTCCTCGCCCTTCTCTATTTTGACTAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                                          <  2:16323/149‑1 (MQ=255)
                                                            gacagacaTTGAGCAGGAAAAAGTCGCGCATTTTTTGGCTGGTGATGGAATTCGCTATTTCCTCGCCCGTCTCGTTTTTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                                          <  2:5586743/150‑1 (MQ=255)
                                                              cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                        <  2:1770780/150‑1 (MQ=255)
                                                              cagacaGGGAGCAGCAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                        <  2:7012854/150‑1 (MQ=255)
                                                              cagaaaTTGAGCAGTAAAAAGTCGCGCAATTTTTGTCACGTGAAGTAATTCGAAATTTCCTCGCCCGTCTCTATTTTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                        <  2:6520754/150‑1 (MQ=255)
                                                              cacacaTTGAGCAGGAAAAAGCCGCGCATTTTTGGTCACGTGACGTAATCCTCATTTTCCTCGCCCTTCTCGATTTTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                        <  1:164406/147‑1 (MQ=255)
                                                               agacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTCAGTActgct                                                                                       >  2:6689248/1‑150 (MQ=255)
                                                                   aTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTGTATCTTTTTACACGCTAATTAATGCTGAtt                                                                                   >  2:4728471/1‑150 (MQ=255)
                                                                    tttagaaGGAAAAATTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCTATATTGACTAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                                                  <  2:3901131/144‑1 (MQ=255)
                                                                      gAGCAGGAAAAAGGCGCGCAATTTTGGTCACGTGAAGGAATTCGCAAGTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                                                <  2:7804742/150‑1 (MQ=255)
                                                                       aGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                                               <  2:3334467/150‑1 (MQ=255)
                                                                       aGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTGTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                                               >  1:9099713/1‑150 (MQ=255)
                                                                       aGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                                               <  2:2148673/150‑1 (MQ=255)
                                                                        gccggAAAATGTCGCGCAATTTGTGTCACGTGACGGAAGTCGCAATTGCCTCGCCCGTCTCGATATTGACGAATCCGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                                              <  1:5797437/147‑1 (MQ=255)
                                                                         cAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAATCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTAAACGCTAATTACTGCTGATTTTcgcg                                                                             >  1:154831/1‑150 (MQ=255)
                                                                          aGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                                            >  1:2373986/1‑150 (MQ=255)
                                                                           ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                                           >  1:8092867/1‑150 (MQ=255)
                                                                             aaaaaGTCGCGAAATTTTTGTCACTTGACGGAATTCGAAATTTCCTCGCCCGTCTCGATATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                                         <  1:9572446/150‑1 (MQ=255)
                                                                                aaGTCGCCCAATTTTTGCCACGTGAAGGAACTCGCAATTTCCCCGCCCGGCTCTATATTGACGAATCAGCATAAAAAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                                                      <  1:3396924/150‑1 (MQ=255)
                                                                                 aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAATCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                                                     >  2:6358324/1‑150 (MQ=255)
                                                                                 aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                                                     >  2:3562061/1‑150 (MQ=255)
                                                                                     gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCTTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                                                 <  1:1246293/150‑1 (MQ=255)
                                                                                        cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                                              >  1:6983903/1‑150 (MQ=255)
                                                                                           tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                                           <  2:9009396/150‑1 (MQ=255)
                                                                                           tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                                           >  2:4598561/1‑150 (MQ=255)
                                                                                            ttttGTCACGTGATGGAATCCGCATTGTCCTCGCCCGTCTCTATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                                          <  2:9467812/150‑1 (MQ=255)
                                                                                             tttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                                                         >  1:4398656/1‑150 (MQ=255)
                                                                                               tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                                                       <  2:10490882/150‑1 (MQ=255)
                                                                                                 tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                                                     >  1:9893459/1‑150 (MQ=255)
                                                                                                  cACGTTACGGATTTCGCATTTTCCGCGCCCGCCTCGATATTGCCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                                                    <  1:3758783/150‑1 (MQ=255)
                                                                                                  cACGTGAAGGAATTCGCAATTTCCTCGCCCTTCTCGATATTGACTAATAAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                                                    <  1:853468/150‑1 (MQ=255)
                                                                                                    cGTGAATTAATTCGAAATTTCCTCGCCCTAATCTATATTGACTAATCACAAAAAAAAAAAAACGGCAGGATATTATCTCCTGCCGTTTAACTTATTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                  <  2:4832485/150‑1 (MQ=255)
                                                                                                    cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                                                  >  2:6080105/1‑150 (MQ=255)
                                                                                                      tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTGCGGCAATAAAATCCAGCGCGTTGTTGa                                                >  2:102884/1‑150 (MQ=255)
                                                                                                       gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATACAGCGCTTTGTTGAt                                               >  1:8744772/1‑150 (MQ=255)
                                                                                                        aaGTAATCCGCATTTTCCTCGCCCTTCTCTATTTTGACTAACAACAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                                              <  1:9317330/150‑1 (MQ=255)
                                                                                                           gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTGTATCTTTTTACACGCGAATTACTGCTGATTTTCGCGTTCCGCAATAAAATCCAGCGCGTTGTTGATACGc                                           >  1:2705714/1‑150 (MQ=255)
                                                                                                               tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                                       >  1:10368542/1‑150 (MQ=255)
                                                                                                               tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                                       <  1:5257418/150‑1 (MQ=255)
                                                                                                               tCGCAATTTCCTCGCCCGCCTCGATATTGACGAATCAGCAGAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                                       <  2:6853922/150‑1 (MQ=255)
                                                                                                                    aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTg                                  <  2:86941/150‑1 (MQ=255)
                                                                                                                    aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTg                                  <  2:9799934/150‑1 (MQ=255)
                                                                                                                    aTTGCCTCGCCCGTCTCGATATTGACGAAAAAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTg                                  <  1:3128681/150‑1 (MQ=255)
                                                                                                                     tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTGTATATTTTTAAACGCTAATTACTGCTGATTTTCGAGTTCAGCAAAAAAATCCAGCGCATTGTTGATAAGCTCGATACTGc                                 >  1:7622746/1‑150 (MQ=255)
                                                                                                                      ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg                                <  2:2578555/150‑1 (MQ=255)
                                                                                                                           cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGTCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt                           <  1:3591652/150‑1 (MQ=255)
                                                                                                                           cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACAGCGGGTCt                           >  2:20264/1‑150 (MQ=255)
                                                                                                                           cGCCCGTCTCGATATTGACGAACCAGCATAAATAAAAAACGGCACGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt                           <  2:8541019/150‑1 (MQ=255)
                                                                                                                            gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATGTTTTTAAAAGCTAATTACTTCTGATTTTAGCGTTGAGCAATAAAATACAGCGCTTTTTTGTTACTCTCGATAATGCGGGTCtt                          >  1:2965865/1‑150 (MQ=255)
                                                                                                                            gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCtt                          <  1:10823640/150‑1 (MQ=255)
                                                                                                                             cccGTCTCTATATTGACTAATCACCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa                         <  1:8881009/150‑1 (MQ=255)
                                                                                                                             cccGTCTCGATATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa                         <  2:3053803/150‑1 (MQ=255)
                                                                                                                               cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                       <  2:3658650/150‑1 (MQ=255)
                                                                                                                                  ctcGATTTTGACGAATCAGCATAAATAAAAAACGGCTGGATTTTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAAt                    <  2:8175574/150‑1 (MQ=255)
                                                                                                                                   tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAAtt                   <  1:5562782/150‑1 (MQ=255)
                                                                                                                                     gATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGc                 <  2:10821004/150‑1 (MQ=255)
                                                                                                                                       tatTGACTAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGt               <  2:2249980/150‑1 (MQ=255)
                                                                                                                                         ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGa             <  1:6349898/150‑1 (MQ=255)
                                                                                                                                          tGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGaa            <  1:6141597/150‑1 (MQ=255)
                                                                                                                                          tGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTAAACGCTAAATAATTCTGATTTTGGCGTTCAGCAATAAAAAACAGCGCGTTTTTGTTACGCGCGATAATGCGGGTCTTAACAATTGCGTGaa            >  1:2806969/1‑150 (MQ=255)
                                                                                                                                           gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAc           >  2:7983662/1‑150 (MQ=255)
                                                                                                                                              gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGt        <  2:1411679/150‑1 (MQ=255)
                                                                                                                                                    gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATc  <  2:7191260/150‑1 (MQ=255)
                                                                                                                                                    gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATc  >  2:5321326/1‑150 (MQ=255)
                                                                                                                                                    gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATATTTTTACACGCTAATTACTGGTGATTTTCGCGTTCTGCAATAAATTCCAGCGCTTTGTTGTTTCGGTCGATTCTTCGGGTGTTAGGAATTGCGTGTACGGTAACATc  >  2:10153892/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATACGCTCCTGGAGCACGCGCCCCGGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATC  >  NC_000913/2519052‑2519349

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: