Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,145,792 G→A 5.0% A285T (GCA→ACA)  rluC → 23S rRNA pseudouridine(955/2504/2580) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,145,7920GA 5.0% 257.9 / 2.0 79A285T (GCA→ACA) rluC23S rRNA pseudouridine(955/2504/2580) synthase
Reads supporting (aligned to +/- strand):  ref base G (41/34);  new base A (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CCCGGTAACAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGC  >  NC_000913/1145644‑1145933
                                                                                                                                                    |                                                                                                                                             
cccGGTAACAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGc                                                                                                                                              <  1:1276632/150‑1 (MQ=255)
     tAACAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGt                                                                                                                                         <  2:59375/150‑1 (MQ=255)
        cAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGa                                                                                                                                      >  2:9303/1‑150 (MQ=255)
        cAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGa                                                                                                                                      >  1:1172436/1‑150 (MQ=255)
        cAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTACAGCGTTGa                                                                                                                                      >  2:353588/1‑150 (MQ=255)
        cAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTACAGCGTTGa                                                                                                                                      >  2:359177/1‑150 (MQ=255)
          gggCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAg                                                                                                                                    <  1:1272231/150‑1 (MQ=255)
                aCTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGGCCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACt                                                                                                                              <  1:961233/150‑1 (MQ=255)
                   cACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCAt                                                                                                                           >  2:567212/1‑150 (MQ=255)
                       aGATCCGTGTGCATAAACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCgg                                                                                                                       <  2:236736/150‑1 (MQ=255)
                        gATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCggg                                                                                                                      >  2:774324/1‑150 (MQ=255)
                        gATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCggg                                                                                                                      >  2:484960/1‑150 (MQ=255)
                          tCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGa                                                                                                                                      <  2:987511/132‑1 (MQ=255)
                          tCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGa                                                                                                                                      >  1:987511/1‑132 (MQ=255)
                            cGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGAcc                                                                                                                  >  2:804042/1‑150 (MQ=255)
                            cGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGAcc                                                                                                                  <  1:748655/150‑1 (MQ=255)
                                tgCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtg                                                                                                              >  2:586808/1‑150 (MQ=255)
                                 gCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtga                                                                                                             >  2:687247/1‑150 (MQ=255)
                                 gCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtga                                                                                                             >  2:739784/1‑150 (MQ=255)
                                 gCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtga                                                                                                             <  1:1438996/150‑1 (MQ=255)
                                    tACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGGAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtgaggt                                                                                                          <  2:419833/150‑1 (MQ=255)
                                    tACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtgaggt                                                                                                          <  1:1010861/150‑1 (MQ=255)
                                    tACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtgaggt                                                                                                          >  1:235843/1‑150 (MQ=255)
                                      cacaATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCggtgaggtga                                                                                                        >  2:1414055/1‑150 (MQ=255)
                                       acaATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGAt                                                                                                       >  1:1009593/1‑150 (MQ=255)
                                              gCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATc                                                                                                <  1:519348/150‑1 (MQ=255)
                                                gggTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGa                                                                                              >  1:252146/1‑150 (MQ=255)
                                                   tCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGACAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGc                                                                                           <  2:789056/150‑1 (MQ=255)
                                                     aTCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTc                                                                                         <  1:9303/150‑1 (MQ=255)
                                                       ccGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCg                                                                                       <  2:107054/150‑1 (MQ=255)
                                                         gATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCgat                                                                                     <  1:774324/150‑1 (MQ=255)
                                                           ttGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCgatgg                                                                                   <  2:400566/150‑1 (MQ=255)
                                                            tGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCgatgga                                                                                  >  1:296782/1‑150 (MQ=255)
                                                              ccTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCgatggatg                                                                                >  2:363793/1‑150 (MQ=255)
                                                                      gATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTg                                                                        <  2:1111218/150‑1 (MQ=255)
                                                                                gTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCgttgtt                                                              >  1:212223/1‑150 (MQ=255)
                                                                                 tGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGttt                                                             >  1:292380/1‑150 (MQ=255)
                                                                                    ccGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCa                                                          <  2:296782/150‑1 (MQ=255)
                                                                                     cGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCaa                                                         >  1:64450/1‑150 (MQ=255)
                                                                                     cGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCaa                                                         >  2:1118199/1‑150 (MQ=255)
                                                                                       tGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCaaaa                                                       >  2:337173/1‑150 (MQ=255)
                                                                                             tGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCg                                                 >  2:553044/1‑150 (MQ=255)
                                                                                             tGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCg                                                 >  1:930444/1‑150 (MQ=255)
                                                                                                cAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTaa                                              >  1:1021307/1‑150 (MQ=255)
                                                                                                cAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTaa                                              >  2:616417/1‑150 (MQ=255)
                                                                                                 aGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAAc                                             >  2:187166/1‑150 (MQ=255)
                                                                                                  ggCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAAcg                                            >  2:746631/1‑150 (MQ=255)
                                                                                                   gCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAAcgc                                           >  1:1312693/1‑150 (MQ=255)
                                                                                                     aGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAAcgcgc                                         >  1:699552/1‑150 (MQ=255)
                                                                                                      gCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAAcgcgcg                                        <  1:215246/150‑1 (MQ=255)
                                                                                                       cTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACgcgcgc                                       <  1:1513754/150‑1 (MQ=255)
                                                                                                        tCACTGAAGCAGGCACGGGATCAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCt                                      <  2:1041557/150‑1 (MQ=255)
                                                                                                         cACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTa                                     >  2:391078/1‑150 (MQ=255)
                                                                                                           cTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAAt                                   >  2:563571/1‑150 (MQ=255)
                                                                                                           cTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAAt                                   >  1:879972/1‑150 (MQ=255)
                                                                                                           cTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAAt                                   <  2:428984/150‑1 (MQ=255)
                                                                                                             gAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATAc                                 >  1:286286/1‑150 (MQ=255)
                                                                                                               aGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACaa                               <  2:212223/150‑1 (MQ=255)
                                                                                                                  aGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCa                            <  1:586808/150‑1 (MQ=255)
                                                                                                                   ggCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCa                            >  2:1506971/1‑149 (MQ=255)
                                                                                                                   ggCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCa                            <  1:1506971/149‑1 (MQ=255)
                                                                                                                     cACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGGTTGAAGCGTTTTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCatat                         >  1:819565/1‑150 (MQ=255)
                                                                                                                      aCGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCAtata                        >  2:736376/1‑150 (MQ=255)
                                                                                                                       cGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATaa                       >  2:1112610/1‑150 (MQ=255)
                                                                                                                        gggATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATaaa                      >  2:481268/1‑150 (MQ=255)
                                                                                                                         ggATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAAc                     <  1:337173/150‑1 (MQ=255)
                                                                                                                           aTTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCt                   <  1:363793/150‑1 (MQ=255)
                                                                                                                             tAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGa                 >  1:1262978/1‑150 (MQ=255)
                                                                                                                                 tCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAAt             <  1:563571/150‑1 (MQ=255)
                                                                                                                                   gCCTGTTCCTGTACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCGCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATgg           <  2:776886/150‑1 (MQ=255)
                                                                                                                                   gCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATgg           >  1:589402/1‑150 (MQ=255)
                                                                                                                                     cTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGc         <  1:687247/150‑1 (MQ=255)
                                                                                                                                     cTGTTCCTGCACGCTACAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGc         <  1:353588/150‑1 (MQ=255)
                                                                                                                                     cTGTTCCTGCACGCTACAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGc         <  1:359177/150‑1 (MQ=255)
                                                                                                                                          ccTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTg    >  2:911185/1‑150 (MQ=255)
                                                                                                                                           cTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTgg   <  1:553044/150‑1 (MQ=255)
                                                                                                                                           cTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTgg   <  1:391078/150‑1 (MQ=255)
                                                                                                                                            tGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGc  <  2:879972/150‑1 (MQ=255)
                                                                                                                                            cgCACGCTGCAGCGTTGAAGGTTACTCATCCGGGGACCGGTGAGGTGACGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGc  <  2:233680/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CCCGGTAACAGGGCGTACTCACCAGATCCGTGTGCATACACAATATGCGGGTCATCCGATTGCCTTTGACGATCGCTACGGTGACCGTGAATTTGACAGGCAGCTCACTGAAGCAGGCACGGGATTAAATCGCCTGTTCCTGCACGCTGCAGCGTTGAAGTTTACTCATCCGGGGACCGGTGAGGTGATGCGTATCGAAGCTCCGATGGATGAAGGTTTGAAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGC  >  NC_000913/1145644‑1145933

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