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,5400TA97.8% 242.8 / ‑2.5 89intergenic (+93/‑101)gltT/rrlBtRNA‑Glu/23S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (47/40);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168443‑4168636
                                                                                                 |                                                                                                
cGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACa                                                                                                >  2:1487569/1‑100 (MQ=255)
 gCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACaa                                                                                               >  2:721367/1‑100 (MQ=255)
 gCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGAAGTTTGAGATACTTGCTCTTTAAAAATCTGGAACa                                                                                                >  1:481509/1‑99 (MQ=255)
  ccACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAg                                                                                              >  1:768127/1‑100 (MQ=255)
  ccACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAg                                                                                              >  1:880777/1‑100 (MQ=255)
   cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                             <  2:98760/100‑1 (MQ=255)
   cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                             >  1:249219/1‑100 (MQ=255)
   cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                             <  2:392692/100‑1 (MQ=255)
     cTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTg                                                                                           <  2:725606/100‑1 (MQ=255)
       tGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaa                                                                                         >  2:643471/1‑100 (MQ=255)
       tGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCGCTTTAAAAATCTGGAACAAGCTGaa                                                                                         >  2:1538510/1‑100 (MQ=255)
        gCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaa                                                                                        >  1:1366216/1‑100 (MQ=255)
         cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       >  2:793809/1‑100 (MQ=255)
          tGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAt                                                                                      <  1:573575/100‑1 (MQ=255)
          tGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAt                                                                                      <  1:332813/100‑1 (MQ=255)
          tGGGTTGTGAGTGAAAGTGGGTGATCTTAATGTCTCAAAACTCATCTTTGGGGGATGTTTGAGATATTTGCTCTTTAAAAATCTGGGACAAGCTg                                                                                           <  2:874193/95‑1 (MQ=14)
           ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAtt                                                                                     >  1:720171/1‑100 (MQ=255)
           ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAtt                                                                                     <  1:713994/100‑1 (MQ=255)
            gTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTg                                                                                    <  1:310324/100‑1 (MQ=255)
             tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                   <  2:1366209/100‑1 (MQ=255)
               tgtgAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaaa                                                                                 <  1:1024416/100‑1 (MQ=255)
                gtgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAGATCTGGAACAAGCTGAAAATTGAAac                                                                                <  1:964021/100‑1 (MQ=255)
                gtgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCGCTTTAAAAATCTGGATCAAGCTTAAAATTGAAac                                                                                <  2:321932/100‑1 (MQ=255)
                  gagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAacac                                                                              >  1:727080/1‑100 (MQ=255)
                  gagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAacac                                                                              >  1:147724/1‑100 (MQ=255)
                   agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAacac                                                                              <  1:992472/99‑1 (MQ=255)
                   agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACt                                                                             <  2:231300/100‑1 (MQ=255)
                   agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACt                                                                             >  2:1452417/1‑100 (MQ=255)
                    gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            >  1:317634/1‑100 (MQ=255)
                    gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            <  1:1487578/100‑1 (MQ=255)
                    gtgaAAGTCGCCGACCTGAATATCTTAATACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            <  2:1324078/100‑1 (MQ=14)
                     tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           <  1:838795/100‑1 (MQ=255)
                       aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGCAAATTGAAACACTGaac                                                                         <  2:393947/100‑1 (MQ=255)
                       aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         <  1:74297/100‑1 (MQ=255)
                         aGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaacaa                                                                       >  2:418618/1‑100 (MQ=255)
                         aGTCGCCGACCTTAATATCTCAAAACCCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaacaa                                                                       <  1:44561/100‑1 (MQ=255)
                            cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    >  2:1103315/1‑100 (MQ=25)
                             gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   >  1:855545/1‑100 (MQ=25)
                             gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   <  1:536952/100‑1 (MQ=25)
                             gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   >  1:1220169/1‑100 (MQ=25)
                             gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   >  1:178501/1‑100 (MQ=25)
                             gCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgag                                                                   <  1:1174617/100‑1 (MQ=25)
                                   cTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgt                                                             <  1:1163199/100‑1 (MQ=17)
                                    ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgtt                                                            >  1:639793/1‑100 (MQ=18)
                                        tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  2:880772/100‑1 (MQ=18)
                                          tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAATACTGAACAACGAGAGTTGTTCg                                                          >  2:858081/1‑96 (MQ=12)
                                          tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  2:508341/100‑1 (MQ=18)
                                           ctcAAAACTCATCTTCGGGTGTTGTTTGAGATATTTGCTCTTTAAAACTCTGGCACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  2:854412/1‑99 (MQ=11)
                                           ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                     <  1:1374109/100‑1 (MQ=18)
                                           ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTCGTTCGTGAGt                                                     >  2:1035183/1‑100 (MQ=255)
                                              aaaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTACAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctc                                                  >  2:865041/1‑100 (MQ=255)
                                               aaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctct                                                 >  2:1018676/1‑100 (MQ=18)
                                                   tCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaaa                                             <  2:1536721/100‑1 (MQ=17)
                                                    cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAt                                            >  2:332934/1‑100 (MQ=17)
                                                    cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAt                                            >  2:82940/1‑100 (MQ=17)
                                                      tCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAttt                                          <  1:643473/100‑1 (MQ=17)
                                                          cGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGc                                      <  1:979749/100‑1 (MQ=18)
                                                           gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                     >  2:491387/1‑100 (MQ=18)
                                                           gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                     <  1:1126577/100‑1 (MQ=18)
                                                            ggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCaa                                    >  2:1282246/1‑100 (MQ=17)
                                                             gTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAac                                   >  1:1413896/1‑100 (MQ=18)
                                                              tGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                                  <  1:882285/100‑1 (MQ=18)
                                                              tGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                                  <  1:865046/100‑1 (MQ=18)
                                                                aTGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACg                                <  2:123529/100‑1 (MQ=18)
                                                                 tGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga                               >  1:1056185/1‑100 (MQ=18)
                                                                 tGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga                               >  1:432219/1‑100 (MQ=18)
                                                                  gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  1:68946/1‑100 (MQ=18)
                                                                  gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  2:359639/1‑100 (MQ=18)
                                                                   tttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatg                             <  2:534143/100‑1 (MQ=18)
                                                                    ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga                            >  1:1424985/1‑100 (MQ=18)
                                                                    ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  1:986312/1‑98 (MQ=18)
                                                                      gagaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGatgatg                          <  1:1100231/100‑1 (MQ=18)
                                                                        gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                        >  2:193570/1‑100 (MQ=18)
                                                                          tatTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                      >  2:1551586/1‑100 (MQ=17)
                                                                            tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTGGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGCATCga                    <  2:481508/100‑1 (MQ=255)
                                                                             ttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaa                   >  1:8625/1‑100 (MQ=18)
                                                                                 tctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa               <  2:146046/100‑1 (MQ=18)
                                                                                 tctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaaga                <  2:720169/99‑1 (MQ=17)
                                                                                  ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  1:656729/1‑100 (MQ=18)
                                                                                  ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  2:1137070/1‑100 (MQ=18)
                                                                                   tttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACCACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaatc             <  2:111749/100‑3 (MQ=255)
                                                                                     tAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACAt           >  1:1400957/1‑100 (MQ=17)
                                                                                       aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt         >  2:1464398/1‑100 (MQ=17)
                                                                                       aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt         >  1:1014807/1‑100 (MQ=18)
                                                                                           tCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     <  2:412023/100‑1 (MQ=32)
                                                                                           tCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     <  1:1494431/100‑1 (MQ=18)
                                                                                             tGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGt   <  2:1144565/100‑1 (MQ=18)
                                                                                             tGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTAGCAACGCGATGATGGATCGAAAGAAACATCTTCGGGt   >  2:1064820/1‑100 (MQ=11)
                                                                                              ggAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGtt  >  1:869138/1‑100 (MQ=18)
                                                                                                 |                                                                                                
CGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168443‑4168636

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