Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,760 T→C intergenic (+81/‑10) 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,7600TC99.3% 498.0 / ‑5.0 137intergenic (+81/‑10)ygaH/mprAputative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (64/72);  total (64/73)
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.76e-01

Only 100 of 137 total aligned reads displayed.

GATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTG  >  NC_000913/2810664‑2810856
                                                                                                |                                                                                                
gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCAt                                                                                               <  2:1689399/100‑1 (MQ=255)
gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCAt                                                                                               >  2:359261/1‑100 (MQ=255)
  tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                             >  1:187916/1‑100 (MQ=255)
   ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAcc                                                                                            >  1:912210/1‑100 (MQ=255)
   ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAcc                                                                                            <  2:3462350/100‑1 (MQ=255)
         tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                      <  1:2977023/100‑1 (MQ=255)
         tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                      <  1:1554477/100‑1 (MQ=255)
         tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                      <  2:3258801/100‑1 (MQ=255)
          atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATgg                                                                                     >  2:2843486/1‑100 (MQ=255)
           tataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGa                                                                                    >  2:3178625/1‑100 (MQ=255)
            ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGAt                                                                                   >  2:2199812/1‑100 (MQ=255)
            ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGAt                                                                                   <  1:2877663/100‑1 (MQ=255)
             taACCCACATGAATCGTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGAt                                                                                   <  2:3431459/99‑1 (MQ=255)
             taACCCACAAGGATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAGATGGATa                                                                                  <  1:2688324/100‑1 (MQ=255)
             taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATa                                                                                  >  1:238272/1‑100 (MQ=255)
              aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATCTCGGCTGAAATTAATGAGGCCATACCCAAATGGATAg                                                                                 <  1:1893228/100‑1 (MQ=255)
              aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAg                                                                                 <  2:2137595/100‑1 (MQ=255)
               aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGt                                                                                >  1:3184851/1‑100 (MQ=255)
                cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGtt                                                                               >  2:1797581/1‑100 (MQ=255)
                cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGtt                                                                               >  2:1847274/1‑100 (MQ=255)
                 ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGttn                                                                              <  1:258501/100‑2 (MQ=255)
                  cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                             >  2:436219/1‑100 (MQ=255)
                   acaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGt                                                                            >  1:2927007/1‑100 (MQ=255)
                    caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCAGTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGtt                                                                           <  1:2257464/100‑1 (MQ=255)
                    caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGtt                                                                           <  2:1354580/100‑1 (MQ=255)
                    caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGtt                                                                           <  1:1847264/100‑1 (MQ=255)
                     aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGttt                                                                          <  1:888258/100‑1 (MQ=255)
                     aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGttt                                                                          <  1:3097370/100‑1 (MQ=255)
                      ataaTCATTTTTCTAAAACAATACATTTACTGTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGt                                                                            <  1:2939587/95‑1 (MQ=255)
                      aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTa                                                                         <  2:3231960/100‑1 (MQ=255)
                         aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGc                                                                      <  1:66821/100‑1 (MQ=255)
                          tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGcc                                                                     <  1:447138/100‑1 (MQ=255)
                            aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCa                                                                   >  2:898461/1‑100 (MQ=255)
                            aTTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCa                                                                   >  1:2101167/1‑100 (MQ=255)
                             tttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCAt                                                                  >  1:2722709/1‑100 (MQ=255)
                              ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                                 >  1:661626/1‑100 (MQ=255)
                              ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                                 <  2:2017864/100‑1 (MQ=255)
                               tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                                <  2:1917048/100‑1 (MQ=255)
                               tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                                <  2:2452204/100‑1 (MQ=255)
                               tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                                <  2:511418/100‑1 (MQ=255)
                                  cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             >  2:3295280/1‑100 (MQ=255)
                                  cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             >  2:3434969/1‑100 (MQ=255)
                                  cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             <  1:2082178/100‑1 (MQ=255)
                                  cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                             >  1:2609420/1‑100 (MQ=255)
                                   tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                            >  1:3252908/1‑100 (MQ=255)
                                    aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                           <  1:1160627/100‑1 (MQ=255)
                                     aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                          <  1:2669347/100‑1 (MQ=255)
                                     aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                          >  1:2976327/1‑100 (MQ=255)
                                      aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                         <  2:1009995/100‑1 (MQ=255)
                                       aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                        >  2:771406/1‑100 (MQ=255)
                                            aCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaa                                                   <  1:3350893/100‑1 (MQ=255)
                                             cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                  >  1:3345262/1‑100 (MQ=255)
                                              aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                                 <  2:2990543/100‑1 (MQ=255)
                                              aTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAAt                                                 <  2:132817/100‑1 (MQ=255)
                                                 tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                              >  2:3151230/1‑100 (MQ=255)
                                                  aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                             >  1:859314/1‑100 (MQ=255)
                                                  aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                             <  2:3497536/100‑1 (MQ=255)
                                                   cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                            <  2:575435/100‑1 (MQ=255)
                                                   cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                            <  2:385891/100‑1 (MQ=255)
                                                   cTTTATTTGTCACTGTCGTTACTATATAGGCTGAAATTAATGAGGCCATACACAAATCGATAGTTCGTTTAAGCCCATTGAACAAAGGCTAAAATTTcgc                                            >  1:855520/1‑100 (MQ=255)
                                                    tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                           <  2:199434/100‑1 (MQ=255)
                                                      tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                         >  1:37498/1‑100 (MQ=255)
                                                          tGTCACTGTCGTTACTATATCGGCTGAAATTAACGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                     >  1:2119921/1‑100 (MQ=255)
                                                             cACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                  <  2:2096754/100‑1 (MQ=255)
                                                             cACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                  >  2:2778205/1‑100 (MQ=255)
                                                               cTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                >  1:2349002/1‑100 (MQ=255)
                                                               cTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                                >  1:1637668/1‑100 (MQ=255)
                                                                 gTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                              >  2:152207/1‑100 (MQ=255)
                                                                   cGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                            >  2:1038300/1‑100 (MQ=255)
                                                                    gTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                           >  2:169291/1‑100 (MQ=255)
                                                                      tACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAtt                         >  1:122571/1‑100 (MQ=255)
                                                                         tataTCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                      <  2:661629/100‑1 (MQ=255)
                                                                          cgatctgCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     <  1:3470158/94‑1 (MQ=255)
                                                                          atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     >  2:1489908/1‑100 (MQ=255)
                                                                          atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     <  2:1066102/100‑1 (MQ=255)
                                                                          atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     >  2:1052620/1‑100 (MQ=255)
                                                                          atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                     <  1:792140/100‑1 (MQ=255)
                                                                            atCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                   <  2:912214/100‑1 (MQ=255)
                                                                               ggCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc                >  1:34203/1‑100 (MQ=255)
                                                                               ggCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc                >  2:3374081/1‑100 (MQ=255)
                                                                               ggCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATACTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc                <  1:3178605/100‑1 (MQ=255)
                                                                                tcTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa               <  2:1715802/99‑1 (MQ=255)
                                                                                 cTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCaa              >  1:1715796/1‑99 (MQ=255)
                                                                                 cTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg              <  1:2258303/100‑1 (MQ=255)
                                                                                  tGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg             >  1:1974355/1‑100 (MQ=255)
                                                                                   gAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga            <  1:905727/100‑1 (MQ=255)
                                                                                   gAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga            <  1:436217/100‑1 (MQ=255)
                                                                                   gAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga            <  1:3307680/100‑1 (MQ=255)
                                                                                     aaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga          <  1:1359086/100‑1 (MQ=255)
                                                                                     aaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga          >  1:1191200/1‑100 (MQ=255)
                                                                                     aaTTAATGAGGCCATACCCAAATGGATAGATCGTTTACGCCCATTGAACAAAGGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga          >  2:3377839/1‑100 (MQ=255)
                                                                                      aTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt         <  2:1019063/100‑1 (MQ=255)
                                                                                        tAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATcc       <  1:400239/100‑1 (MQ=255)
                                                                                         aaTGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt      <  2:1468315/100‑1 (MQ=255)
                                                                                         aaTGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt      <  2:1225741/100‑1 (MQ=255)
                                                                                         aaTGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt      <  1:2427921/100‑1 (MQ=255)
                                                                                           tGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc    <  2:1184613/100‑1 (MQ=255)
                                                                                           tGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc    <  1:3151210/100‑1 (MQ=255)
                                                                                            gAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCt   <  1:2085497/100‑1 (MQ=255)
                                                                                             aGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTg  <  2:1637674/100‑1 (MQ=255)
                                                                                                |                                                                                                
GATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTG  >  NC_000913/2810664‑2810856

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