Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 518,334 T→G 5.8% intergenic (‑55/+6) ybbN ← / ← ybbO DnaK co‑chaperone, thioredoxin‑like protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,3340TG 5.8% 1662.8 / 7.9 568intergenic (‑55/+6)ybbN/ybbODnaK co‑chaperone, thioredoxin‑like protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (226/283);  new base G (47/2);  total (283/285)
Fisher's exact test for biased strand distribution p-value = 1.37e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 589 total aligned reads displayed.

ACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTC  >  NC_000913/518186‑518483
                                                                                                                                                    |                                                                                                                                                      
  cAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCa                                                                                                                                                     <  1:248773/150‑1 (MQ=255)
  cAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGGCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCGTCa                                                                                                                                                     >  2:7173622/1‑150 (MQ=255)
  cAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGGTGTCGATATTTTCTTTTACATGGGGGCTTAAGCGCGCGCGTcc                                                                                                                                                     >  1:3687369/1‑149 (MQ=255)
   aaaaaTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCaa                                                                                                                                                    >  2:4936940/1‑150 (MQ=255)
    aaaaTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAAc                                                                                                                                                   <  2:6474459/150‑1 (MQ=255)
      aaTAGAACAGCACCGGAGTGGTCATCGACAGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTc                                                                                                                                                 >  1:7513599/1‑150 (MQ=255)
       aTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTGGGTATTTTGTTTTACATGGGGGCTTAAGCGCGCGCGTCAACTcc                                                                                                                                                >  1:6346881/1‑149 (MQ=255)
          gAACAGCACCGGAGTGGTCATCGACTGTTCTAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑Accc                                                                                                                                            <  2:2673176/150‑1 (MQ=255)
            aCAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTGGGTATTTTGTTTGTCATGGGGGCTTAAGCGCGCGCTTCAACTC‑Accccg                                                                                                                                          >  1:1172706/1‑149 (MQ=255)
             cAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTATTTTACATGGGGGCTTAAGCGCGCGCTTAAACAC‑ACCCCTg                                                                                                                                         >  2:1844400/1‑150 (MQ=255)
                cACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGGTGTCGATATTTTCTTTGACATTGGGGCGTAAGCGCGCGCTTCAACAC‑ACCCCGGCaa                                                                                                                                      >  1:9204536/1‑150 (MQ=255)
                cACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCGGTTGTCGATATTTTCTTTTAAATGGGGGCGTAAGCGCGCGCTTCAAAGC‑ACCCCAGCac                                                                                                                                      >  1:2464444/1‑149 (MQ=255)
                        tGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGt                                                                                                                              <  2:6972028/150‑1 (MQ=255)
                        tGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGt                                                                                                                              <  1:4264075/150‑1 (MQ=255)
                        tGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGt                                                                                                                              <  1:8776927/150‑1 (MQ=255)
                          gTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTGTATGTTTGCGATATTTTAATTGTCAATGGGGCTTAAGCGCGCCCATAAACTC‑ACCCCGGCAAAAatttttac                                                                                                                            >  1:1745164/1‑146 (MQ=255)
                          gTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTcc                                                                                                                            <  2:2031542/150‑1 (MQ=255)
                             aTCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACAAGGGGGCTTAAGCGCGCGCGTCAACAC‑ACCCCAGCAAAATTTTGTCCaac                                                                                                                         >  1:337279/1‑148 (MQ=255)
                              tCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTTACATGGGGGCTTAAGCGCGCGCTTAAACAC‑ACCCCTGCAAAATTTTGTCCAtaa                                                                                                                        >  1:4999343/1‑148 (MQ=255)
                              tCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCa                                                                                                                        <  2:9283580/150‑1 (MQ=255)
                                gACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCAcg                                                                                                                      >  2:9086423/1‑150 (MQ=255)
                                gACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTAGATATTTTCTTTGAAATGGGGGCTTAATCGCGCGCTTCAAAAA‑ACCCCAGCAAAATTTTTTACATCaag                                                                                                                      >  2:3495505/1‑148 (MQ=255)
                                 aCTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCAcgc                                                                                                                     <  1:26291/150‑1 (MQ=255)
                                       cAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCg                                                                                                               <  2:5356223/150‑1 (MQ=255)
                                        aaGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTGTCTGTTGTCGATATTTTATTTTTCATTGGGGCTTAATCGCCCGCGTCAACAC‑ACCCCAGCCAAATTTTTTACAACACGcccccccg                                                                                                              >  1:2059311/1‑148 (MQ=255)
                                        aaGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGGTGTGGATATTTTGTTTTTCATGGGGGCGTAAGCGCGCCCGTCAACTC‑ACCCCGGCAATATTTTGTCCATCACGCGcccccg                                                                                                              >  1:2295642/1‑148 (MQ=255)
                                         aGAAACTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCTATATTTTGTCCATCACGCGCCCAGgc                                                                                                             <  2:2936445/150‑1 (MQ=255)
                                          gAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑AACCCTGCAATATATTGTACATCACGCGCCCCCgca                                                                                                            >  2:7378995/1‑150 (MQ=255)
                                                 ttGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGgcg                                                                                                     <  2:5181060/150‑1 (MQ=255)
                                                  tGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTGTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGgcgc                                                                                                    <  1:7044552/150‑1 (MQ=255)
                                                  tGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGgcgc                                                                                                    <  2:5925687/150‑1 (MQ=255)
                                                   gCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCt                                                                                                   <  2:387932/150‑1 (MQ=255)
                                                   gCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCt                                                                                                   <  1:4082861/150‑1 (MQ=255)
                                                    cAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCtt                                                                                                  >  2:4535335/1‑150 (MQ=255)
                                                     aGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCATTTAAATTGGGGCTTAAGCGCGCGCATCAAATC‑AACCCTGCAAAATTTTTTACATCACGCGCCCCCGCAGCAGGCGCTTa                                                                                                 >  2:1772893/1‑150 (MQ=255)
                                                      ggTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTaa                                                                                                <  1:6200337/150‑1 (MQ=255)
                                                       gTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGAAATTTTATTTTTCAATGGGGGTTAATCGCGCGCGTCAAAAA‑ACCCCAGCAAAAATTTTTACAACAACCGCCCCCCCACCAGGCGCATAAg                                                                                               >  1:5114086/1‑150 (MQ=255)
                                                        ttAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACAC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGc                                                                                              >  2:5765852/1‑150 (MQ=255)
                                                            aTTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCa                                                                                          <  1:4890428/150‑1 (MQ=255)
                                                                gTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTGTGTTGTCGATATTTTATTTTACATGGGGGCGTTATAGCCCGCATCAACACAACCCCGCCAAAATTTTTACCA‑CAAGCGCCCCCGCCCCACGCGCGTAAGCAcccaaac                                                                                      >  1:5664061/1‑145 (MQ=255)
                                                                gTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTATCGGTTGTCGATATTTTCTTTGACAAGGGGGCTTAAGCGCCCGCGTCAACAC‑AACCCTGCAATATTTTGTCCATAACGCGCCCCGGCACCAGGCGCGTAATCAcccatac                                                                                      >  2:7895420/1‑145 (MQ=255)
                                                                 ttAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCGGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTAcc                                                                                     >  2:2112848/1‑150 (MQ=255)
                                                                    aTGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGcc                                                                                  <  2:3726939/150‑1 (MQ=255)
                                                                      gTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCa                                                                                <  2:4648181/150‑1 (MQ=255)
                                                                      gTTGACAATATTTTCTACGGACATGGAGTCGCTCTATTTTTTCGATATTTTATTTTAAAATGGGGCTTAATCGCGCGCGTAAAAAC‑AACCCAGCAAAATATTTTGCAAAAACCCCCCCCGCCGCAGGCGCGTAAGCACCAATAACGcccc                                                                                >  1:6694838/1‑149 (MQ=255)
                                                                        tGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAgg                                                                              <  1:8955371/150‑1 (MQ=255)
                                                                           cAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCa                                                                           >  2:7711671/1‑150 (MQ=255)
                                                                            aaTATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCAc                                                                          <  2:8492817/150‑1 (MQ=255)
                                                                              tatTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCa                                                                        <  2:861356/150‑1 (MQ=255)
                                                                              tatTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCa                                                                        <  1:1284685/150‑1 (MQ=255)
                                                                                ttttCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGc                                                                      <  2:7786763/150‑1 (MQ=255)
                                                                                     tACGGACATGGAGTCGCTCGCTGTTGTCGATATTTTCGTTGACATGGGGGCTTAAGCGCGCCCTTCCACGC‑ACCCCGGCAATAATTTGTCCATCACCCCCCCCGGCAGCAGGCGCTTAAGCACCAATACCCCCCACGGCACCcgcggccc                                                                 >  2:8828917/1‑146 (MQ=255)
                                                                                           cATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATa                                                           <  2:8001590/150‑1 (MQ=255)
                                                                                              ggAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAgcg                                                        <  1:9043160/150‑1 (MQ=255)
                                                                                              ggAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAgcg                                                        >  2:2329169/1‑150 (MQ=255)
                                                                                                  tCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCCTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATc                                                    >  2:5664644/1‑150 (MQ=255)
                                                                                                     ctctctGTTGTCGATATTTTCTTTGACATGGGGGATTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTc                                                 <  1:3901658/150‑1 (MQ=255)
                                                                                                      tctctGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCg                                                <  2:4081614/150‑1 (MQ=255)
                                                                                                       ctctGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCgg                                               >  2:6591940/1‑150 (MQ=255)
                                                                                                          tgttgtCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCtt                                            <  2:8094867/150‑1 (MQ=255)
                                                                                                            ttgtCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTtct                                          >  2:4249027/1‑150 (MQ=255)
                                                                                                              gtCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCg                                        <  1:1018979/150‑1 (MQ=255)
                                                                                                               tCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGc                                       >  2:1366069/1‑150 (MQ=255)
                                                                                                               tCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGc                                       <  2:3687369/150‑1 (MQ=255)
                                                                                                               tCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCAGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCAATAACGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGc                                       >  1:7909553/1‑150 (MQ=255)
                                                                                                                  atatTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaa                                    >  2:8302778/1‑150 (MQ=255)
                                                                                                                   tatTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaat                                   >  1:64529/1‑150 (MQ=255)
                                                                                                                     ttttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaataa                                 >  1:3221691/1‑150 (MQ=255)
                                                                                                                      tttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATaaa                                >  1:8652009/1‑150 (MQ=255)
                                                                                                                         cTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGc                             <  1:9266622/150‑1 (MQ=255)
                                                                                                                         cTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGc                             <  1:8296699/150‑1 (MQ=255)
                                                                                                                          tttGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCa                            >  1:5764619/1‑150 (MQ=255)
                                                                                                                          tttGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCa                            >  1:1899674/1‑150 (MQ=255)
                                                                                                                          tttGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCa                            >  2:6207516/1‑150 (MQ=255)
                                                                                                                           ttGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCAt                           <  1:5222282/150‑1 (MQ=255)
                                                                                                                           ttGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCAt                           >  1:8520615/1‑150 (MQ=255)
                                                                                                                             gACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGc                                                                      >  1:2152551/1‑105 (MQ=255)
                                                                                                                              aCATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGc                        >  2:334304/1‑150 (MQ=255)
                                                                                                                                 tGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTa                     <  2:1138322/150‑1 (MQ=255)
                                                                                                                                  gggggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTAc                    <  2:3120280/150‑1 (MQ=255)
                                                                                                                                   ggggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACt                   >  1:284325/1‑150 (MQ=255)
                                                                                                                                   ggggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACt                   >  1:5364189/1‑150 (MQ=255)
                                                                                                                                   ggggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACt                   >  1:4498730/1‑150 (MQ=255)
                                                                                                                                   ggggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACt                   >  2:215663/1‑150 (MQ=255)
                                                                                                                                    gggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACtt                  <  1:1040038/150‑1 (MQ=255)
                                                                                                                                    gggCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACtt                  >  1:6877122/1‑150 (MQ=255)
                                                                                                                                       cTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGt               >  2:3181702/1‑150 (MQ=255)
                                                                                                                                        ttAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAgcgc                                                       <  2:681122/109‑1 (MQ=255)
                                                                                                                                        ttAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTc              <  1:7781041/150‑1 (MQ=255)
                                                                                                                                         tAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGTCAGCAGGCTCTTAAGCACCATTACCGCCCAGGTCAACATCGTCACCGGATAGCGCATCTTGGGCTTCTCGCTAATAAAAGCTTGGCGTACTTTGTcc             <  2:2260476/150‑1 (MQ=255)
                                                                                                                                            gcgcgcgcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTccacc          >  2:7063387/1‑150 (MQ=255)
                                                                                                                                             cgcgcgcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTccacca         >  1:5961183/1‑150 (MQ=255)
                                                                                                                                             cgcgcgcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTccacca         <  1:4301645/150‑1 (MQ=255)
                                                                                                                                              gcgcgcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCcaccac        >  1:3008069/1‑150 (MQ=255)
                                                                                                                                              gcgcgcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCcaccac        >  2:1584323/1‑150 (MQ=255)
                                                                                                                                               cgcgcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCaccacc       <  2:6698281/150‑1 (MQ=255)
                                                                                                                                                  gcTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCt    >  2:4433284/1‑150 (MQ=255)
                                                                                                                                                    ttCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTc  >  2:4825611/1‑150 (MQ=255)
                                                                                                                                                    ttCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTc  >  2:2500036/1‑150 (MQ=255)
                                                                                                                                                    ttCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTc  >  2:6328283/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTC‑ACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTC  >  NC_000913/518186‑518483

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