Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,756 G→A 100% 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,7560GA99.1% 364.4 / ‑4.5 110intergenic (+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 (0/1);  new base A (59/50);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 110 total aligned reads displayed.

AGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCT  >  NC_000913/2810661‑2810852
                                                                                               |                                                                                                
aGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGt                                                                                              <  1:1072945/100‑1 (MQ=255)
aGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGt                                                                                              <  1:298022/100‑1 (MQ=255)
 gTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTAATTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGt                                                                                              <  2:390586/99‑1 (MQ=255)
     tGGCGATTATATAACCCACAGGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAc                                                                                         >  2:1564318/1‑100 (MQ=255)
     tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATAc                                                                                         >  2:24506/1‑100 (MQ=255)
        cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCa                                                                                      >  1:57291/1‑100 (MQ=255)
             atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATgg                                                                                 <  1:1693973/100‑1 (MQ=255)
                taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATa                                                                              >  2:2063963/1‑100 (MQ=255)
                taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATa                                                                              >  1:1322728/1‑100 (MQ=255)
                 aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAg                                                                             <  1:1398371/100‑1 (MQ=255)
                  aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGt                                                                            >  1:1178294/1‑100 (MQ=255)
                   cccACAAGAATCATTTTTTTTTTTCTTTTCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAg                                                                             >  2:380825/1‑98 (MQ=255)
                   cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGtt                                                                           >  1:435950/1‑100 (MQ=255)
                    ccccaAGAATCATTTTTCTAAAAAAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGtt                                                                           <  1:383508/96‑1 (MQ=255)
                    ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTc                                                                          <  1:521492/100‑1 (MQ=255)
                    ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTc                                                                          <  1:220180/100‑1 (MQ=255)
                     cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCg                                                                         >  2:2074669/1‑100 (MQ=255)
                      acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGt                                                                        <  2:1355416/100‑1 (MQ=255)
                       caAGAATCATTTTTTTTTATTATTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGtt                                                                       >  2:322259/1‑100 (MQ=255)
                        aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGttt                                                                      >  1:1846649/1‑100 (MQ=255)
                         aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGttt                                                                      <  2:458176/99‑1 (MQ=255)
                         aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTa                                                                     <  1:706469/100‑1 (MQ=255)
                         aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTa                                                                     <  1:523391/100‑1 (MQ=255)
                         aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTAGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTa                                                                     >  1:1526311/1‑100 (MQ=255)
                           aaTCATTTTTATAAAACAATAAATTTACTTTATTTGTTACTGTCGTTACTATATCGGCTGAAATTAATAAGTTCATACGCAAATGGATAGTTAGttttc                                                                    >  2:617441/1‑97 (MQ=255)
                                tttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCAt                                                              >  1:516563/1‑100 (MQ=255)
                                  tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                            <  2:699040/100‑1 (MQ=255)
                                  tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                            >  1:883586/1‑100 (MQ=255)
                                  tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                            <  1:1823923/100‑1 (MQ=255)
                                   ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                           <  2:496490/100‑1 (MQ=255)
                                   ttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACAATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGa                                                           <  1:637059/100‑1 (MQ=255)
                                    tCTAAAACAATACATTTACTTTATTTGTCACTGTCGTCACTATATCGGCTGAAATTAATAAGGTCATACACAAATGAATAGTTCGGTTACTCCCATTg                                                            >  2:144993/1‑98 (MQ=255)
                                     cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                         <  1:473877/100‑1 (MQ=255)
                                     cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                         <  1:792103/100‑1 (MQ=255)
                                     cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                         >  1:2443489/1‑100 (MQ=255)
                                      tAAACCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                         <  1:975299/99‑1 (MQ=255)
                                      tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                        >  1:1901486/1‑100 (MQ=255)
                                        aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                      >  1:2434827/1‑100 (MQ=255)
                                          aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                    >  2:636125/1‑100 (MQ=255)
                                           cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                                   <  2:2379768/100‑1 (MQ=255)
                                           cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                                   <  1:1099715/100‑1 (MQ=255)
                                                cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                              >  1:1629770/1‑100 (MQ=255)
                                                 aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                             >  1:1327581/1‑100 (MQ=255)
                                                 aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                             >  2:1895122/1‑100 (MQ=255)
                                                  tttACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAtt                                            >  2:178733/1‑100 (MQ=255)
                                                     attttATTTGTCACTGTCGTTACTAGTTCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  2:2093880/98‑1 (MQ=255)
                                                     aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                         <  2:2305197/100‑1 (MQ=255)
                                                      cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  1:2109790/1‑100 (MQ=255)
                                                      cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  1:1032298/1‑100 (MQ=255)
                                                      cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  2:1654371/1‑100 (MQ=255)
                                                      cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                        >  1:2328042/1‑100 (MQ=255)
                                                       tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                       >  1:2328476/1‑100 (MQ=255)
                                                       tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                       >  1:424557/1‑100 (MQ=255)
                                                        ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCgcgc                                      <  2:424553/100‑1 (MQ=255)
                                                        ttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCgcgc                                      >  2:2290666/1‑100 (MQ=255)
                                                         tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                     <  1:1714902/100‑1 (MQ=255)
                                                           tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                   <  1:161954/100‑1 (MQ=255)
                                                           tttGTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                   <  1:160679/100‑1 (MQ=255)
                                                              gTCACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                <  1:1512315/100‑1 (MQ=255)
                                                               ncACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgc                               >  2:1693602/2‑100 (MQ=255)
                                                                cACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                              <  2:585320/100‑1 (MQ=255)
                                                                cACTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                              >  1:1307108/1‑100 (MQ=255)
                                                                 aCTGTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCa                             >  1:346877/1‑100 (MQ=255)
                                                                  cTGTCGTTACAATATGGGCTGAAATTAATAAGGTCCTACCCCAATGGATAGTTCGCTTACGCCCATTGAACAACTGGTAAAATTTCGCCCCAGCCGCTAc                            >  2:502496/1‑100 (MQ=255)
                                                                    gTCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                          >  2:515004/1‑100 (MQ=255)
                                                                     tCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                         >  2:1831514/1‑100 (MQ=255)
                                                                     tCGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                         <  1:6526/100‑1 (MQ=255)
                                                                      cGTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                        >  2:1389383/1‑100 (MQ=255)
                                                                       gTTACTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                       <  2:222317/100‑1 (MQ=255)
                                                                           cTATATCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtttt                   <  1:771033/100‑1 (MQ=255)
                                                                            tataTCGGCTGAAATCAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                  >  1:2312481/1‑100 (MQ=255)
                                                                             atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:361773/100‑1 (MQ=255)
                                                                             atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  1:199014/100‑1 (MQ=255)
                                                                             atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  2:880966/100‑1 (MQ=255)
                                                                             atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  1:2063977/100‑1 (MQ=255)
                                                                             atatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                 <  1:2068247/100‑1 (MQ=255)
                                                                              tatCGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                >  2:374237/1‑100 (MQ=255)
                                                                                tccgatctaaTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTa              <  1:309746/92‑1 (MQ=255)
                                                                                 cGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt             >  1:1263476/1‑100 (MQ=255)
                                                                                 cGGCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt             >  1:1181624/1‑100 (MQ=255)
                                                                                  ggCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc            >  1:1329122/1‑100 (MQ=255)
                                                                                  ggCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc            >  1:1425430/1‑100 (MQ=255)
                                                                                  ggCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc            >  1:744243/1‑100 (MQ=255)
                                                                                   gCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa           <  1:2027661/100‑1 (MQ=255)
                                                                                   gCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa           <  1:1831528/100‑1 (MQ=255)
                                                                                   gCTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa           >  2:433542/1‑100 (MQ=255)
                                                                                    cTGAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATccg          <  1:1736071/100‑3 (MQ=255)
                                                                                      gAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga        <  2:1298768/100‑1 (MQ=255)
                                                                                      gAAATTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga        >  1:29402/1‑100 (MQ=255)
                                                                                       aaaTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag       <  2:822129/100‑1 (MQ=255)
                                                                                        aaTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTaagatcgga      >  2:309742/1‑92 (MQ=255)
                                                                                        aaTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga      >  2:23662/1‑100 (MQ=255)
                                                                                         aTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt     <  2:2261512/100‑1 (MQ=255)
                                                                                         aTTAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt     <  2:80355/100‑1 (MQ=255)
                                                                                          ttAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATc    >  1:537987/1‑100 (MQ=255)
                                                                                          ttAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATc    <  2:423499/100‑1 (MQ=255)
                                                                                          ttAATAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATc    >  2:2459951/1‑100 (MQ=255)
                                                                                            aaTAAGGTCATACGCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCTCCAGCCGCCACGAAGATTTTCCTTATCAGGCGATCCt  >  2:1581037/1‑100 (MQ=255)
                                                                                            aaTAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGCTCCt  >  2:2399581/1‑100 (MQ=255)
                                                                                            aaTAAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt  >  2:2198083/1‑100 (MQ=255)
                                                                                               |                                                                                                
AGTGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCT  >  NC_000913/2810661‑2810852

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: