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.

AAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCA  >  NC_000913/1330248‑1330544
                                                                                                                                                     |                                                                                                                                                    
aaGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAgg                                                                                                                                                       <  1:399728/149‑1 (MQ=255)
   gCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTa                                                                                                                                                   >  2:1282673/1‑149 (MQ=255)
    ccGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTaa                                                                                                                                                  >  4:227632/1‑149 (MQ=255)
    ccGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTaa                                                                                                                                                  <  2:929426/149‑1 (MQ=255)
    ccGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTaa                                                                                                                                                  >  4:923706/1‑149 (MQ=255)
     cGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAg                                                                                                                                                 >  2:751902/1‑149 (MQ=255)
      gTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGa                                                                                                                                                >  4:1348243/1‑149 (MQ=255)
      gTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGa                                                                                                                                                <  1:736188/149‑1 (MQ=255)
              gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGCAAAACGCGAGGTCTTAG‑CCTCGCGt                                                                                                                                        >  3:1459253/1‑149 (MQ=255)
              gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGt                                                                                                                                        >  4:685532/1‑149 (MQ=255)
              gTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGt                                                                                                                                        >  4:1389578/1‑149 (MQ=255)
                cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAG‑CCTCGCGttt                                                                                                                                      <  3:1343046/149‑1 (MQ=255)
                cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAG‑CCTCGCGttt                                                                                                                                      <  1:435995/149‑1 (MQ=255)
                cAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGttt                                                                                                                                      >  3:1413241/1‑149 (MQ=255)
                     tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGCAAAACGCGAGGTCTTAG‑CCTCGCGtttgtctt                                                                                                                                 <  3:692275/149‑6 (MQ=38)
                     tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCtt                                                                                                                                 >  4:530018/1‑149 (MQ=255)
                     tACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCtt                                                                                                                                 <  4:1224122/149‑1 (MQ=255)
                        gCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAG‑CCTCGCGTTTGTCtttta                                                                                                                              <  1:787174/149‑3 (MQ=255)
                        gCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTaa                                                                                                                              >  4:863809/1‑149 (MQ=255)
                           atatGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCa                                                                                                                           >  1:641326/1‑149 (MQ=255)
                              tGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAg                                                                                                                                                 <  4:1192173/124‑1 (MQ=255)
                              tGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAAcc                                                                                                                        >  1:174661/1‑149 (MQ=255)
                               gCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCa                                                                                                                       <  3:923706/149‑1 (MQ=255)
                                 aGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAga                                                                                                                     >  3:200830/1‑149 (MQ=255)
                                     tAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAgatgat                                                                                                                 >  3:1190327/1‑149 (MQ=255)
                                      aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGata                                                                                                                <  2:976705/149‑1 (MQ=255)
                                      aaaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGata                                                                                                                >  1:598414/1‑149 (MQ=255)
                                       gaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGAGCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                               >  2:1351458/2‑149 (MQ=255)
                                       aaCGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGatat                                                                                                               >  2:1509691/1‑149 (MQ=255)
                                           aCTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATAttttt                                                                                                           >  2:186648/1‑149 (MQ=255)
                                             tCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCt                                                                                                         >  3:421456/1‑149 (MQ=255)
                                             tCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCt                                                                                                         <  1:1509691/149‑1 (MQ=255)
                                               aTTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGCCCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGa                                                                                                       >  3:482539/1‑149 (MQ=255)
                                                 tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     >  2:1330389/1‑149 (MQ=255)
                                                 tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     <  4:469211/149‑1 (MQ=255)
                                                 tGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     <  4:375804/149‑1 (MQ=255)
                                                  gACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAg                                                                                                    >  2:1172298/1‑149 (MQ=255)
                                                     cGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGcac                                                                                                 >  2:1500823/1‑149 (MQ=255)
                                                     cGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGcac                                                                                                 >  2:650693/1‑149 (MQ=255)
                                                        ttCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATg                                                                                              <  1:941659/149‑1 (MQ=255)
                                                          cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAACGACAAAGGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATTAACCAGATGATATTTTTCTGAAAGCATATggg                                                                                            >  4:1128047/1‑149 (MQ=255)
                                                          cACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATggg                                                                                            >  2:628357/1‑149 (MQ=255)
                                                           aCCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                           >  4:251077/1‑149 (MQ=255)
                                                           aCCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGc                                                                                           >  2:378125/1‑149 (MQ=255)
                                                             cGGCAGCGCGGCAGATAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCa                                                                                         <  1:10349/149‑1 (MQ=255)
                                                              ggCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                        <  2:598414/149‑1 (MQ=255)
                                                                    gcgGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGttt                                                                                  <  3:863809/149‑1 (MQ=255)
                                                                        cAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              <  3:836383/149‑1 (MQ=255)
                                                                        cAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              <  2:318254/149‑1 (MQ=255)
                                                                         agagagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACa                                                                             <  3:1383827/149‑1 (MQ=255)
                                                                           agagCGCCGATCGATTGTCGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTTAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACata                                                                           <  1:785172/149‑1 (MQ=255)
                                                                            gagCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACatat                                                                          <  3:980757/149‑1 (MQ=255)
                                                                                gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTaaa                                                                      <  4:680627/149‑1 (MQ=255)
                                                                                gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATAGGCCAGGTGTTTGAACATATTaaa                                                                      >  4:1390858/1‑149 (MQ=255)
                                                                                   gatcgatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATAtt                                                                         >  2:1572663/1‑143 (MQ=255)
                                                                                       gatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATgg                                                                                             >  3:1200342/1‑119 (MQ=255)
                                                                                       gatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATgg                                                                                             <  4:1200342/119‑1 (MQ=255)
                                                                                       gatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCAATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACcgcg                                                               >  2:1238691/1‑149 (MQ=255)
                                                                                           gttgCTACGCTGGTGGCAGCGGGGTCAATAGACATTGATGTAAAAGACAAACGCGAGGTCTTAG‑CCTCGCGttt                                                                                                                                      >  4:963955/1‑74 (MQ=25)
                                                                                           gttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGTCTTAG‑CCTCGCGttt                                                                                                                                      <  3:963955/74‑1 (MQ=37)
                                                                                           gttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGc                                                           <  1:1477015/149‑1 (MQ=255)
                                                                                              gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAAt                                                                                                                             >  4:1530565/1‑80 (MQ=255)
                                                                                              gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCttt                                                        <  1:1282673/149‑1 (MQ=255)
                                                                                              gCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGT‑AAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                       <  4:933115/149‑1 (MQ=255)
                                                                                               cTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                       <  1:1199571/149‑1 (MQ=255)
                                                                                                 aCGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTgg                                                     <  1:1172298/149‑1 (MQ=255)
                                                                                                   gCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              >  2:1461048/1‑122 (MQ=255)
                                                                                                   gCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAAc                                                                              <  1:1461048/122‑1 (MQ=255)
                                                                                                     tggtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGt                                                 <  1:757820/149‑1 (MQ=255)
                                                                                                       gtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                               <  4:1413241/149‑1 (MQ=255)
                                                                                                       gtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                               <  2:549429/149‑1 (MQ=255)
                                                                                                        tggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTAGCTGTTgg                                              >  3:19379/1‑149 (MQ=255)
                                                                                                              cGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCATTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATcc                                        <  2:450284/149‑1 (MQ=255)
                                                                                                              cGGGGTCAAAAGACATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATAATATTTATCTGAAAGCACATGGGCCAGGTATTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATcc                                        <  1:378125/149‑1 (MQ=255)
                                                                                                                gggTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCtt                                      <  4:1190327/149‑1 (MQ=255)
                                                                                                                  gTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGt                                    <  3:934556/149‑1 (MQ=255)
                                                                                                                   tCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                   >  2:671994/1‑149 (MQ=255)
                                                                                                                   tCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                   >  1:98875/1‑149 (MQ=255)
                                                                                                                    cAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTc                                  <  3:615912/149‑1 (MQ=255)
                                                                                                                     aaaaGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCa                                 >  4:1457429/1‑149 (MQ=255)
                                                                                                                      aaaGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                        <  3:897176/93‑1 (MQ=255)
                                                                                                                      aaaGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAg                                                                                        >  4:897176/1‑93 (MQ=255)
                                                                                                                          ccATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGaaa                                                                                                     <  1:1082386/76‑1 (MQ=255)
                                                                                                                            aTTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                          >  4:506931/1‑149 (MQ=255)
                                                                                                                            aTTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTaaaa                          >  1:1230752/1‑149 (MQ=255)
                                                                                                                             ttGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCtttt                                                       <  2:1112717/119‑1 (MQ=255)
                                                                                                                                aTGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAAATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCTTCTAAAAagta                      <  4:482539/149‑1 (MQ=255)
                                                                                                                                   tAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACTAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTc                   >  1:500164/1‑149 (MQ=255)
                                                                                                                                    aaaaGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                  >  2:872853/1‑149 (MQ=255)
                                                                                                                                     aaaGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGc                 >  4:1280908/1‑149 (MQ=255)
                                                                                                                                       aGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCg               <  2:775132/149‑1 (MQ=255)
                                                                                                                                         aCAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGt             >  3:147523/1‑149 (MQ=255)
                                                                                                                                           aaaCGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTaa           >  4:1568367/1‑149 (MQ=255)
                                                                                                                                           aaaCGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTaa           >  3:1034553/1‑149 (MQ=255)
                                                                                                                                               gcgAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       >  1:1350734/1‑149 (MQ=255)
                                                                                                                                               gcgAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       <  2:500164/149‑1 (MQ=255)
                                                                                                                                               gcgAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       <  4:200830/149‑1 (MQ=255)
                                                                                                                                               gcgAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATAc       <  3:476026/149‑1 (MQ=255)
                                                                                                                                                   gg‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGc   >  3:557451/1‑149 (MQ=255)
                                                                                                                                                    g‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCa  <  1:671994/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGCCGTGAAGTGGTCAATGTACGCTATATGCAGCGTAAACGACTCATTGACCGATTCACCGGCAGCGCGGCAGAGAGCGCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGG‑CTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCA  >  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: ‑