Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,810,760 T→C 100% 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.4% 482.0 / ‑6.2 140intergenic (+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 (2/1);  new base C (65/72);  total (67/73)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810663‑2810854
                                                                                                 |                                                                                              
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCa                                                                                              >  2:1065205/1‑100 (MQ=255)
 gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCAt                                                                                             >  1:638541/1‑100 (MQ=255)
 gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCAt                                                                                             >  2:1310447/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                           >  2:2300073/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                           >  1:935108/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                           >  1:1408393/1‑100 (MQ=255)
   tGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAc                                                                                           <  1:1255375/100‑1 (MQ=255)
     gCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATAccc                                                                                         >  2:1475085/1‑100 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCa                                                                                        <  2:2431973/100‑1 (MQ=255)
        aTTATATAACCCACAAGAATCATTTTTCTAAAACAATATATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCaaa                                                                                      <  2:1251668/100‑1 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAATCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                    <  1:733853/100‑1 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAAt                                                                                     >  2:666989/1‑99 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                    <  2:1171792/100‑1 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATg                                                                                    <  1:538957/100‑1 (MQ=255)
            tataGCCCACACGAATCATTTTTACAAAACAAAACATATACTTTATTTGTCACTTTCCGCACTAAATCGCCGGAAATTAATAAGTTAACACCCaatt                                                                                     >  2:1304690/1‑95 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGAt                                                                                 <  1:2207851/100‑1 (MQ=255)
              taACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATa                                                                                <  1:1714391/100‑1 (MQ=255)
               aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAg                                                                               <  1:1606187/100‑1 (MQ=255)
                aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGt                                                                              <  2:3002542/100‑1 (MQ=255)
                aCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGt                                                                              >  1:1732962/1‑100 (MQ=255)
                 cccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGAATAGtt                                                                             >  1:379176/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                           >  2:2421406/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                           <  2:2362656/100‑1 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATAACGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCg                                                                           >  1:610877/1‑100 (MQ=255)
                     caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGtt                                                                         <  2:1557139/100‑1 (MQ=255)
                      aaGAATCATTTTTCTAAACCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGttt                                                                        <  1:1048093/100‑1 (MQ=255)
                      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGttt                                                                        <  1:2408157/100‑1 (MQ=255)
                      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGttt                                                                        <  2:384065/100‑1 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTa                                                                       <  1:2583538/100‑1 (MQ=255)
                        gAATCATTTTTCTAAAACAATACATTTACTTTAGTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTAc                                                                      >  2:2024196/1‑100 (MQ=255)
                         aaTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACg                                                                     >  2:2822056/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGccc                                                                  >  1:638827/1‑100 (MQ=255)
                               ttttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCAtt                                                               <  2:1124688/100‑1 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                              >  2:1399858/1‑100 (MQ=255)
                                tttCTAAAACAACACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTg                                                              <  1:731942/100‑1 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           >  2:1498669/1‑100 (MQ=255)
                                    ttccgatctACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                          <  1:1851419/92‑1 (MQ=255)
                                    tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATCTCCGCTGCAATTAATGCGGTCCTCCCACAATGGATAGTTCGTCTACGACCACTAAACa                                                          <  2:2596063/100‑1 (MQ=255)
                                    tAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGTCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACa                                                          >  2:972946/1‑100 (MQ=255)
                                     aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                         >  1:2792240/1‑100 (MQ=255)
                                     aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                         >  2:2965819/1‑100 (MQ=255)
                                      aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                        >  2:526049/1‑100 (MQ=255)
                                      aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                        <  2:2576543/100‑1 (MQ=255)
                                       aaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAAt                                                       >  1:1300861/1‑100 (MQ=255)
                                        aCAATACATTTACTTTATTTTTCACTTTCGTTACTATATCGTCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  2:1072859/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  1:1987579/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  2:923600/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  1:2707049/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      >  1:1469943/1‑100 (MQ=255)
                                         cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACACAAATGGATAGTTCGTTTACGCCCATTGAACAAATGc                                                     >  1:2352513/1‑100 (MQ=255)
                                           aTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTa                                                   >  1:2038721/1‑100 (MQ=255)
                                            tACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaa                                                  <  1:770741/100‑1 (MQ=255)
                                            tACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaa                                                  <  2:2218002/100‑1 (MQ=255)
                                             aCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaagatcggaa                                                 >  2:1851404/1‑92 (MQ=255)
                                             aCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaa                                                 <  2:1584799/100‑1 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                <  1:233574/100‑1 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                            >  2:1788254/1‑100 (MQ=255)
                                                   aCTTTATTTATCACTGTCGTTAATATATCGGATGAATTCAATGCGACCATACTCAAATAGCTAGTTCGTTTACGCCCATTGACCAAAGGCTAAAATTTcg                                           >  2:946806/1‑100 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          <  1:1498680/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          >  1:1679819/1‑100 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                         <  1:2939520/100‑1 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                         >  1:698931/1‑100 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                         <  2:671285/100‑1 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                         <  1:1806329/100‑1 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                         >  2:939265/1‑100 (MQ=255)
                                                          ttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgc                                    >  1:1658860/1‑100 (MQ=255)
                                                           tGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                   >  2:1994897/1‑100 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  <  1:2779206/100‑1 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  >  2:211278/1‑100 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgcc                                   >  2:2141751/1‑99 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCACATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  >  2:507031/1‑100 (MQ=255)
                                                             tCACTGTCGTTACTATATCGGCTGAAATTAATGAGGCGATACCCAAATGAATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgc                                 >  2:759559/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                >  2:1041171/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                >  1:1777450/1‑100 (MQ=255)
                                                                cTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAGAAATGCTAAAATTTCGCGCCAGCCGCCAc                              <  2:2510378/100‑1 (MQ=255)
                                                                cTGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCAc                              <  2:2510502/100‑1 (MQ=255)
                                                                 tGTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACg                             <  1:2993924/100‑1 (MQ=255)
                                                                  gTCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGa                            >  1:1188487/1‑100 (MQ=255)
                                                                   tCGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGaa                           <  2:651097/100‑1 (MQ=255)
                                                                    cGTTACTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                          >  1:1020795/1‑100 (MQ=255)
                                                                        aCTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                      <  2:1433337/100‑1 (MQ=255)
                                                                        aCTATATCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGAttt                      >  1:1934850/1‑100 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  2:1300852/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  2:987697/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  1:1421885/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  2:760893/100‑1 (MQ=255)
                                                                              tCGGCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                 <  1:1994914/99‑1 (MQ=255)
                                                                               cagCTTAAATTAATTATTCCAGACCCAAATGTATATTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt               <  1:571521/98‑1 (MQ=255)
                                                                                 gCTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCa             <  1:387738/100‑1 (MQ=255)
                                                                                  cTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg            <  1:1337412/100‑1 (MQ=255)
                                                                                  cTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg            >  1:2418904/1‑100 (MQ=255)
                                                                                  cTGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAg            >  1:2268009/1‑100 (MQ=255)
                                                                                   tGAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAgg           <  2:2613278/100‑1 (MQ=255)
                                                                                    gAAATTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGga          <  2:171276/100‑1 (MQ=255)
                                                                                     aaaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag         <  1:2336147/100‑1 (MQ=255)
                                                                                      aaTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga        >  1:1770753/1‑100 (MQ=255)
                                                                                       aTTAATGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       <  1:2827898/100‑1 (MQ=255)
                                                                                        ttAATAAGGCCATACTTAATTGGATAGTTCGTTTACGCCCATTGAACAAATTCTAAAATCTCGCGCTAGCCGCCACGAAGATTTTCCTTATCAGGAGATc      <  2:2849477/100‑1 (MQ=255)
                                                                                            tGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCttn  <  1:2677916/100‑2 (MQ=255)
                                                                                            tGAGGCCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTc  >  2:2066150/1‑100 (MQ=255)
                                                                                                 |                                                                                              
TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810663‑2810854

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: