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

GCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAG  >  NC_000913/3922395‑3922684
                                                                                                                                                       |                                                                                                                                              
gCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTgcgc                                                                                                                                               >  1:482324/1‑151 (MQ=255)
gCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTgcgc                                                                                                                                               >  2:615369/1‑151 (MQ=255)
  aGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGccc                                                                                                                                            >  2:308890/1‑151 (MQ=255)
  aGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGccc                                                                                                                                            >  2:727042/1‑151 (MQ=255)
   gAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGcccc                                                                                                                                           >  1:488160/1‑151 (MQ=255)
      gCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGccccccg                                                                                                                                       >  1:1377019/1‑149 (MQ=255)
          ccGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGAt                                                                                                                                    >  1:110349/1‑151 (MQ=255)
           cGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGAt                                                                                                                                    <  2:1074681/150‑1 (MQ=255)
            gAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGAttt                                                                                                                                  <  1:371981/151‑1 (MQ=255)
              gCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCa                                                                                                                                <  1:833097/151‑1 (MQ=255)
               cAACGTTTCGACGCACGAGCGACACAGTCATCACGTTTTTCACTCCTGCGCCCTTCGAGGTATGCCGCGTCTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAGAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTc                                                                                                                                 <  1:17992/149‑1 (MQ=255)
                  cGTTTCTACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCaaa                                                                                                                              >  1:1245144/1‑149 (MQ=255)
                  cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAAc                                                                                                                             >  2:113596/1‑150 (MQ=255)
                  cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACa                                                                                                                            <  2:493675/151‑1 (MQ=255)
                  cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACa                                                                                                                            <  1:94095/151‑1 (MQ=255)
                     ttCGACTCACGAGCGACACAGACATCACGTTTTTCACTCATGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAAt                                                                                                                          <  2:505394/150‑1 (MQ=255)
                     tccGAGTCACGAGCGACACCGACATCTCTTTTTTCACTTCTGCAACCTTCGACGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTATAGTCAAAGATCAAAAAGCGTACAAATTATACGGTGCATCCCACGTTATTTCAAGCAAt                                                                                                                          <  1:207085/148‑1 (MQ=255)
                      tCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATa                                                                                                                         <  2:453806/150‑1 (MQ=255)
                        gACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGt                                                                                                                      <  1:1166848/151‑1 (MQ=255)
                         aCTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGt                                                                                                                      >  2:470599/1‑150 (MQ=255)
                            cgcGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGc                                                                                                                   <  1:144320/148‑1 (MQ=255)
                            cACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGcc                                                                                                                  >  1:342313/1‑151 (MQ=255)
                            cACGAGCGACACAGACATCACGTTTTTCACGCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGcc                                                                                                                  <  2:815936/151‑1 (MQ=255)
                              cGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCa                                                                                                                 <  1:252362/150‑1 (MQ=255)
                                 gCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACCTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGCCAAATTATACCGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCaaa                                                                                                               >  1:1047737/1‑149 (MQ=255)
                                   gACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGt                                                                                                           >  2:1000364/1‑151 (MQ=255)
                                   gACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGt                                                                                                           <  2:831819/151‑1 (MQ=255)
                                      acaGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGa                                                                                                         >  2:299768/1‑150 (MQ=255)
                                       caGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCACACACTAAGTAGCCAAAAGGTGAAt                                                                                                       >  1:1005234/1‑151 (MQ=255)
                                         gACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaa                                                                                                     <  2:225553/151‑1 (MQ=255)
                                          aCATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaaa                                                                                                    <  2:881304/151‑1 (MQ=255)
                                              cACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGtt                                                                                                <  1:42157/151‑1 (MQ=255)
                                                cGTTTTTCACTCCTGCTCCCTTGGAGGTATCCCGCGTGTCGTATAAATCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATGATACGGTGCGCGCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTGTa                                                                                              <  1:753527/151‑1 (MQ=255)
                                                cGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTa                                                                                              <  2:1187360/151‑1 (MQ=255)
                                                 gTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTaa                                                                                             <  1:701230/151‑1 (MQ=255)
                                                  tttttCACTCCTGCTCCCTTCTAGGTATGCCGCGTGTCGCATAAAACTTTCTTTAAGGCTTAGAGTCATGCATCAAAAAGCGGTCAAATTATACGCTGCGCCCCCCGTGATTTCAAACATTAAGCAGCCAAAAGTTGAATAACTGttt                                                                                               >  1:532086/1‑148 (MQ=255)
                                                    tttCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAat                                                                                           <  1:903115/150‑1 (MQ=255)
                                                       cACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAtt                                                                                        <  1:1166047/150‑1 (MQ=255)
                                                        aCTCCTGCTCCCTTGGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATGCGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAGGGTGAATAAATGTTTAAata                                                                                          <  1:70893/147‑1 (MQ=255)
                                                        aCTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAtttt                                                                                      >  2:929323/1‑151 (MQ=255)
                                                         cTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt                                                                                     >  1:1137744/1‑151 (MQ=255)
                                                           ccTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc                                                                                    <  2:733187/150‑1 (MQ=255)
                                                            cTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt                                                                                     <  2:1268172/148‑1 (MQ=255)
                                                            cTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTcc                                                                                   >  1:1404603/1‑150 (MQ=255)
                                                             tGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTcc                                                                                   >  2:549523/1‑149 (MQ=255)
                                                              gCTCCCTTCGAGGTATGCGGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc                                                                                    <  1:216987/147‑1 (MQ=255)
                                                                 cccTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGc                                                                              <  1:31628/150‑1 (MQ=255)
                                                                     tagAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGAGAATAAATGTTTAAATATTTTTCCAGTGCATAc                                                                          >  1:1184838/3‑150 (MQ=255)
                                                                     tCGAGGTATGCCGCGTGTCGTGTAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATAc                                                                          <  2:1369608/150‑1 (MQ=255)
                                                                     tCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATAc                                                                          >  1:18811/1‑150 (MQ=255)
                                                                       gAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAt                                                                       >  1:789478/1‑151 (MQ=255)
                                                                       gAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAt                                                                       >  1:1406521/1‑151 (MQ=255)
                                                                        aGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAtt                                                                      <  1:1355866/151‑1 (MQ=255)
                                                                         ggTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTg                                                                     >  1:57799/1‑151 (MQ=255)
                                                                          gTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGc                                                                    >  1:505186/1‑151 (MQ=255)
                                                                            aTGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCg                                                                   <  2:1246865/150‑1 (MQ=255)
                                                                               ccGCGTGTCGTATAAAACTATCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCTCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAGTTGCGACtt                                                               <  1:1343598/151‑1 (MQ=255)
                                                                                 gcgTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACtttt                                                             <  2:1212255/151‑1 (MQ=255)
                                                                                   gtgtCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCt                                                           <  2:1159068/151‑1 (MQ=255)
                                                                                       cGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCt                                                        <  2:902634/150‑1 (MQ=255)
                                                                                       cGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCt                                                        >  1:1037663/1‑150 (MQ=255)
                                                                                        gTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGc                                                         <  2:725745/148‑1 (MQ=255)
                                                                                             aaaCTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTg                                                 >  2:875035/1‑151 (MQ=255)
                                                                                                   tctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGa                                            <  2:190156/150‑1 (MQ=255)
                                                                                                    ctttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATc                                          >  2:1178361/1‑151 (MQ=255)
                                                                                                      ttAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGAt                                           <  2:1058552/148‑1 (MQ=255)
                                                                                                       tAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCAt                                        <  1:1010595/150‑1 (MQ=255)
                                                                                                        aaGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAATAGGCGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCAt                                        >  1:430596/1‑149 (MQ=255)
                                                                                                         aGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGa                                      >  1:1419143/1‑150 (MQ=255)
                                                                                                         aGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGa                                      >  1:691150/1‑150 (MQ=255)
                                                                                                                gagTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGt                              <  1:327013/151‑1 (MQ=255)
                                                                                                                 agTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTa                             <  2:10743/151‑1 (MQ=255)
                                                                                                                 agTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTa                             <  2:342325/151‑1 (MQ=255)
                                                                                                                  gTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTa                             >  1:85160/1‑150 (MQ=255)
                                                                                                                    cAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTaaa                           <  2:1377111/150‑1 (MQ=255)
                                                                                                                    cAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAAt                          <  2:211252/151‑1 (MQ=255)
                                                                                                                     aaGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATa                         <  2:537823/151‑1 (MQ=255)
                                                                                                                     aaGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATa                         <  2:1310097/151‑1 (MQ=255)
                                                                                                                      aGCATCGAAAAGCGGTCAGAGTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCGAAAAGGTGAATAAATGTGTAAATATTTTTCCAGTACAGTCAATTGCGACTTTTCTGCTAGCCCTGGGCGGTCATGAAAAGCTGTaa                            <  1:890687/147‑1 (MQ=255)
                                                                                                                      aGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATa                         >  2:1027384/1‑150 (MQ=255)
                                                                                                                       gCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAc                       <  1:108803/151‑1 (MQ=255)
                                                                                                                        cATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcg                      >  2:298872/1‑151 (MQ=255)
                                                                                                                         aTCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcgc                     >  2:707498/1‑151 (MQ=255)
                                                                                                                         aTCAAAAAGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAAcg                      <  1:1275006/150‑1 (MQ=255)
                                                                                                                            aaaaaGCGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTAGATCATGAAAAACTGTAAATAACGCGt                   >  1:681921/1‑150 (MQ=255)
                                                                                                                                  cGGTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACt            <  1:1318904/151‑1 (MQ=255)
                                                                                                                                   ggTCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTg           >  2:445512/1‑151 (MQ=255)
                                                                                                                                     tCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGt         >  2:1051029/1‑151 (MQ=255)
                                                                                                                                     tCAAATTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGt         >  2:828272/1‑151 (MQ=255)
                                                                                                                                       aaaTTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTg        >  2:29689/1‑150 (MQ=255)
                                                                                                                                        aaTTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGAt      >  2:348009/1‑151 (MQ=255)
                                                                                                                                        aaTTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGAt      >  1:571221/1‑151 (MQ=255)
                                                                                                                                        aaTTATACGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGAt      >  2:1305639/1‑151 (MQ=255)
                                                                                                                                            ataCGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATTAAAAACTGTAAATACCGCGTAAACACTGGTGATAAAg  <  2:909294/151‑1 (MQ=255)
                                                                                                                                            ataCGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGGGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACACTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATaa    >  1:706783/1‑149 (MQ=255)
                                                                                                                                             taCGGTGCGCCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTGAATAACGCGTAAACACTGGTGATAAAg  <  1:693536/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                              
GCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCTAACCCTGTTCGATCATGAAAAACTGTAAATAACGCGTAAACACTGGTGATAAAG  >  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: