Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,168,540 T→A 100% 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.5% 500.3 / 1.5 196intergenic (+93/‑101)gltT/rrlBtRNA‑Glu/23S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base A (93/96);  total (96/100)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 196 total aligned reads displayed.

CGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGT  >  NC_000913/4168443‑4168635
                                                                                                 |                                                                                               
   tcttccgatctttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                            <  2:2984510/90‑1 (MQ=14)
   cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAgn                                                                                            <  2:1012347/100‑2 (MQ=255)
    aCTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCt                                                                                           <  2:1257846/100‑1 (MQ=255)
        gCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaa                                                                                       >  2:651681/1‑100 (MQ=255)
         cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                      <  1:1261366/100‑1 (MQ=255)
         cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                      <  1:3772981/100‑1 (MQ=255)
         cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                      >  1:2226776/1‑100 (MQ=255)
           ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTCATGTTTGAGGTATTTACTCTTTAAGAATCTGGAACAAGCTGAAAAtt                                                                                    >  2:2810544/1‑100 (MQ=14)
            gTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTg                                                                                   <  1:443180/100‑1 (MQ=40)
            gTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTg                                                                                   <  2:3585614/100‑1 (MQ=255)
             tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                  <  2:968279/100‑1 (MQ=255)
             tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                  >  1:3297318/1‑100 (MQ=255)
             tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                  <  1:157737/100‑1 (MQ=255)
              ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGGTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                 >  1:687327/1‑100 (MQ=255)
              ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                 <  1:4045466/100‑1 (MQ=255)
               tgtgAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaaa                                                                                <  2:3732280/100‑1 (MQ=255)
               tgtgAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaaa                                                                                <  2:2160032/100‑1 (MQ=255)
                gtgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAac                                                                               <  1:509115/100‑1 (MQ=255)
                 tgagtgaACGTCGCCGACCTTAATATCTCAAAACTCAGCTGGGGGTGATGTTTTAGATATTTGCACTTTGAAAATCTGGAACAAGCTGAAAATTGAAaca                                                                              >  2:2064849/1‑100 (MQ=14)
                 tgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAaca                                                                              >  1:3115298/1‑100 (MQ=255)
                  gagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAatac                                                                             >  2:3702991/1‑100 (MQ=255)
                  gagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAacac                                                                             <  1:184236/100‑1 (MQ=255)
                   agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGATTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACt                                                                            >  2:4156089/1‑100 (MQ=255)
                    gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                           >  2:2847277/1‑100 (MQ=255)
                    gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                           >  2:701330/1‑100 (MQ=255)
                    gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                           <  1:1097785/100‑1 (MQ=255)
                    gtcaaaGTCGCCGACCTTAATATCTCAACCCTCATCTTCGGGTCATCTTGCAGATATTTGCCCTTTAAAAATCTGGAAAAAGCTGAATATTGCAACACtc                                                                           <  2:3430460/97‑2 (MQ=14)
                     tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                          >  1:3350043/1‑100 (MQ=255)
                      gaAAGTCGCCGATCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                         <  2:822684/100‑1 (MQ=255)
                      gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                         >  1:3926528/1‑100 (MQ=255)
                        aaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaaca                                                                       >  2:4169244/1‑100 (MQ=255)
                         aGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaacaa                                                                      >  2:3488261/1‑100 (MQ=255)
                            cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                   >  2:501300/1‑100 (MQ=25)
                            cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAAAAACga                                                                   >  2:4086226/1‑100 (MQ=255)
                              ccGACCTTAATATCTCAAAACTCATCTTCGGGTGATTTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgaga                                                                 <  1:2923558/100‑1 (MQ=255)
                               cGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGagag                                                                <  1:958126/100‑1 (MQ=25)
                                gACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgt                                                               >  1:3711182/1‑100 (MQ=25)
                                  ccTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttg                                                             >  1:2581235/1‑100 (MQ=17)
                                    ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTCGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgtt                                                           >  1:1204468/1‑100 (MQ=12)
                                     tAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTACACATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTc                                                          >  2:1970278/1‑100 (MQ=255)
                                     tAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTc                                                          >  2:1555564/1‑100 (MQ=17)
                                     tAATATCTCAAAACTCATCTACGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgtt                                                           >  2:2546539/1‑99 (MQ=12)
                                        tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                       <  2:1487274/100‑1 (MQ=18)
                                        tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                       <  1:3971465/100‑1 (MQ=18)
                                        tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                       <  1:2987317/100‑1 (MQ=17)
                                          tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                     <  1:696620/100‑1 (MQ=18)
                                          tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                     <  1:3488287/100‑1 (MQ=18)
                                          tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                     >  2:2783526/1‑100 (MQ=17)
                                           ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                    >  1:1798066/1‑100 (MQ=18)
                                           ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                    >  1:2584320/1‑100 (MQ=18)
                                            tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGATAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtc                                                   <  2:2226762/100‑1 (MQ=255)
                                              aaaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctc                                                 >  1:4115721/1‑100 (MQ=18)
                                               aaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTACAAATCTGGAACAAGCTGAACATTGAAACACTGAACAACGAGCGTTGTTCGTGAGtctct                                                <  1:3911663/100‑1 (MQ=11)
                                                 aCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAACTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCa                                              <  2:1229912/100‑1 (MQ=12)
                                                 aCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCa                                              <  1:886809/100‑1 (MQ=18)
                                                  cTCATCTTTGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGGTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                             <  1:3438617/100‑1 (MQ=255)
                                                  cTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                             >  1:2729637/1‑100 (MQ=18)
                                                       cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAtttt                                        >  2:1132066/1‑100 (MQ=18)
                                                        ttCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTc                                       <  2:1378244/100‑1 (MQ=18)
                                                         tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                      <  1:234709/100‑1 (MQ=17)
                                                         tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                      <  2:954667/100‑1 (MQ=17)
                                                         tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                      <  1:33615/100‑1 (MQ=18)
                                                         tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                      <  1:2499555/100‑1 (MQ=18)
                                                          cGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGc                                     <  1:2871812/100‑1 (MQ=18)
                                                           gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTcaca                                    >  2:88042/1‑100 (MQ=12)
                                                           gggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                    <  2:1786535/100‑1 (MQ=18)
                                                            ggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCaa                                   >  2:379868/1‑100 (MQ=18)
                                                            ggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCaa                                   >  2:4044889/1‑100 (MQ=17)
                                                            ggTGATGTTTGAGATATTTGCTCTCTAAAAGTCTGGAACAAGCTGAAAATTGCAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCa                                    >  2:2682367/1‑99 (MQ=11)
                                                                  gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                             >  1:3678680/1‑100 (MQ=18)
                                                                    ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga                           <  2:1263452/100‑1 (MQ=18)
                                                                     tGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat                          >  2:1634407/1‑100 (MQ=18)
                                                                        gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACATAATGATGaa                       >  2:2475723/1‑100 (MQ=255)
                                                                        gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                       <  1:1973485/100‑1 (MQ=17)
                                                                        gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                       >  2:2628341/1‑100 (MQ=18)
                                                                          tatTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                     <  2:2091324/100‑1 (MQ=18)
                                                                           atTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAACCACTGAACAACGAGAGTTGTTCGTGAGTCTCACAAATTTTCGCAACACGATGATGAATCg                    >  2:1888563/1‑100 (MQ=11)
                                                                            tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCga                   >  1:2333214/1‑100 (MQ=18)
                                                                            tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg                    >  1:1676193/1‑99 (MQ=17)
                                                                             ttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaa                  >  2:1401536/1‑100 (MQ=17)
                                                                              tGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCCACACGATGATGAATCgaaa                 <  2:3910040/100‑1 (MQ=255)
                                                                               gCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaag                <  1:2189736/100‑1 (MQ=18)
                                                                               gCTCGTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAAAGAGAGTTGTTCGTGAGTCTCTAAAATTTTCGCAACACGATGATGAATCgaaag                >  2:1948585/1‑100 (MQ=11)
                                                                                ctctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaaga               >  1:795252/1‑100 (MQ=17)
                                                                                ctctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaaga               <  1:2550618/100‑1 (MQ=18)
                                                                                ctctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaaga               <  2:2418402/100‑1 (MQ=18)
                                                                                  ctTTAAGAATCTGGATCAGGCTAAAAATTGAAACACTTAACAACGAGAGTTGTACGTGCGCCTTGGAAATTTTCGCAACACGATGATGACTcataag                >  2:558976/1‑92 (MQ=2)
                                                                                  ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa             >  1:1094199/1‑100 (MQ=18)
                                                                                  ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa             <  1:64153/100‑1 (MQ=18)
                                                                                   tttAAAAATCTGGACCAAGCTGAAAATTCAAACCCTGAACAACGAGAGTTGGTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAAc            <  2:1443421/100‑1 (MQ=11)
                                                                                   tttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAAc            <  2:3036377/100‑1 (MQ=17)
                                                                                     tAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCACATTTTCGCAACACGATGATGAATCGAAAGAAACAt          >  2:3668436/1‑100 (MQ=255)
                                                                                     tAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGACTCGAAAGAAACa           >  2:1060361/1‑99 (MQ=255)
                                                                                       aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt        >  2:1129067/1‑100 (MQ=17)
                                                                                       aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt        >  2:2118117/1‑100 (MQ=18)
                                                                                         aaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTc      <  1:3036508/100‑1 (MQ=18)
                                                                                         aaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTc      >  2:2469303/1‑100 (MQ=18)
                                                                                          aTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCg     >  1:1657360/1‑100 (MQ=18)
                                                                                           tCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg    >  2:1938465/1‑100 (MQ=18)
                                                                                            cTGGAACAAGCTGAAAATTGAAACACTTAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCg     <  2:1349887/98‑1 (MQ=12)
                                                                                                 |                                                                                               
CGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGT  >  NC_000913/4168443‑4168635

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