Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,021 T→C 8.5% intergenic (+25/+11) cysS → / ← ybcI cysteine‑‑tRNA ligase/conserved inner membrane protein YbcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,0210TC 8.5% 439.1 / 23.7 152intergenic (+25/+11)cysS/ybcIcysteine‑‑tRNA ligase/conserved inner membrane protein YbcI
Reads supporting (aligned to +/- strand):  ref base T (61/78);  new base C (5/8);  total (66/86)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 155 total aligned reads displayed.

AACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCACCTG  >  NC_000913/555872‑556169
                                                                                                                                                     |                                                                                                                                                    
aaCAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGttt                                                                                                                                                      >  2:1601860/1‑150 (MQ=255)
   aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCg                                                                                                                                                   >  1:1283263/1‑150 (MQ=255)
   aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCg                                                                                                                                                   <  2:881372/150‑1 (MQ=255)
   aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCg                                                                                                                                                   >  1:410008/1‑150 (MQ=255)
        cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATcc                                                                                                                                              <  1:2352385/150‑1 (MQ=255)
        cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATcc                                                                                                                                              >  2:1166196/1‑150 (MQ=255)
        cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATcc                                                                                                                                              <  2:910744/150‑1 (MQ=255)
         tGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCg                                                                                                                                             <  2:577494/150‑1 (MQ=255)
          ggATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCgg                                                                                                                                            >  2:1377982/1‑150 (MQ=255)
           gATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGt                                                                                                                                           <  1:2275660/150‑1 (MQ=255)
           gATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATCCggc                                                                                                                                           <  1:154773/150‑2 (MQ=255)
           gATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATCCggc                                                                                                                                           <  1:2045277/150‑2 (MQ=255)
           cacgCCCGTAACGCGACAGACTGGGGTGCGGCGGATGCGGCGCGTGAGCGTCTTAACGAGACGGGGATCGTGGTGGAAGCTGGCCCGCAAGGGACCACCCGGCGCCGTAAGGAATTGCGCTATTGCCGGAGGCGAAAACTCGCATCCggc                                                                                                                                           <  2:446372/147‑2 (MQ=255)
            ttGCCCGTAAAGCGAAAGACTTGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACGTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATCCggca                                                                                                                                          <  2:722585/149‑3 (MQ=255)
              gCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTAt                                                                                                                                        >  2:1812139/1‑150 (MQ=255)
               cccGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATc                                                                                                                                       >  1:710747/1‑150 (MQ=255)
                  gTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTc                                                                                                                                    <  2:2420209/150‑1 (MQ=255)
                   tAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCt                                                                                                                                   >  1:711621/1‑150 (MQ=255)
                      aGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTgcg                                                                                                                                <  1:1900746/150‑1 (MQ=255)
                       gCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTgcgc                                                                                                                               >  2:1823845/1‑150 (MQ=255)
                        cGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGcc                                                                                                                              <  1:880779/150‑1 (MQ=255)
                        cGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGcc                                                                                                                              <  1:2242460/150‑1 (MQ=255)
                        cGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTAATGCCGGATGCGAGTTTTAGCATCCGGTTATCGTCTGCGcc                                                                                                                              >  2:1947691/1‑150 (MQ=255)
                         gAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGccc                                                                                                                             >  2:1747047/1‑149 (MQ=255)
                          aaaGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGccac                                                                                                                            >  1:535863/1‑150 (MQ=255)
                          aaaGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGccac                                                                                                                            >  1:2314013/1‑150 (MQ=255)
                          aaaGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGccac                                                                                                                            >  1:1020022/1‑150 (MQ=255)
                             gACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCcaccac                                                                                                                         >  2:660780/1‑150 (MQ=255)
                                 gggCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCacaaca                                                                                                                     >  1:1925809/1‑150 (MQ=255)
                                  ggcggcggcGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACAt                                                                                                                    <  2:535863/150‑1 (MQ=255)
                                    cggcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGATCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTc                                                                                                                  <  2:1239665/150‑1 (MQ=255)
                                      gcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTccc                                                                                                                <  2:1510073/150‑1 (MQ=255)
                                      gcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTccc                                                                                                                <  2:2223394/150‑1 (MQ=255)
                                       gggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCa                                                                                                               <  1:272919/149‑1 (MQ=255)
                                       cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCa                                                                                                               <  1:1479054/150‑1 (MQ=255)
                                       cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCa                                                                                                               <  2:2093672/150‑1 (MQ=255)
                                        ggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCAt                                                                                                              <  1:2342751/150‑1 (MQ=255)
                                         gcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCGGGCGTCGTAAGTACTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCCTc                                                                                                             <  2:1954597/150‑1 (MQ=255)
                                          cggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCa                                                                                                            <  2:323936/150‑1 (MQ=255)
                                           ggATGCGGGGGGTGATGGTCTTAACTAGATGGGGATCGTGCTGGAAGTTGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTTTTGCCGTATGCGAGTTTTCACATCCGGTTATCGTCTGCGCCACCACAACATTCGCATCAg                                                                                                           <  2:2444210/150‑1 (MQ=255)
                                           ggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAg                                                                                                           <  1:386365/150‑1 (MQ=255)
                                            gATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGt                                                                                                          <  2:429517/150‑1 (MQ=255)
                                             aTGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATCGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTa                                                                                                         <  2:2023170/150‑1 (MQ=255)
                                              tGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAg                                                                                                        <  1:487494/150‑1 (MQ=255)
                                               gcggcgCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGc                                                                                                       >  2:898473/1‑150 (MQ=255)
                                                 ggcgCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCAt                                                                                                     <  1:1979381/150‑1 (MQ=255)
                                                      gTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCg                                                                                                >  2:196968/1‑150 (MQ=255)
                                                        gATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCTTTATTTTATTTTTAGGATGAGGTTATAGTTTTAGACACTAAAACATTCTCATCATTATCATccctgta                                                                                              <  2:2384640/150‑5 (MQ=255)
                                                         aTCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCa                                                                                             <  2:1105404/150‑1 (MQ=255)
                                                           cGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAc                                                                                           <  1:371746/150‑1 (MQ=255)
                                                            gTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAcc                                                                                          <  1:597602/150‑1 (MQ=255)
                                                              cTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACcac                                                                                        <  1:225105/150‑1 (MQ=255)
                                                                 aaCGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACAcc                                                                                     >  1:1408769/1‑150 (MQ=255)
                                                                 aaCGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACAcc                                                                                     >  1:2201292/1‑150 (MQ=255)
                                                                       aTGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATc                                                                               <  2:1687382/150‑1 (MQ=255)
                                                                       aTGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATc                                                                               <  1:1823845/150‑1 (MQ=255)
                                                                          gggATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAAt                                                                            <  2:281127/150‑1 (MQ=255)
                                                                              tCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAg                                                                        <  1:1602072/150‑1 (MQ=255)
                                                                              tCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAg                                                                        <  2:1786328/150‑1 (MQ=255)
                                                                               cGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGa                                                                       <  1:1328275/150‑1 (MQ=255)
                                                                                gTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGaa                                                                      >  2:228809/1‑150 (MQ=255)
                                                                                 tGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGaaa                                                                     <  2:410008/150‑1 (MQ=255)
                                                                                     ggAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAcc                                                                                          <  1:1051852/125‑1 (MQ=255)
                                                                                     ggAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAcc                                                                                          >  2:1051852/1‑125 (MQ=255)
                                                                                      gAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATc                                                                >  1:317572/1‑150 (MQ=255)
                                                                                        aGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCAc                                                              <  1:1166927/150‑1 (MQ=255)
                                                                                           tGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTg                                                           >  2:2253458/1‑150 (MQ=255)
                                                                                           tGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGGTTTCGCTTCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTg                                                           <  2:1319317/150‑1 (MQ=255)
                                                                                             gCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGAt                                                         <  1:170082/150‑1 (MQ=255)
                                                                                                 gCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGccc                                                     <  2:916558/150‑1 (MQ=255)
                                                                                                   aaGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGt                                                   <  1:1166196/150‑1 (MQ=255)
                                                                                                    aGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTa                                                  >  2:463174/1‑150 (MQ=255)
                                                                                                     gggACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTAc                                                 <  2:1531153/150‑1 (MQ=255)
                                                                                                     gggACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTAc                                                 <  1:326435/150‑1 (MQ=255)
                                                                                                      ggACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACg                                                <  2:1684103/150‑1 (MQ=255)
                                                                                                       gACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACgg                                               >  1:1438502/1‑150 (MQ=255)
                                                                                                        accaccTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGc                                              <  1:432588/150‑1 (MQ=255)
                                                                                                             cTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggt                                         <  2:599870/150‑1 (MQ=255)
                                                                                                             cTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggt                                         <  1:1812139/150‑1 (MQ=255)
                                                                                                              tGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggtg                                        <  1:660780/150‑1 (MQ=255)
                                                                                                              tGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggtg                                        >  2:81458/1‑150 (MQ=255)
                                                                                                                gCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTa                                      <  2:2366519/150‑1 (MQ=255)
                                                                                                                 cgtcgtAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTaa                                     <  1:1285713/150‑1 (MQ=255)
                                                                                                                  gtcgtAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAAc                                    <  1:196968/150‑1 (MQ=255)
                                                                                                                  gtcgtAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATACCCGGCAACCACACCCACATCAATTAAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAAc                                    >  2:907334/1‑150 (MQ=255)
                                                                                                                     gtaagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgag                                 >  1:1758583/1‑150 (MQ=255)
                                                                                                                     gtaagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACATGATGCCCGTACGGCGTGGTGTAACgag                                 >  1:1901722/1‑150 (MQ=255)
                                                                                                                       aagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGAc                               >  2:819247/1‑150 (MQ=255)
                                                                                                                       aagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGAc                               <  2:2205906/150‑1 (MQ=255)
                                                                                                                        agtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACa                              <  1:898473/150‑1 (MQ=255)
                                                                                                                        agtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATAGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACa                              >  1:647738/1‑150 (MQ=255)
                                                                                                                                  cTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACgg                    >  2:126686/1‑150 (MQ=255)
                                                                                                                                   tATTGCCGGATGCTAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGgc                   <  1:1695755/150‑1 (MQ=255)
                                                                                                                                       accGGATGCGAAAACTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGcgcgc               >  2:692560/2‑149 (MQ=255)
                                                                                                                                          ggATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACtt            >  1:972755/1‑150 (MQ=255)
                                                                                                                                           gATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACttt           <  2:1042450/150‑1 (MQ=255)
                                                                                                                                              gCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTAt        >  2:377255/1‑150 (MQ=255)
                                                                                                                                               cGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATc       >  2:349495/1‑150 (MQ=255)
                                                                                                                                                  gTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACGTTTATCAcc    >  2:1001245/1‑150 (MQ=255)
                                                                                                                                                   ttttCGCATCCGGTTATCGTCTGCGCCCCCCCCCCATTCCCATCAGTAGCATCCCCGGCCACCACCCCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGCCAATGCAAACGGCGCGACTTTTATCCCCt   >  2:1183533/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCACCTG  >  NC_000913/555872‑556169

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