Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,922,549 (C)5→6 intergenic (‑109/+508) atpI ← / ← rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,922,5441.C96.7% 308.1 / 2.2 90intergenic (‑104/+513)atpI/rsmGATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (39/48);  total (39/51)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑C‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAG  >  NC_000913/3922395‑3922684
                                                                                                                                                       |                                                                                                                                              
gCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTgc‑g‑c                                                                                                                                               >  1:482324/1‑151 (MQ=255)
gCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTgc‑g‑c                                                                                                                                               >  2:615369/1‑151 (MQ=255)
  aGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑c‑cc                                                                                                                                            >  2:308890/1‑151 (MQ=255)
  aGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑c‑cc                                                                                                                                            >  2:727042/1‑151 (MQ=255)
   gAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑c‑ccc                                                                                                                                           >  1:488160/1‑151 (MQ=255)
      gCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑c‑ccc‑ccg                                                                                                                                       >  1:1377019/1‑149 (MQ=255)
          ccGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGAt                                                                                                                                    >  1:110349/1‑151 (MQ=255)
           cGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGAt                                                                                                                                    <  2:1074681/150‑1 (MQ=255)
            gAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGAttt                                                                                                                                  <  1:371981/151‑1 (MQ=255)
              gCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCa                                                                                                                                <  1:833097/151‑1 (MQ=255)
               cAACGTTTCGACGCACGAGCGACACAGTCATCACGTTTTTCACTCCTGCGCCCTTCGAGGTATGCCGCGTCTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAGAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTc                                                                                                                                 <  1:17992/149‑1 (MQ=255)
                  cGTTTCTACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCaaa                                                                                                                              >  1:1245144/1‑149 (MQ=255)
                  cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAAc                                                                                                                             >  2:113596/1‑150 (MQ=255)
                  cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACa                                                                                                                            <  2:493675/151‑1 (MQ=255)
                  cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACa                                                                                                                            <  1:94095/151‑1 (MQ=255)
                     ttCGACTCACGAGCGACACAGACATCACGTTTTTCACTCATGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAAt                                                                                                                          <  2:505394/150‑1 (MQ=255)
                     tccGAGTCACGAGCGACACCGACATCTCTTTTTTCACTTCTGCAACCTTCGACGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTATAGTCAAAGATCAAAAAGCGTACAAATTATACGGTGCAT‑C‑CCA‑CGTTATTTCAAGCAAt                                                                                                                          <  1:207085/148‑1 (MQ=255)
                      tCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATa                                                                                                                         <  2:453806/150‑1 (MQ=255)
                        gACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGt                                                                                                                      <  1:1166848/151‑1 (MQ=255)
                         aCTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGt                                                                                                                      >  2:470599/1‑150 (MQ=255)
                            cgcGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGc                                                                                                                   <  1:144320/148‑1 (MQ=255)
                            cACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGcc                                                                                                                  >  1:342313/1‑151 (MQ=255)
                            cACGAGCGACACAGACATCACGTTTTTCACGCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGcc                                                                                                                  <  2:815936/151‑1 (MQ=255)
                              cGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCa                                                                                                                 <  1:252362/150‑1 (MQ=255)
                                 gCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACCTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGCCAAATTATACCGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCaaa                                                                                                               >  1:1047737/1‑149 (MQ=255)
                                   gACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGt                                                                                                           >  2:1000364/1‑151 (MQ=255)
                                   gACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGt                                                                                                           <  2:831819/151‑1 (MQ=255)
                                      acaGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGa                                                                                                         >  2:299768/1‑150 (MQ=255)
                                       caGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCACACACTAAGTAGCCAAAAGGTGAAt                                                                                                       >  1:1005234/1‑151 (MQ=255)
                                         gACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaa                                                                                                     <  2:225553/151‑1 (MQ=255)
                                          aCATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaaa                                                                                                    <  2:881304/151‑1 (MQ=255)
                                              cACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGtt                                                                                                <  1:42157/151‑1 (MQ=255)
                                                cGTTTTTCACTCCTGCTCCCTTGGAGGTATCCCGCGTGTCGTATAAATCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATGATACGGTGC‑G‑CGCCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTGTa                                                                                              <  1:753527/151‑1 (MQ=255)
                                                cGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTa                                                                                              <  2:1187360/151‑1 (MQ=255)
                                                 gTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTaa                                                                                             <  1:701230/151‑1 (MQ=255)
                                                  tttttCACTCCTGCTCCCTTCTAGGTATGCCGCGTGTCGCATAAAACTTTCTTTAAGGCTTAGAGTCATGCATCAAAAAGCGGTCAAATTATACGCTGC‑GCC‑CCC‑CGTGATTTCAAACATTAAGCAGCCAAAAGTTGAATAACTGttt                                                                                               >  1:532086/1‑148 (MQ=255)
                                                    tttCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAat                                                                                           <  1:903115/150‑1 (MQ=255)
                                                       cACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAtt                                                                                        <  1:1166047/150‑1 (MQ=255)
                                                        aCTCCTGCTCCCTTGGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATGCGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAGGGTGAATAAATGTTTAAata                                                                                          <  1:70893/147‑1 (MQ=255)
                                                        aCTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAtttt                                                                                      >  2:929323/1‑151 (MQ=255)
                                                         cTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt                                                                                     >  1:1137744/1‑151 (MQ=255)
                                                           ccTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc                                                                                    <  2:733187/150‑1 (MQ=255)
                                                            cTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt                                                                                     <  2:1268172/148‑1 (MQ=255)
                                                            cTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTcc                                                                                   >  1:1404603/1‑150 (MQ=255)
                                                             tGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTcc                                                                                   >  2:549523/1‑149 (MQ=255)
                                                              gCTCCCTTCGAGGTATGCGGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc                                                                                    <  1:216987/147‑1 (MQ=255)
                                                                 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGc                                                                              <  1:31628/150‑1 (MQ=255)
                                                                     tagAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGAGAATAAATGTTTAAATATTTTTCCAGTGCATAc                                                                          >  1:1184838/3‑150 (MQ=255)
                                                                     tCGAGGTATGCCGCGTGTCGTGTAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATAc                                                                          <  2:1369608/150‑1 (MQ=255)
                                                                     tCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATAc                                                                          >  1:18811/1‑150 (MQ=255)
                                                                       gAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAt                                                                       >  1:789478/1‑151 (MQ=255)
                                                                       gAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAt                                                                       >  1:1406521/1‑151 (MQ=255)
                                                                        aGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAtt                                                                      <  1:1355866/151‑1 (MQ=255)
                                                                         ggTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTg                                                                     >  1:57799/1‑151 (MQ=255)
                                                                          gTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGc                                                                    >  1:505186/1‑151 (MQ=255)
                                                                            aTGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCg                                                                   <  2:1246865/150‑1 (MQ=255)
                                                                               ccGCGTGTCGTATAAAACTATCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑C‑CCCTCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAGTTGCGACtt                                                               <  1:1343598/151‑1 (MQ=255)
                                                                                 gcgTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtttt                                                             <  2:1212255/151‑1 (MQ=255)
                                                                                   gtgtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCt                                                           <  2:1159068/151‑1 (MQ=255)
                                                                                       cGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCt                                                        <  2:902634/150‑1 (MQ=255)
                                                                                       cGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCt                                                        >  1:1037663/1‑150 (MQ=255)
                                                                                        gTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGc                                                         <  2:725745/148‑1 (MQ=255)
                                                                                             aaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTg                                                 >  2:875035/1‑151 (MQ=255)
                                                                                                   tctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGa                                            <  2:190156/150‑1 (MQ=255)
                                                                                                    ctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATc                                          >  2:1178361/1‑151 (MQ=255)
                                                                                                      ttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGAt                                           <  2:1058552/148‑1 (MQ=255)
                                                                                                       tAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCAt                                        <  1:1010595/150‑1 (MQ=255)
                                                                                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAATAGGCGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCAt                                        >  1:430596/1‑149 (MQ=255)
                                                                                                         aGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGa                                      >  1:1419143/1‑150 (MQ=255)
                                                                                                         aGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGa                                      >  1:691150/1‑150 (MQ=255)
                                                                                                                gagTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGt                              <  1:327013/151‑1 (MQ=255)
                                                                                                                 agTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTa                             <  2:10743/151‑1 (MQ=255)
                                                                                                                 agTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTa                             <  2:342325/151‑1 (MQ=255)
                                                                                                                  gTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTa                             >  1:85160/1‑150 (MQ=255)
                                                                                                                    cAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTaaa                           <  2:1377111/150‑1 (MQ=255)
                                                                                                                    cAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAAt                          <  2:211252/151‑1 (MQ=255)
                                                                                                                     aaGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATa                         <  2:537823/151‑1 (MQ=255)
                                                                                                                     aaGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATa                         <  2:1310097/151‑1 (MQ=255)
                                                                                                                      aGCATCGAAAAGCGGTCAGAGTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCGAAAAGGTGAATAAATGTGTAAATATTTTTCCAGTACAGTCAATTGCGACTTTTCTGCTAGCCCTGGGCGGTCATGAAAAGCTGTaa                            <  1:890687/147‑1 (MQ=255)
                                                                                                                      aGCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATa                         >  2:1027384/1‑150 (MQ=255)
                                                                                                                       gCATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAc                       <  1:108803/151‑1 (MQ=255)
                                                                                                                        cATCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcg                      >  2:298872/1‑151 (MQ=255)
                                                                                                                         aTCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcgc                     >  2:707498/1‑151 (MQ=255)
                                                                                                                         aTCAAAAAGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcg                      <  1:1275006/150‑1 (MQ=255)
                                                                                                                            aaaaaGCGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTAGATCATGAAAAACTGTAAATAACGCGt                   >  1:681921/1‑150 (MQ=255)
                                                                                                                                  cGGTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACt            <  1:1318904/151‑1 (MQ=255)
                                                                                                                                   ggTCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTg           >  2:445512/1‑151 (MQ=255)
                                                                                                                                     tCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGt         >  2:1051029/1‑151 (MQ=255)
                                                                                                                                     tCAAATTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGt         >  2:828272/1‑151 (MQ=255)
                                                                                                                                       aaaTTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTg        >  2:29689/1‑150 (MQ=255)
                                                                                                                                        aaTTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGAt      >  2:348009/1‑151 (MQ=255)
                                                                                                                                        aaTTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGAt      >  1:571221/1‑151 (MQ=255)
                                                                                                                                        aaTTATACGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGAt      >  2:1305639/1‑151 (MQ=255)
                                                                                                                                            ataCGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATTAAAAACTGTAAATACCGCGTAAACACTGGTGATAAAg  <  2:909294/151‑1 (MQ=255)
                                                                                                                                            ataCGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGGGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACACTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATaa    >  1:706783/1‑149 (MQ=255)
                                                                                                                                             taCGGTGC‑GCC‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTGAATAACGCGTAAACACTGGTGATAAAg  <  1:693536/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                              
GCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGC‑G‑C‑CCC‑CGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAG  >  NC_000913/3922395‑3922684

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