Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,756 G→A intergenic (+77/‑14) ygaH → / → mprA putative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,7560GA96.3% 366.1 / 1.1 108intergenic (+77/‑14)ygaH/mprAputative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (57/47);  total (59/49)
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.99e-01

Only 100 of 108 total aligned reads displayed.

TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGA  >  NC_000913/2810663‑2810848
                                                                                             |                                                                                            
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCa                                                                                        <  1:34867/100‑1 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAc                                                                                     >  1:1094462/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCGAAAACAATACATTTACTTTATTTGTCACTGTCGTTAGTATATCGGCTGAAATTAATAAAGTCATAc                                                                                     >  2:1403073/1‑100 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAcc                                                                                    >  2:261357/1‑100 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATAccc                                                                                   <  1:1451962/100‑1 (MQ=255)
      cttccgatctaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCa                                                                                  <  1:58085/90‑1 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCa                                                                                  >  1:1969889/1‑100 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCa                                                                                  >  1:1300891/1‑100 (MQ=255)
        aTTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCaaa                                                                                <  2:26997/100‑1 (MQ=255)
         ttATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAAt                                                                               <  1:498813/100‑1 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATg                                                                              <  2:819628/100‑1 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATg                                                                              <  1:456805/100‑1 (MQ=255)
           ctataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATgg                                                                             <  2:1131083/99‑1 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATgg                                                                             >  2:431088/1‑100 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATgg                                                                             >  1:624721/1‑100 (MQ=255)
            tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGa                                                                            >  1:1224411/1‑100 (MQ=255)
            tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGa                                                                            <  2:2295668/100‑1 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGag                                                                           >  1:1131085/1‑99 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATa                                                                          >  2:1782433/1‑100 (MQ=255)
                aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCaagatcggaag                                                                        >  2:58083/1‑91 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCg                                                                     >  2:1857251/1‑100 (MQ=255)
                    acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGt                                                                    <  1:431090/100‑1 (MQ=255)
                    acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGt                                                                    <  1:1674124/100‑1 (MQ=255)
                      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGttt                                                                  >  2:592365/1‑100 (MQ=255)
                      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGttt                                                                  >  2:591664/1‑100 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTa                                                                 <  2:547366/100‑1 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTa                                                                 >  2:5580/1‑100 (MQ=255)
                         aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACg                                                               >  1:1783122/1‑100 (MQ=255)
                         aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACg                                                               >  1:1591023/1‑100 (MQ=255)
                         aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCAATGTAGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACg                                                               >  1:1365113/1‑100 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                              >  2:1144904/1‑100 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGc                                                              <  1:1530016/100‑1 (MQ=255)
                           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGcc                                                             <  1:2007090/100‑1 (MQ=255)
                             gatcttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCa                                                           <  2:2148373/96‑1 (MQ=255)
                             aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCa                                                           >  1:847523/1‑100 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                        <  1:186222/100‑1 (MQ=255)
                                nttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                        >  2:406034/2‑100 (MQ=255)
                                 ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTGTACGCCCATTGa                                                       >  2:1148719/1‑100 (MQ=255)
                                 ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTGTACGCCCATTGa                                                       >  2:1148575/1‑100 (MQ=255)
                                  tCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCaagatc                                                      >  1:2148380/1‑95 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                     >  2:275279/1‑100 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTGTACGCCCATTGAAc                                                     <  1:1148721/100‑1 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTGTACGCCCATTGAAc                                                     <  1:1148577/100‑1 (MQ=255)
                                    tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                    >  2:1289557/1‑100 (MQ=255)
                                      aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                  <  2:1591018/100‑1 (MQ=255)
                                      aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCAAAATTGGTAGTGCGTTTAAGCCAATAGAACaaa                                                  >  2:679748/1‑100 (MQ=255)
                                       aaCAATACATTGACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                 <  1:615968/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                <  1:971386/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACACAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                 >  2:2322571/1‑99 (MQ=255)
                                         cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                               >  2:468402/1‑100 (MQ=255)
                                         cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                               <  1:1782438/100‑1 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACagatcggaaga                                          >  1:1971549/1‑90 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                          >  2:1759804/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                          >  2:257361/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                          >  1:2036568/1‑100 (MQ=255)
                                               aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                         >  1:2003055/1‑100 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                      <  1:2279928/100‑1 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                      <  1:2078744/100‑1 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                      >  1:2019724/1‑100 (MQ=255)
                                                   aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                     <  1:5580/100‑1 (MQ=255)
                                                   aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                     <  1:1149709/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                    >  1:347464/1‑100 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                    <  1:1144906/100‑1 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                   <  1:209823/100‑1 (MQ=255)
                                                       tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                 >  1:1138097/1‑100 (MQ=255)
                                                        aTTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCa                                >  1:1130657/1‑100 (MQ=255)
                                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                               <  2:1340992/100‑1 (MQ=255)
                                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                               <  2:2051425/100‑1 (MQ=255)
                                                          ttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgc                              >  2:1130117/1‑100 (MQ=255)
                                                          ttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgc                              <  1:1025782/100‑1 (MQ=255)
                                                          ttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATCGAACAAATGCTAAAATTTcg                                     >  2:2237558/1‑93 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                             >  1:1886447/1‑100 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                             <  2:1985299/100‑1 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                            <  2:2031215/100‑1 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                          <  1:942090/100‑1 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                          >  1:1914270/1‑100 (MQ=255)
                                                               aCTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCa                         >  2:1862412/1‑100 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                      >  2:730658/1‑100 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                      >  2:1461575/1‑100 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCAGCCACGa                      >  2:1912043/1‑100 (MQ=255)
                                                                    cGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                    >  1:2198584/1‑100 (MQ=255)
                                                                      ttACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAt                  <  2:1997782/100‑1 (MQ=255)
                                                                       tACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                 >  2:557021/1‑100 (MQ=255)
                                                                        aCTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                <  1:1461579/100‑1 (MQ=255)
                                                                         cTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt               >  2:87010/1‑100 (MQ=255)
                                                                         cTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt               >  1:1738827/1‑100 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcn             <  1:406036/100‑2 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  2:1528751/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc             <  2:847521/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt            >  1:1642086/1‑100 (MQ=255)
                                                                            tatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt            <  2:1886440/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt            <  1:1130119/100‑1 (MQ=255)
                                                                             atCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt           <  1:150328/100‑1 (MQ=255)
                                                                             atCGGCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTACt            >  2:2213690/1‑99 (MQ=255)
                                                                                ggCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc        >  1:549346/1‑100 (MQ=255)
                                                                                 gCTGAAATTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa       <  1:1872444/100‑1 (MQ=255)
                                                                                     aaaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag   <  1:730659/100‑1 (MQ=255)
                                                                                      aaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga  >  2:710084/1‑100 (MQ=255)
                                                                                      aaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga  >  2:1214973/1‑100 (MQ=255)
                                                                                      aaTTAATAAAGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga  >  2:1218041/1‑100 (MQ=255)
                                                                                             |                                                                                            
TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGA  >  NC_000913/2810663‑2810848

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