Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,168,540 T→A intergenic (+93/‑101) gltT → / → rrlB tRNA‑Glu/23S ribosomal RNA of rrnB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,168,5400TA95.4% 231.4 / 4.0 87intergenic (+93/‑101)gltT/rrlBtRNA‑Glu/23S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base A (51/32);  total (54/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168445‑4168636
                                                                                               |                                                                                                
ccACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAg                                                                                              >  2:1588661/1‑100 (MQ=255)
       cTGGTTTGTGAGTTAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGaaaa                                                                                       <  1:1809754/100‑1 (MQ=255)
       cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       >  1:396747/1‑100 (MQ=255)
       cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       >  1:1764931/1‑100 (MQ=255)
         ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAtt                                                                                     >  2:137879/1‑100 (MQ=255)
         ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAtt                                                                                     <  1:367551/100‑1 (MQ=255)
         ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAtt                                                                                     >  2:819859/1‑100 (MQ=255)
          gTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTg                                                                                    <  1:643153/100‑1 (MQ=255)
           tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                   <  1:1660292/100‑1 (MQ=255)
            ttGTTAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  2:1621695/100‑1 (MQ=255)
            ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  2:262971/100‑1 (MQ=255)
            ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  2:1377424/100‑1 (MQ=255)
             tgtgAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaaa                                                                                 <  1:662167/100‑1 (MQ=255)
               tgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAac                                                                                <  1:1855754/99‑1 (MQ=255)
                gagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAacac                                                                              <  1:283434/100‑1 (MQ=255)
                 agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACt                                                                             >  2:71991/1‑100 (MQ=255)
                   tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           <  1:1252245/100‑1 (MQ=255)
                    gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                          >  2:988579/1‑100 (MQ=255)
                    gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                          >  2:689180/1‑100 (MQ=255)
                    gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                          <  2:954818/100‑1 (MQ=255)
                    gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                          >  1:1377017/1‑100 (MQ=255)
                     aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         <  2:1554789/100‑1 (MQ=255)
                      aaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaaca                                                                        >  1:227962/1‑100 (MQ=255)
                      aaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaaca                                                                        <  1:291145/100‑1 (MQ=255)
                          cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    >  2:1791694/1‑100 (MQ=25)
                           gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   >  2:1352200/1‑100 (MQ=25)
                           gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   <  2:1356452/100‑1 (MQ=25)
                           gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   <  2:1836698/100‑1 (MQ=25)
                              gACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgt                                                                >  2:840799/1‑100 (MQ=25)
                                  ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgtt                                                            >  2:727195/1‑100 (MQ=17)
                                      tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        >  1:1523456/1‑100 (MQ=18)
                                      tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  1:995493/100‑1 (MQ=18)
                                      tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  1:1561805/100‑1 (MQ=18)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAAGGAGAGTTGTTCGTGAg                                                      >  1:1136615/1‑100 (MQ=12)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  1:39298/100‑1 (MQ=18)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  2:356013/1‑100 (MQ=18)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  2:917250/1‑100 (MQ=18)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  2:932857/100‑1 (MQ=18)
                                         ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTGGTGAGt                                                     >  2:481933/1‑100 (MQ=12)
                                         ctcAAAACTCACCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                     >  1:1758380/1‑100 (MQ=255)
                                          tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtc                                                    >  1:348632/1‑100 (MQ=32)
                                           cAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtct                                                   >  2:1524341/1‑100 (MQ=33)
                                           cAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtct                                                   >  2:1625228/1‑100 (MQ=18)
                                           cAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtct                                                   >  1:202631/1‑100 (MQ=17)
                                            aaaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctc                                                  <  2:1238754/100‑1 (MQ=17)
                                                 tCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaaa                                             >  1:1222677/1‑100 (MQ=17)
                                                     cTTCGGGTGATTTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAtttt                                         <  2:396749/100‑1 (MQ=255)
                                                     cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAtttt                                         <  2:1764942/100‑1 (MQ=18)
                                                      ttCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTc                                        >  1:317497/1‑100 (MQ=18)
                                                         gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                     >  2:481989/1‑100 (MQ=18)
                                                         gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                     <  2:620234/100‑1 (MQ=18)
                                                         gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                     >  2:1929865/1‑100 (MQ=17)
                                                          ggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCaa                                    >  2:1735003/1‑100 (MQ=18)
                                                            tGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAACCACTGAACAACGAGAGTTGTTCGTGAGTCTCTCCAATCTTcccccc                                   >  1:1561547/1‑94 (MQ=11)
                                                              aTGTTTGAGAGATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGACTTCTTCGTGAGTCTCTCAAATTTTCGCAACACg                                <  1:1279483/100‑1 (MQ=11)
                                                               tGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga                               >  2:1463208/1‑100 (MQ=18)
                                                                gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAGTTTTCGCAACACgat                              >  2:1912839/1‑100 (MQ=255)
                                                                 tttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg                             >  2:1086135/1‑100 (MQ=18)
                                                                 tttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg                             >  2:1390071/1‑100 (MQ=18)
                                                                  ttGTGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg                             <  1:1362572/99‑1 (MQ=12)
                                                                  ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGATGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga                            >  2:1363851/1‑100 (MQ=255)
                                                                  ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGATGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga                            >  2:1364130/1‑100 (MQ=255)
                                                                    gagaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg                          >  1:966267/1‑100 (MQ=18)
                                                                      gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                        >  1:510201/1‑100 (MQ=17)
                                                                        tatTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                      <  2:1808168/100‑1 (MQ=18)
                                                                        tatTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                      <  1:1352192/100‑1 (MQ=17)
                                                                         atTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg                     >  1:1754114/1‑100 (MQ=17)
                                                                            tGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                  >  2:1806714/1‑100 (MQ=18)
                                                                            tGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                  <  2:1377025/100‑1 (MQ=18)
                                                                             gCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaag                 >  1:484120/1‑100 (MQ=18)
                                                                              ctctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaaga                >  2:1652300/1‑100 (MQ=18)
                                                                               tctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa               >  2:7679/1‑100 (MQ=18)
                                                                                ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  2:250156/1‑100 (MQ=18)
                                                                                ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  2:327162/1‑100 (MQ=18)
                                                                                ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  1:185353/1‑100 (MQ=17)
                                                                                 tttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAAc             <  1:1189333/100‑1 (MQ=17)
                                                                                 tttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTAGCAACACGATGATGAATCGAAAGAAAc             >  2:336861/1‑100 (MQ=255)
                                                                                  ttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACa            >  2:1233307/1‑100 (MQ=33)
                                                                                  ttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACa            <  2:1631711/100‑1 (MQ=18)
                                                                                    aaaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATc          >  1:1418132/1‑100 (MQ=17)
                                                                                        aTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCg      <  2:1523464/100‑1 (MQ=18)
                                                                                        aTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCg      <  1:71991/100‑1 (MQ=18)
                                                                                         tCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     >  2:1800987/1‑100 (MQ=18)
                                                                                         tCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     >  2:1127028/1‑100 (MQ=18)
                                                                                            ggAACAAGCTGAAAATTGAAACACTTAACAACGAGAGTTGTTCGTGAGTCTCTCAACTTTTCGCAACACGATGATGAATCGAAAGAAACATCtccagg    >  1:1270172/1‑93 (MQ=11)
                                                                                            ggAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     >  1:1268757/1‑97 (MQ=17)
                                                                                            ggAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCACAACACGATGATGAATCGAAAGAAACATCTTCgggat  >  1:1397203/1‑98 (MQ=255)
                                                                                               |                                                                                                
CCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168445‑4168636

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