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,5400TA99.2% 311.1 / ‑4.2 124intergenic (+93/‑101)gltT/rrlBtRNA‑Glu/23S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (65/58);  total (66/58)
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.21e-01

Only 100 of 124 total aligned reads displayed.

CCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168445‑4168636
                                                                                               |                                                                                                
ccACTTGCTGGTTTGTGAGTGAAAGTCGTCGACCTTAATATCTCAAGTTTCATCTTTGGGTGATGTTTGAGATATTTATTATATAAAAATCTGGAACAAg                                                                                              <  2:664646/100‑1 (MQ=14)
ccACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGAAACAAg                                                                                              >  1:2860579/1‑100 (MQ=255)
 cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                             <  1:998053/100‑1 (MQ=255)
 cACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGc                                                                                             <  1:2825156/100‑1 (MQ=255)
       cTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTAGGGTGATGTTTGAGATATTTGATCTTTAAAAATCTGGAACAAGCTGaaaa                                                                                       >  1:1694890/1‑100 (MQ=255)
         ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAAt                                                                                      >  2:3014752/1‑99 (MQ=255)
         ggTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATATGGAACAAGCTGAAAAtt                                                                                     >  1:525547/1‑100 (MQ=255)
         ggTTTGTCAGTGAAAGTCGCAGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGAATTTTGCTCTTTAAAAATTTGGAACAAGCTGAAAAtt                                                                                     >  2:673791/1‑100 (MQ=14)
           tttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGa                                                                                   <  1:402387/100‑1 (MQ=255)
            ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  2:933075/100‑1 (MQ=255)
            ttGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaa                                                                                  <  2:2243542/100‑1 (MQ=255)
             tgtgAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATTTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaaa                                                                                 <  1:2757391/100‑1 (MQ=255)
             tgtgAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGaaa                                                                                 <  2:2361122/100‑1 (MQ=255)
               tgcgtgaGAGTGGCCTACGTTGGTATTTCAAAACTCATCTTCGTGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAATCTGAAAATTGAAac                                                                                <  2:2873526/96‑1 (MQ=11)
               tgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAaca                                                                               >  2:188742/1‑100 (MQ=255)
               tgagtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAaca                                                                               >  2:2885396/1‑100 (MQ=255)
                  gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            <  2:852515/100‑1 (MQ=255)
                  gtgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTg                                                                            >  1:2943349/1‑100 (MQ=255)
                   tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           <  1:1847822/100‑1 (MQ=255)
                   tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           >  2:1867279/1‑100 (MQ=255)
                   tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           <  1:1116729/100‑1 (MQ=255)
                   tgaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGa                                                                           <  1:1847726/100‑1 (MQ=255)
                    gaAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaa                                                                          <  1:1410584/100‑1 (MQ=255)
                     aaaGTCGCCTACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         <  1:1481781/100‑1 (MQ=12)
                     aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         <  2:503465/100‑1 (MQ=255)
                     aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         >  1:730426/1‑100 (MQ=255)
                     aaaGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGaac                                                                         <  1:1555724/100‑1 (MQ=255)
                         tCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACg                                                                     >  2:2299564/1‑100 (MQ=255)
                          cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    >  1:2549223/1‑100 (MQ=25)
                          cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    >  1:812842/1‑100 (MQ=25)
                          cGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACga                                                                    >  2:1748168/1‑100 (MQ=25)
                            ccGACCTTAATATCTCAAAACTTATCTTCGGCTGATGCTTGAGATATGTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACgaga                                                                  <  2:2048128/100‑1 (MQ=11)
                             cGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATGTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACg                                                                     <  2:2528256/96‑1 (MQ=12)
                                ccTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttg                                                              >  1:1831832/1‑100 (MQ=18)
                                ccTTAATATCTCAAAACTCATATTCGGATGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGagag                                                                 >  2:1486656/1‑97 (MQ=255)
                                  ttAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAgttgtt                                                            >  1:2584925/1‑100 (MQ=18)
                                    aaTATCTCAAAACTCATCTTCGGGTTATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCg                                                          <  2:2653995/100‑1 (MQ=12)
                                    aaTATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCg                                                          <  2:8275/100‑1 (MQ=18)
                                      tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTTCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  2:133753/100‑1 (MQ=255)
                                      tatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTg                                                        <  1:2010379/100‑1 (MQ=17)
                                       atCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAACTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGa                                                       >  2:2807681/1‑100 (MQ=12)
                                       atCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGa                                                       >  2:2664881/1‑100 (MQ=17)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  1:2973266/100‑1 (MQ=18)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  1:2236407/100‑1 (MQ=18)
                                        tctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAg                                                      <  1:2236054/100‑1 (MQ=18)
                                         ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                     >  1:1780792/1‑100 (MQ=17)
                                         ctcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGt                                                     >  1:2347177/1‑100 (MQ=18)
                                          tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAGTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAgac                                                    >  2:919672/1‑98 (MQ=255)
                                          tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtc                                                    <  2:160594/100‑1 (MQ=18)
                                          tcAAAACTCATCTTCGGGTGATGGATGAGATATTTGCTCTTAAAAAATCTGTAACAAGCTGAAAATTGAAAAACTGAACAATGAGAGTTGTTCGTGAg                                                      >  2:1070340/1‑98 (MQ=11)
                                           cAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtct                                                   >  2:2146510/1‑100 (MQ=18)
                                             aaaaTCATCTTCGGGTGATGTTTGTGATATTTGCTTTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGtctct                                                 <  1:218600/100‑1 (MQ=11)
                                               aCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCa                                               >  2:2640420/1‑100 (MQ=17)
                                                tctaTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                              <  1:1918017/97‑1 (MQ=16)
                                                cTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                              <  2:1942814/100‑1 (MQ=17)
                                                cTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                              <  2:685863/100‑1 (MQ=18)
                                                 tCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAACATTGAAACCCTGAACAACGAGAGTTGTTCGTGAGTCTCTCaa                                              <  1:1505732/99‑1 (MQ=255)
                                                 tCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaaa                                             <  2:1949306/100‑1 (MQ=18)
                                                  cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAt                                            <  2:1225307/100‑1 (MQ=17)
                                                   aTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCaaaga                                           >  2:1917994/1‑98 (MQ=16)
                                                   aTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAAtt                                           <  1:1954854/100‑1 (MQ=18)
                                                      ttCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGCAACTTGAAACACTGAACAACCAGAGTTGTTCGTGAGTCTCTCAAATTTTc                                        <  2:1805645/100‑1 (MQ=11)
                                                      ttCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTc                                        <  1:1953522/100‑1 (MQ=18)
                                                       tccgatcttGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCg                                       <  1:2425641/92‑1 (MQ=11)
                                                        cGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAATCTGAAAATTGAAACACTGTACAACGAGAGTTGTTCTTGAGTCTCTCAAATTTgcgc                                      <  2:1400258/100‑1 (MQ=11)
                                                        cGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGc                                      <  1:1224615/100‑1 (MQ=17)
                                                          ggTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCaa                                    <  2:426681/100‑1 (MQ=18)
                                                               tGTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCgagatcgga                               >  2:2425612/1‑92 (MQ=11)
                                                                gTTTGAGATATTTGCTCTTTAAAAATCTGGAACGAGCTGAAAATTGAAACACTGAACAACGAGATTTGTTCGTGAGTCTCTCAAATTTTTGCAACACgat                              >  2:1999575/1‑100 (MQ=11)
                                                                gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  2:1048791/1‑100 (MQ=18)
                                                                gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  1:1382604/1‑100 (MQ=17)
                                                                gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  1:468074/1‑100 (MQ=17)
                                                                gTTTGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat                              >  2:2033046/1‑100 (MQ=17)
                                                                  ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTAGTGAGTCTCTCAAATTTTCGCAACACgatga                            >  2:2725357/1‑100 (MQ=12)
                                                                  ttGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAAAAACGAGAATTGTTAGTGAGTATATCAAATTTTCGCAACACgatga                            >  1:954563/1‑100 (MQ=11)
                                                                   tGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat                           >  2:2829544/1‑100 (MQ=18)
                                                                   tGAGATATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat                           >  1:1225489/1‑100 (MQ=17)
                                                                     agaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTTTCTCAAATTTTAGTAACACgataatta                         >  1:596956/1‑98 (MQ=11)
                                                                     agaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgatn                          <  2:2620152/99‑2 (MQ=18)
                                                                      gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                        >  2:122677/1‑100 (MQ=17)
                                                                      gaTATTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa                        >  2:1184202/1‑100 (MQ=18)
                                                                         atTTGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg                     >  2:1416221/1‑100 (MQ=17)
                                                                          tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCga                    >  2:171343/1‑100 (MQ=18)
                                                                          tttGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCga                    >  2:1155431/1‑100 (MQ=17)
                                                                            tGCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                  <  1:869666/100‑1 (MQ=17)
                                                                             gCTCTTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaag                 <  2:1953771/100‑1 (MQ=18)
                                                                              ctctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaaga                >  2:2291579/1‑100 (MQ=18)
                                                                               tctTTTAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGCCTCTCAAATTGTCGAAACACGATGATGAATCgaaaga                >  2:2518327/1‑99 (MQ=11)
                                                                               tctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa               >  2:299989/1‑100 (MQ=17)
                                                                                ntttAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa               >  2:1561575/2‑99 (MQ=17)
                                                                                ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  2:2519723/1‑100 (MQ=17)
                                                                                ctTTAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa              >  2:1883132/1‑100 (MQ=18)
                                                                                   tAAAAATCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACAt           >  2:1206869/1‑100 (MQ=18)
                                                                                    aaaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACaac          >  2:1842047/1‑98 (MQ=16)
                                                                                     aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt         >  2:1159143/1‑100 (MQ=18)
                                                                                        aTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCg      >  2:557601/1‑100 (MQ=17)
                                                                                         tCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg     <  2:2347150/100‑1 (MQ=18)
                                                                                          cTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCggg    <  1:999751/100‑1 (MQ=17)
                                                                                          cTGGAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCggg    <  2:1279255/100‑1 (MQ=17)
                                                                                            ggAACAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGtt  >  1:2359038/1‑100 (MQ=18)
                                                                                               |                                                                                                
CCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAGAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTT  >  NC_000913/4168445‑4168636

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