Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,019 T→A 5.3% intergenic (+23/+13) cysS → / ← ybcI cysteinyl‑tRNA synthetase/DUF457 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,0190TA 5.3% 657.3 / 18.1 184intergenic (+23/+13)cysS/ybcIcysteinyl‑tRNA synthetase/DUF457 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (69/105);  new base A (2/8);  total (71/113)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

AACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCACCT  >  NC_000913/555872‑556168
                                                                                                                                                   |                                                                                                                                                     
aaCAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGttt                                                                                                                                                     >  2:1597167/1‑150 (MQ=255)
aaCAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGttt                                                                                                                                                     >  1:1845174/1‑150 (MQ=255)
 aCAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGtttt                                                                                                                                                    >  1:3227074/1‑150 (MQ=255)
 aCAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGtttt                                                                                                                                                    >  1:3227103/1‑150 (MQ=255)
  cAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTc                                                                                                                                                   <  1:2586195/150‑1 (MQ=255)
  cAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTc                                                                                                                                                   <  1:2070956/150‑1 (MQ=255)
   aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCg                                                                                                                                                  >  1:2452092/1‑150 (MQ=255)
   aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCg                                                                                                                                                  <  2:808765/150‑1 (MQ=255)
   aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCg                                                                                                                                                  >  1:767623/1‑150 (MQ=255)
      gTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCAt                                                                                                                                               >  1:366409/1‑150 (MQ=255)
       tCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATc                                                                                                                                              <  2:221293/150‑1 (MQ=255)
       tCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATc                                                                                                                                              <  1:2917453/150‑1 (MQ=255)
        cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTACCGGATGCGAAAACTCGCATcc                                                                                                                                             <  1:227153/150‑1 (MQ=255)
         tGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCg                                                                                                                                            <  1:2284877/150‑1 (MQ=255)
         tGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATCCg                                                                                                                                            <  1:2554054/150‑1 (MQ=255)
          ggATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATCCgg                                                                                                                                           <  2:1705461/150‑1 (MQ=255)
           gATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGt                                                                                                                                          <  1:43410/150‑1 (MQ=255)
               cccGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATc                                                                                                                                      >  2:1606288/1‑150 (MQ=255)
                ccGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCg                                                                                                                                     >  2:720512/1‑150 (MQ=255)
                  gTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTc                                                                                                                                   >  2:2407221/1‑150 (MQ=255)
                  gTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTc                                                                                                                                   >  2:2819201/1‑150 (MQ=255)
                              aCTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACcaca                                                                                                                       <  1:2521245/150‑1 (MQ=255)
                              aCTGGGCGGCGGCGGATGCGGCGCGTGATAGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTTAGTAATTGCGATATTTCCGGATGCGAGTTTTCGCATCCGTTTATAGTCTGCGCCACcaca                                                                                                                       <  2:1446658/150‑1 (MQ=255)
                                         gcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCacaa                                                                                                                      >  1:1279679/1‑140 (MQ=255)
                                         gcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATc                                                                                                            <  2:1890954/150‑1 (MQ=255)
                                          cggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCa                                                                                                           <  2:1167747/150‑1 (MQ=255)
                                           ggctGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGGCATTGCGCTATTGCCGGATGCGCGTTTTCGCATCCGGGTATCGGCTGCGCCACCACAACAGTCCCATccg                                                                                                          >  2:870321/4‑148 (MQ=255)
                                             aTGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTa                                                                                                        <  2:1834279/150‑1 (MQ=255)
                                             aTGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTa                                                                                                        <  2:1416411/150‑1 (MQ=255)
                                              tGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAg                                                                                                       <  2:1706186/150‑1 (MQ=255)
                                              tGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAg                                                                                                       <  2:2697193/150‑1 (MQ=255)
                                                cggcgCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCa                                                                                                     <  2:2850782/150‑1 (MQ=255)
                                                 ggcgCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTCATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCAt                                                                                                    >  2:2282213/1‑150 (MQ=255)
                                                     cgTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATcccc                                                                                                <  2:558270/150‑1 (MQ=255)
                                                      gTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCCATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCg                                                                                               <  1:2093539/150‑1 (MQ=255)
                                                       tGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCgg                                                                                              <  1:2130254/150‑1 (MQ=255)
                                                       tGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCgg                                                                                              <  1:1763752/150‑1 (MQ=255)
                                                       tGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCgg                                                                                              <  2:1596265/150‑1 (MQ=255)
                                                       tGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCgg                                                                                              <  1:469267/150‑1 (MQ=255)
                                                          tCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCaa                                                                                           <  2:734636/150‑1 (MQ=255)
                                                          tCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCaa                                                                                           <  1:2324014/150‑1 (MQ=255)
                                                           cGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAc                                                                                          <  1:90371/150‑1 (MQ=255)
                                                               ttAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGATTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACcaca                                                                                      <  2:991664/150‑1 (MQ=255)
                                                                      gaTGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACAt                                                                               <  1:573189/150‑1 (MQ=255)
                                                                       aTGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATc                                                                              >  1:1220775/1‑150 (MQ=255)
                                                                         ggggATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTAGTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCac                                                                            >  2:2103051/1‑149 (MQ=255)
                                                                          gggATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCCCATCAAt                                                                           >  2:2774235/1‑150 (MQ=255)
                                                                            gATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGTGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTc                                                                         >  1:3239545/1‑150 (MQ=255)
                                                                            gATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTc                                                                         >  2:887837/1‑150 (MQ=255)
                                                                             aTCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCa                                                                        <  1:2434046/150‑1 (MQ=255)
                                                                             aTCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCa                                                                        <  1:666658/150‑1 (MQ=255)
                                                                             aTCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCa                                                                        <  1:2866667/150‑1 (MQ=255)
                                                                               cGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGa                                                                      <  1:1214549/150‑1 (MQ=255)
                                                                               cGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGa                                                                      <  2:1086129/150‑1 (MQ=255)
                                                                                gTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGaa                                                                     <  2:2260151/150‑1 (MQ=255)
                                                                                gTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTAGTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGaa                                                                     <  2:3074434/150‑1 (MQ=255)
                                                                                  gCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaat                                                                   >  1:2597861/1‑150 (MQ=255)
                                                                                  gCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaat                                                                   >  1:1183382/1‑150 (MQ=255)
                                                                                  gCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaat                                                                   >  1:1012005/1‑150 (MQ=255)
                                                                                    tGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaataa                                                                 >  1:276611/1‑150 (MQ=255)
                                                                                         gATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCAcc                                                            >  1:1027000/1‑150 (MQ=255)
                                                                                           tGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTg                                                          >  2:2942799/1‑150 (MQ=255)
                                                                                              cccGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATg                                                       >  1:1363611/1‑150 (MQ=255)
                                                                                               ccGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGc                                                      <  2:2621711/150‑1 (MQ=255)
                                                                                               ccGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGc                                                      <  2:2640191/150‑1 (MQ=255)
                                                                                               ccGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGc                                                      <  2:637915/150‑1 (MQ=255)
                                                                                                cGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGcc                                                     >  2:2568720/1‑150 (MQ=255)
                                                                                                cGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGcc                                                     >  1:3134754/1‑150 (MQ=255)
                                                                                                        accaccTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGc                                             <  2:2346289/150‑1 (MQ=255)
                                                                                                          caccTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgt                                           <  1:1749376/150‑1 (MQ=255)
                                                                                                            ccTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtgg                                         >  2:2879537/1‑150 (MQ=255)
                                                                                                             cTGGCGTCGTAAGTAATTTCGCTATTGCCGGATGCGAGTTTTAGCATCCGGTTATAGTCTGCGCCACCACAACATTCCCATCAGTAGCAACCCCGGCAACCACACCCACAAAAAATCAGAAATAATAACCTGATGCCCGTACGGCgtggt                                        >  1:1503209/1‑150 (MQ=255)
                                                                                                                gCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTa                                     >  1:493181/1‑150 (MQ=255)
                                                                                                                gCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTa                                     <  2:1536695/150‑1 (MQ=255)
                                                                                                                  gtcgtAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACCTCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAAc                                   <  2:1224777/150‑1 (MQ=255)
                                                                                                                  gtcgtAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAAc                                   <  2:1376406/150‑1 (MQ=255)
                                                                                                                    cgtAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACga                                 <  1:386921/150‑1 (MQ=255)
                                                                                                                     gtaagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgag                                >  1:3196256/1‑150 (MQ=255)
                                                                                                                     gtaagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgag                                >  1:853729/1‑150 (MQ=255)
                                                                                                                      taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                               >  1:1113501/1‑150 (MQ=255)
                                                                                                                      taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                               >  2:1455340/1‑150 (MQ=255)
                                                                                                                      taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                               >  2:1036478/1‑150 (MQ=255)
                                                                                                                      taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                               <  1:938787/150‑1 (MQ=255)
                                                                                                                         gtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACaa                            <  1:2508937/150‑1 (MQ=255)
                                                                                                                          taaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAAt                           <  1:2973412/150‑1 (MQ=255)
                                                                                                                          taaGTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACCTTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAAt                           <  2:2594764/150‑1 (MQ=255)
                                                                                                                               gcgcTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCaaa                      <  2:2511329/150‑1 (MQ=255)
                                                                                                                                cgcTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAAc                     <  1:1975514/150‑1 (MQ=255)
                                                                                                                                     ttGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGgcgc                <  2:2053500/150‑1 (MQ=255)
                                                                                                                                     ttGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGgcgc                <  2:411116/150‑1 (MQ=255)
                                                                                                                                     ttGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGgcgc                <  2:1175061/150‑1 (MQ=255)
                                                                                                                                         cGGATGCGAAAACTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGACATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACt            >  2:196432/1‑150 (MQ=255)
                                                                                                                                            aTGCTAGTTTTCGCATCCGGTTTTCGTCTGCGCCACCACAACATTCCAATCAGTAGCAACCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGTTGTAACGAGACAATGCAAACGGCGCGACtttt         <  1:1711641/150‑1 (MQ=255)
                                                                                                                                                gtgTTTTCGCATCCGGTTATCGTCTGCGCCACCAAAACATGCCAATCAGTAGCATCCCCGGAAACCACACCCACATCAATTCACAAATAATCACCTGATGCCCGTACGGCGTGGTTTAACGAGACAATGCAAACGGCGCGACTTTTATCa     <  1:1379532/148‑1 (MQ=255)
                                                                                                                                                gAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTTTAACGAGACAATGCAAACGGCGCGACTTTTATCa     >  1:2241589/1‑150 (MQ=255)
                                                                                                                                                 tgTTTTCGCATCCGGTTATCGTCTGCGCCACCACACCATTCCCATCAGTAGCATCCCCGGCAACAACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCAc    <  1:774576/149‑1 (MQ=255)
                                                                                                                                                 aGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATATGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCAc    <  1:2955353/150‑1 (MQ=255)
                                                                                                                                                  gTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCAcc   >  1:2306489/1‑150 (MQ=255)
                                                                                                                                                  gTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCAcc   >  1:2133341/1‑150 (MQ=255)
                                                                                                                                                  gTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCAcc   >  1:2988654/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTTATCACCT  >  NC_000913/555872‑556168

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