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,5400TA94.2% 260.0 / 8.0 103intergenic (+93/‑101)gltT/rrlBtRNA‑Glu/23S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base A (52/45);  total (54/49)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 103 total aligned reads displayed.

GCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168444‑4168636
                                                                                                |                                                                                                
gCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACaa                                                                                               >  1:1193006/1‑100 (MQ=255)
  cACTTGCTTGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAATCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                             <  2:2046834/100‑1 (MQ=255)
  cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGc                                                                                             <  1:1202191/100‑1 (MQ=255)
   aCTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCt                                                                                            <  1:1327960/100‑1 (MQ=255)
    cTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTg                                                                                           <  1:2121248/100‑1 (MQ=255)
    cTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTg                                                                                           <  2:2405418/100‑1 (MQ=255)
      tGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaa                                                                                         >  1:1394371/1‑100 (MQ=255)
      tGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaa                                                                                         >  2:2452149/1‑100 (MQ=255)
        cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       <  1:738432/100‑1 (MQ=255)
        cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       <  1:373126/100‑1 (MQ=255)
        cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       >  1:771074/1‑100 (MQ=255)
           gTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTg                                                                                    <  2:1399572/100‑1 (MQ=255)
            tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                   <  2:1834211/100‑1 (MQ=255)
             ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  >  2:222949/1‑100 (MQ=255)
             ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  1:1763664/100‑1 (MQ=255)
             ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  1:114314/100‑1 (MQ=255)
             ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  1:1200403/100‑1 (MQ=255)
               gtgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGaaan                                                                                <  1:1972525/100‑2 (MQ=255)
               gtgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAac                                                                                <  2:1394362/100‑1 (MQ=255)
               gtgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAac                                                                                >  1:1623987/1‑100 (MQ=255)
               gtcagtgcACGTCGACAACATTAATATAAAAAAACTCATCTACGCGTGAGGTTTGAGAAAATTAATAAAAAAAAATCTGGAACAAGCTGAAAATTGa                                                                                   <  2:1109663/94‑1 (MQ=2)
                tgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAaca                                                                               >  2:1497199/1‑100 (MQ=255)
                tgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAaca                                                                               >  2:564223/1‑100 (MQ=255)
                 gagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAacac                                                                              >  2:1209000/1‑100 (MQ=255)
                  agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACt                                                                             <  1:1612654/100‑1 (MQ=255)
                  agtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACt                                                                             >  2:55723/1‑100 (MQ=255)
                   gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTCAAAATCTGGAACAAGCTGAAAATTGAACCACTg                                                                            >  2:1224189/1‑100 (MQ=255)
                   gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            >  1:1984796/1‑100 (MQ=255)
                   gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            >  1:1692659/1‑100 (MQ=255)
                   gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            >  2:692217/1‑100 (MQ=255)
                     gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCCCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           >  2:2115841/1‑99 (MQ=255)
                      aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaac                                                                         <  1:2145333/100‑1 (MQ=255)
                       aaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGaaca                                                                        >  1:1935832/1‑100 (MQ=40)
                        aGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaacaa                                                                       >  1:42787/1‑100 (MQ=255)
                        aGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         <  2:1355450/98‑1 (MQ=255)
                         gTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaaca                                                                        >  2:1551401/1‑98 (MQ=255)
                           cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    >  1:170045/1‑100 (MQ=25)
                           cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    <  2:2356130/100‑1 (MQ=25)
                             ccGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgaga                                                                  <  1:598175/100‑1 (MQ=25)
                              cGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGagag                                                                 >  1:2174268/1‑100 (MQ=25)
                              cGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGagag                                                                 <  1:1855217/100‑1 (MQ=25)
                              cGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGagag                                                                 <  1:449207/100‑1 (MQ=25)
                              cGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGagag                                                                 <  1:1503653/100‑1 (MQ=25)
                               gACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgt                                                                <  1:532181/100‑1 (MQ=25)
                                aCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaacaactctcgt                                                                <  1:2150863/99‑7 (MQ=255)
                                 ccTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGATATATTTTCTCTTTAAAAATCTGGAACAATCTGAAAATTAAAACATTGAATAACGAGAttt                                                               >  2:1074214/1‑99 (MQ=11)
                                 ccTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttg                                                              <  2:895474/100‑1 (MQ=18)
                                 ccTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttg                                                              >  2:1988733/1‑100 (MQ=18)
                                   ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTTCTCTTTAAAAAGCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttg                                                              <  2:2410227/98‑1 (MQ=255)
                                   ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgtt                                                            <  1:67816/100‑1 (MQ=17)
                                   ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAATTGAAACACTGAACAACGAGAGTTGTTc                                                           >  1:1186112/1‑100 (MQ=255)
                                       tatCTCACAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  2:2214718/100‑1 (MQ=12)
                                       tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  1:199098/100‑1 (MQ=17)
                                       tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  2:855596/100‑1 (MQ=18)
                                       tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  2:565361/100‑1 (MQ=17)
                                       tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  1:1961059/100‑1 (MQ=18)
                                         tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  2:1646837/100‑1 (MQ=18)
                                         tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  1:705759/1‑100 (MQ=18)
                                         tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  2:1815601/1‑100 (MQ=17)
                                         tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  1:2083455/1‑100 (MQ=18)
                                          ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      >  2:1692931/1‑99 (MQ=17)
                                           tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtc                                                    <  1:692224/100‑1 (MQ=18)
                                             aaaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctc                                                  <  2:1408049/100‑1 (MQ=17)
                                               aaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctctc                                                >  1:498111/1‑100 (MQ=18)
                                                aCTCATCTTCGGGTGATGTTTGAGATATTTGCACTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTGGTGAGTCTCTCa                                               >  1:2429843/1‑100 (MQ=12)
                                                 cTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                              >  1:2398635/1‑100 (MQ=18)
                                                      cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAtttt                                         >  2:864675/1‑100 (MQ=18)
                                                        tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                       <  2:1935818/100‑1 (MQ=32)
                                                        tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                       >  1:1000580/1‑100 (MQ=18)
                                                        tCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                       >  1:385131/1‑100 (MQ=17)
                                                         cGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGc                                      <  2:2025602/100‑1 (MQ=18)
                                                          gggTGATGTTTGAGATATTTGCTCTTGAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGc                                      >  2:1759146/1‑99 (MQ=255)
                                                            gTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAac                                   <  1:1761262/100‑1 (MQ=18)
                                                             tGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                                  >  1:1643695/1‑100 (MQ=18)
                                                             tGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAaca                                  <  2:669187/100‑1 (MQ=18)
                                                              gATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac                                 >  1:1671405/1‑100 (MQ=18)
                                                                tGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga                               <  2:1777887/100‑1 (MQ=17)
                                                                      agaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGAtgatga                         <  2:46889/100‑1 (MQ=18)
                                                                       gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                        >  2:2398603/1‑100 (MQ=18)
                                                                       gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                        >  2:305595/1‑100 (MQ=18)
                                                                         tatTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                      <  1:2245133/100‑1 (MQ=18)
                                                                          atTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg                     >  1:2251274/1‑100 (MQ=18)
                                                                           tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCga                    >  1:2321491/1‑100 (MQ=18)
                                                                           tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCGCAAATTTTCGCAACACGATGATGAATCga                    >  2:2058474/1‑100 (MQ=255)
                                                                             tGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                  <  2:1775766/100‑1 (MQ=17)
                                                                             tGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                  <  1:174213/100‑1 (MQ=18)
                                                                                 ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGTGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              <  1:88199/100‑1 (MQ=255)
                                                                                 ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  1:169509/1‑100 (MQ=18)
                                                                                 ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  1:1567807/1‑100 (MQ=18)
                                                                                   ttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACa            >  1:1334727/1‑100 (MQ=17)
                                                                                   ttAAAAATCCGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGCAACa            >  2:488792/1‑100 (MQ=255)
                                                                                    tAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACAt           >  1:157051/1‑100 (MQ=18)
                                                                                    tAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACAt           >  2:844164/1‑100 (MQ=17)
                                                                                      aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt         >  1:287413/1‑100 (MQ=18)
                                                                                       aaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCtt        >  1:315401/1‑100 (MQ=18)
                                                                                        aaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTc       >  1:1478995/1‑100 (MQ=18)
                                                                                         aTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCg      >  2:1510164/1‑100 (MQ=18)
                                                                                          tCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     <  2:1522077/100‑1 (MQ=18)
                                                                                           cTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCggg    <  1:781489/100‑1 (MQ=18)
                                                                                             ggAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGtt  >  2:2295572/1‑100 (MQ=18)
                                                                                                |                                                                                                
GCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168444‑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: