Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,330,396:1 +T 5.3% intergenic (+15/+21) 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,3961.T 5.3% 652.9 / 31.2 152intergenic (+15/+21)sohB/yciNinner membrane protein, S49 peptidase family protein/DUF2498 protein YciN
Reads supporting (aligned to +/- strand):  ref base . (74/70);  new base T (3/5);  total (77/75)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 155 total aligned reads displayed.

AAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCA  >  NC_000913/1330248‑1330544
                                                                                                                                                     |                                                                                                                                                    
aaGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAgg                                                                                                                                                       <  1:399728/149‑1 (MQ=255)
   gCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTa                                                                                                                                                   >  2:1282673/1‑149 (MQ=255)
    ccGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTaa                                                                                                                                                  >  4:227632/1‑149 (MQ=255)
    ccGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTaa                                                                                                                                                  <  2:929426/149‑1 (MQ=255)
    ccGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTaa                                                                                                                                                  >  4:923706/1‑149 (MQ=255)
     cGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAg                                                                                                                                                 >  2:751902/1‑149 (MQ=255)
      gTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGa                                                                                                                                                >  4:1348243/1‑149 (MQ=255)
      gTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGa                                                                                                                                                <  1:736188/149‑1 (MQ=255)
              gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGCAAAACGCGAGGTCTTAGCCTCGCGt                                                                                                                                        >  3:1459253/1‑149 (MQ=255)
              gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGt                                                                                                                                        >  4:685532/1‑149 (MQ=255)
              gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGt                                                                                                                                        >  4:1389578/1‑149 (MQ=255)
                cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGttt                                                                                                                                      <  3:1343046/149‑1 (MQ=255)
                cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGttt                                                                                                                                      <  1:435995/149‑1 (MQ=255)
                cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGttt                                                                                                                                      >  3:1413241/1‑149 (MQ=255)
                     tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGCAAAACGCGAGGTCTTAGCCTCGCGtttgtctt                                                                                                                                 <  3:692275/149‑6 (MQ=38)
                     tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCtt                                                                                                                                 >  4:530018/1‑149 (MQ=255)
                     tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCtt                                                                                                                                 <  4:1224122/149‑1 (MQ=255)
                        gCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGTTTGTCtttta                                                                                                                              <  1:787174/149‑3 (MQ=255)
                        gCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTaa                                                                                                                              >  4:863809/1‑149 (MQ=255)
                           atatGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCa                                                                                                                           >  1:641326/1‑149 (MQ=255)
                              tGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAg                                                                                                                                                 <  4:1192173/124‑1 (MQ=255)
                              tGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAAcc                                                                                                                        >  1:174661/1‑149 (MQ=255)
                               gCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCa                                                                                                                       <  3:923706/149‑1 (MQ=255)
                                 aGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAga                                                                                                                     >  3:200830/1‑149 (MQ=255)
                                     tAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAgatgat                                                                                                                 >  3:1190327/1‑149 (MQ=255)
                                      aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGata                                                                                                                <  2:976705/149‑1 (MQ=255)
                                      aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGata                                                                                                                >  1:598414/1‑149 (MQ=255)
                                       gaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGAGCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                               >  2:1351458/2‑149 (MQ=255)
                                       aaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                               >  2:1509691/1‑149 (MQ=255)
                                           aCTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATAttttt                                                                                                           >  2:186648/1‑149 (MQ=255)
                                             tCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCt                                                                                                         >  3:421456/1‑149 (MQ=255)
                                             tCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCt                                                                                                         <  1:1509691/149‑1 (MQ=255)
                                               aTTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGCCCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGa                                                                                                       >  3:482539/1‑149 (MQ=255)
                                                 tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     >  2:1330389/1‑149 (MQ=255)
                                                 tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     <  4:469211/149‑1 (MQ=255)
                                                 tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     <  4:375804/149‑1 (MQ=255)
                                                  gACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAg                                                                                                    >  2:1172298/1‑149 (MQ=255)
                                                     cGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGcac                                                                                                 >  2:1500823/1‑149 (MQ=255)
                                                     cGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGcac                                                                                                 >  2:650693/1‑149 (MQ=255)
                                                        ttCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATg                                                                                              <  1:941659/149‑1 (MQ=255)
                                                          cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAACGACAAAGGCGAGGCTAAGACCTCGCGTTTTGCTTTAATTAACCAGATGATATTTTTCTGAAAGCATATggg                                                                                            >  4:1128047/1‑149 (MQ=255)
                                                          cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATggg                                                                                            >  2:628357/1‑149 (MQ=255)
                                                           aCCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                           >  4:251077/1‑149 (MQ=255)
                                                           aCCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                           >  2:378125/1‑149 (MQ=255)
                                                             cGGCAGCGCGGCAGATAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCa                                                                                         <  1:10349/149‑1 (MQ=255)
                                                              ggCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                        <  2:598414/149‑1 (MQ=255)
                                                                    gcgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGttt                                                                                  <  3:863809/149‑1 (MQ=255)
                                                                        cAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              <  3:836383/149‑1 (MQ=255)
                                                                        cAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              <  2:318254/149‑1 (MQ=255)
                                                                         agagagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACa                                                                             <  3:1383827/149‑1 (MQ=255)
                                                                           agagCGCCGATCGATTGTCGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTTAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACata                                                                           <  1:785172/149‑1 (MQ=255)
                                                                            gagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACatat                                                                          <  3:980757/149‑1 (MQ=255)
                                                                                gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTaaa                                                                      <  4:680627/149‑1 (MQ=255)
                                                                                gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATAGGCCAGGTGTTTGAACATATTaaa                                                                      >  4:1390858/1‑149 (MQ=255)
                                                                                   gatcgatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATAtt                                                                         >  2:1572663/1‑143 (MQ=255)
                                                                                       gatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATgg                                                                                             >  3:1200342/1‑119 (MQ=255)
                                                                                       gatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATgg                                                                                             <  4:1200342/119‑1 (MQ=255)
                                                                                       gatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCAATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACcgcg                                                               >  2:1238691/1‑149 (MQ=255)
                                                                                           gttgCTACGCTGGTGGCAGCGGGGTCAATAGACATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGttt                                                                                                                                      >  4:963955/1‑74 (MQ=25)
                                                                                           gttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAGCCTCGCGttt                                                                                                                                      <  3:963955/74‑1 (MQ=37)
                                                                                           gttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGc                                                           <  1:1477015/149‑1 (MQ=255)
                                                                                              gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAAt                                                                                                                             >  4:1530565/1‑80 (MQ=255)
                                                                                              gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCttt                                                        <  1:1282673/149‑1 (MQ=255)
                                                                                              gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                       <  4:933115/149‑1 (MQ=255)
                                                                                               cTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                       <  1:1199571/149‑1 (MQ=255)
                                                                                                 aCGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTgg                                                     <  1:1172298/149‑1 (MQ=255)
                                                                                                   gCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              >  2:1461048/1‑122 (MQ=255)
                                                                                                   gCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              <  1:1461048/122‑1 (MQ=255)
                                                                                                     tggtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGt                                                 <  1:757820/149‑1 (MQ=255)
                                                                                                       gtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                               <  4:1413241/149‑1 (MQ=255)
                                                                                                       gtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                               <  2:549429/149‑1 (MQ=255)
                                                                                                        tggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTAGCTGTTgg                                              >  3:19379/1‑149 (MQ=255)
                                                                                                              cGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCATTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATcc                                        <  2:450284/149‑1 (MQ=255)
                                                                                                              cGGGGTCAAAAGACATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATAATATTTATCTGAAAGCACATGGGCCAGGTATTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATcc                                        <  1:378125/149‑1 (MQ=255)
                                                                                                                gggTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCtt                                      <  4:1190327/149‑1 (MQ=255)
                                                                                                                  gTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGt                                    <  3:934556/149‑1 (MQ=255)
                                                                                                                   tCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                   >  2:671994/1‑149 (MQ=255)
                                                                                                                   tCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                   >  1:98875/1‑149 (MQ=255)
                                                                                                                    cAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTc                                  <  3:615912/149‑1 (MQ=255)
                                                                                                                     aaaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCa                                 >  4:1457429/1‑149 (MQ=255)
                                                                                                                      aaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                        <  3:897176/93‑1 (MQ=255)
                                                                                                                      aaaGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                        >  4:897176/1‑93 (MQ=255)
                                                                                                                          ccATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     <  1:1082386/76‑1 (MQ=255)
                                                                                                                            aTTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                          >  4:506931/1‑149 (MQ=255)
                                                                                                                            aTTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                          >  1:1230752/1‑149 (MQ=255)
                                                                                                                             ttGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                       <  2:1112717/119‑1 (MQ=255)
                                                                                                                                aTGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAAATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCTTCTAAAAagta                      <  4:482539/149‑1 (MQ=255)
                                                                                                                                   tAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACTAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTc                   >  1:500164/1‑149 (MQ=255)
                                                                                                                                    aaaaGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                  >  2:872853/1‑149 (MQ=255)
                                                                                                                                     aaaGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGc                 >  4:1280908/1‑149 (MQ=255)
                                                                                                                                       aGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCg               <  2:775132/149‑1 (MQ=255)
                                                                                                                                         aCAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGt             >  3:147523/1‑149 (MQ=255)
                                                                                                                                           aaaCGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTaa           >  4:1568367/1‑149 (MQ=255)
                                                                                                                                           aaaCGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTaa           >  3:1034553/1‑149 (MQ=255)
                                                                                                                                               gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       >  1:1350734/1‑149 (MQ=255)
                                                                                                                                               gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       <  2:500164/149‑1 (MQ=255)
                                                                                                                                               gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       <  4:200830/149‑1 (MQ=255)
                                                                                                                                               gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       <  3:476026/149‑1 (MQ=255)
                                                                                                                                                   ggCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGc   >  3:557451/1‑149 (MQ=255)
                                                                                                                                                    gCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCa  <  1:671994/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCA  >  NC_000913/1330248‑1330544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: