Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,330,399 A→T 7.4% intergenic (+18/+18) sohB → / ← yciN inner membrane protein, S49 peptidase family protein/DUF2498 protein YciN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,330,3990AT 7.4% 409.1 / 14.8 142intergenic (+18/+18)sohB/yciNinner membrane protein, S49 peptidase family protein/DUF2498 protein YciN
Reads supporting (aligned to +/- strand):  ref base A (61/70);  new base T (5/6);  total (66/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

GGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCC  >  NC_000913/1330250‑1330548
                                                                                                                                                      |                                                                                                                                                     
ggCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCTGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTa                                                                                                                                                       <  1:1335921/150‑1 (MQ=255)
 gCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATAGACCGATTCACCGGCAGCGCGGAAGAGAGCGCGGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTaa                                                                                                                                                      >  1:1466755/1‑150 (MQ=255)
       aaGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTc                                                                                                                                                <  2:773044/150‑1 (MQ=255)
        aGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTcg                                                                                                                                               >  2:187986/1‑150 (MQ=255)
         gTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTcgc                                                                                                                                              <  2:1018731/150‑1 (MQ=255)
         gTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTcgc                                                                                                                                              >  1:34828/1‑150 (MQ=255)
           ggTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGt                                                                                                                                            >  2:1911480/1‑150 (MQ=255)
           ggTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGGCAAAACCAATGGATTAAAAATACAAACCCCAGC‑CTAAGAcaacccgt                                                                                                                                            >  1:531858/1‑143 (MQ=255)
              cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGtttt                                                                                                                                         <  2:1056298/150‑1 (MQ=255)
              cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGtttt                                                                                                                                         >  2:1484480/1‑150 (MQ=255)
               aaTGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAG‑CCTCGCGtttgt                                                                                                                                        <  2:141804/150‑3 (MQ=255)
                 tGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCt                                                                                                                                      <  2:368068/150‑1 (MQ=255)
                  gTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCAGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCtt                                                                                                                                     <  2:1380458/150‑1 (MQ=255)
                   tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCttt                                                                                                                                    <  1:661469/150‑1 (MQ=255)
                   tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCttt                                                                                                                                    <  2:125255/150‑1 (MQ=255)
                     cGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCCACAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGCCCTCGCGTTTTGCTTTaa                                                                                                                                  >  2:306518/1‑150 (MQ=255)
                       cTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGCGCGCCGATCGCTGGTTGCTACGCTGGTGGCAGCGGGGTCAAAATCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTGTAATc                                                                                                                                >  2:684621/1‑150 (MQ=255)
                       cTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAACGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATc                                                                                                                                >  1:186330/1‑150 (MQ=255)
                       cTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATc                                                                                                                                <  1:445844/150‑1 (MQ=255)
                                gCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAgatg                                                                                                                       >  1:1382931/1‑150 (MQ=255)
                                  gTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CAAAGACCTCGCGTTTTGCTTTAATAAACCAgatgat                                                                                                                     >  1:765051/1‑150 (MQ=255)
                                  gTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGGGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTGTAATCAACCAgatgat                                                                                                                     >  2:1165280/1‑150 (MQ=255)
                                    aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                                   >  1:1556122/1‑150 (MQ=255)
                                      aCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATAttt                                                                                                                 <  2:859612/150‑1 (MQ=255)
                                       cGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATAtttt                                                                                                                >  1:1087447/1‑150 (MQ=255)
                                         aCTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTc                                                                                                              <  2:1576054/150‑1 (MQ=255)
                                         aCTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTc                                                                                                              >  1:551297/1‑150 (MQ=255)
                                           tCATTGACCGATTCACCGGCAGCGCGGCATAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTg                                                                                                            <  1:1079186/150‑1 (MQ=255)
                                             aTTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaa                                                                                                          <  1:897869/150‑1 (MQ=255)
                                             aTTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGAAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaa                                                                                                          <  2:1259508/150‑1 (MQ=255)
                                             aTTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTCCTACGCTCGTGGCAGCGGGCTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGATTTGCTTTAATCAACCAGATCATATTTTTCTGaa                                                                                                          >  1:634624/1‑150 (MQ=255)
                                                 aCCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGca                                                                                                      <  1:262286/150‑1 (MQ=255)
                                                 aCCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGca                                                                                                      >  1:1743211/1‑150 (MQ=255)
                                                  ccGATTCAACGGCAGCGCGGCAGAGATCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGTGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAAACAACCAGATGATAATTTTCTGAAAGcac                                                                                                     >  1:1178468/1‑150 (MQ=255)
                                                       tCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATggg                                                                                                <  2:1641159/150‑1 (MQ=255)
                                                       tCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATggg                                                                                                <  1:1277687/150‑1 (MQ=255)
                                                       tCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATggg                                                                                                <  1:452541/150‑1 (MQ=255)
                                                        cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                               >  2:883487/1‑150 (MQ=255)
                                                        cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                               >  1:1219650/1‑150 (MQ=255)
                                                         aCCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGcc                                                                                              <  1:1984244/150‑1 (MQ=255)
                                                          ccGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCa                                                                                             >  2:132076/1‑150 (MQ=255)
                                                          ccGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCa                                                                                             <  1:785272/150‑1 (MQ=255)
                                                                  gcgGCAGAGAGCGCCGATCGATTGTTTCTACGCTGGTGGCATCGGTGTCAAAAGCCATTGATGTAAAATACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTg                                                                                     <  2:2000347/150‑1 (MQ=255)
                                                                  gcgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTg                                                                                     >  2:640140/1‑150 (MQ=255)
                                                                  gcgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTg                                                                                     <  2:1992577/150‑1 (MQ=255)
                                                                  gcgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTg                                                                                     <  2:835945/150‑1 (MQ=255)
                                                                    ggCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGaa                                                                                   <  1:706108/150‑1 (MQ=255)
                                                                     gCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                                  >  2:612722/1‑150 (MQ=255)
                                                                      cAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACa                                                                                 <  2:1411424/150‑1 (MQ=255)
                                                                         agagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACatat                                                                              <  1:800751/150‑1 (MQ=255)
                                                                            gcgcCGATCGATTGTGGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTaa                                                                           <  1:1429082/150‑1 (MQ=255)
                                                                            gcgcCCATCGATTTTTGCTACGCTGGTGGCAGCGGTGTCAAAAGCCATTGATGTAAAAGACACACGCGATG‑CTTAGACCTCGCGTTTTGCTTTAATCAACCCGATGTTATTTTTCTGACAGCACATGGGCCAGGTGTTTGAACATATTaa                                                                           <  1:83494/150‑1 (MQ=255)
                                                                             cgcCGATCGATTGTTGCTACGCTGGTGTCAGCGGGGTCAAAAGACATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTaaa                                                                          >  2:517072/1‑150 (MQ=255)
                                                                              gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAac                                                                         <  1:870381/150‑1 (MQ=255)
                                                                              gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAac                                                                         >  1:1184574/1‑150 (MQ=255)
                                                                               ccGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAaca                                                                        >  1:1175198/1‑150 (MQ=255)
                                                                               ccGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATAAACCAGATGATATTTTTCTGAAAGCACATGGGACAGGTGTTTGAACATATTAAaca                                                                        >  2:237849/1‑150 (MQ=255)
                                                                                       ttgttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTg                                                                <  1:1527148/150‑1 (MQ=255)
                                                                                            gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGTGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                           <  1:1622493/150‑1 (MQ=255)
                                                                                             cTACGCTGGTGGCAGCGGTGTCAAAAGCCATTTATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTg                                                          <  2:1169503/150‑1 (MQ=255)
                                                                                              tACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTgg                                                         >  1:964736/1‑150 (MQ=255)
                                                                                                cGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCt                                                       >  1:1055633/1‑150 (MQ=255)
                                                                                                   tggtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGtt                                                    <  2:39918/150‑1 (MQ=255)
                                                                                                    ggtgtCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                                   <  2:524253/150‑1 (MQ=255)
                                                                                                     gtggCAGCGGGGTCAACAGAGATTGATGTAAAATACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTTAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGGGCTTTTGGCTGTTgg                                                  <  2:1265133/150‑1 (MQ=255)
                                                                                                      tggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGc                                                 <  2:967590/150‑1 (MQ=255)
                                                                                                         cAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAAt                                              >  1:839035/1‑150 (MQ=255)
                                                                                                            cGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCt                                           <  1:872046/150‑1 (MQ=255)
                                                                                                            cGGGGGCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCt                                           >  1:416221/1‑150 (MQ=255)
                                                                                                               ggTCAACAGCCATTGATTTAAAAGACAAAATCGAGG‑ATAAGACATCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCATTTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGt                                        <  1:142026/150‑1 (MQ=255)
                                                                                                               ggTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCgtgt                                        >  2:1043558/1‑150 (MQ=255)
                                                                                                                gTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                       <  2:1668598/150‑1 (MQ=255)
                                                                                                                gTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                       <  1:1621544/150‑1 (MQ=255)
                                                                                                                 tCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTc                                      <  2:1708675/150‑1 (MQ=255)
                                                                                                                  cAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCATCCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGGTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCa                                     >  1:25543/1‑150 (MQ=255)
                                                                                                                    aaaGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGTCCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATc                                   <  2:666632/150‑1 (MQ=255)
                                                                                                                      aGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                 <  2:1080358/150‑1 (MQ=255)
                                                                                                                      aGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                 <  2:1087297/150‑1 (MQ=255)
                                                                                                                      aGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTa                                 >  1:1861513/1‑150 (MQ=255)
                                                                                                                         cATTGATGTAAACGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                              <  2:80169/150‑1 (MQ=255)
                                                                                                                         cATTGATGTAAAAGAGAAACTCGAGG‑CTAAGACCTAGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATTGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                              <  2:1870967/150‑1 (MQ=255)
                                                                                                                         cATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                              <  1:38752/150‑1 (MQ=255)
                                                                                                                         cATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                              >  2:886381/1‑150 (MQ=255)
                                                                                                                         cATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATAATTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                              <  1:642033/150‑1 (MQ=255)
                                                                                                                           ttGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAag                            <  1:1412858/150‑1 (MQ=255)
                                                                                                                              aTGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagtag                         <  2:34828/150‑1 (MQ=255)
                                                                                                                                gTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTc                       >  2:942736/1‑150 (MQ=255)
                                                                                                                                 tAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAAAACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                      >  1:1443831/1‑150 (MQ=255)
                                                                                                                                   aaagcaaaaCGCGAGGTCTTAG‑CCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGcc                    >  1:564011/7‑150 (MQ=255)
                                                                                                                                    aaGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCg                   <  1:1206781/150‑1 (MQ=255)
                                                                                                                                      gACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTTAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGt                 >  1:143891/1‑150 (MQ=255)
                                                                                                                                           aCGCGAGG‑CTACGACCTGGCGTTTTGCTTTAATCAACAAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGGGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATa            <  1:1949220/150‑1 (MQ=255)
                                                                                                                                           aCGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTTGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATa            >  1:1451329/1‑150 (MQ=255)
                                                                                                                                           aCGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATa            >  2:619603/1‑150 (MQ=255)
                                                                                                                                              cgAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCa         >  1:310818/1‑150 (MQ=255)
                                                                                                                                               gAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAg        <  1:1804143/150‑1 (MQ=255)
                                                                                                                                                 gg‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGTAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCa      >  1:602790/1‑150 (MQ=255)
                                                                                                                                                     tAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGc   >  2:1713297/1‑150 (MQ=255)
                                                                                                                                                      aaGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTACACACCGCGGTGCTTTTGGCTGTGGGCAATCCTTGTACATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGcc  >  2:239413/1‑150 (MQ=255)
                                                                                                                                                      aaGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAAACCTTGTTCATCAAAAAAGTAGTCGCCGGTAAATACCAGCACGcc  >  1:1864713/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCC  >  NC_000913/1330250‑1330548

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑