Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,505,561 T→C 54.0% M174V (ATG→GTG)  ydcO ← BenE family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,505,5610TC54.0% 34.6 / 144.7 120M174V (ATG→GTG) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (21/33);  new base C (34/30);  total (55/65)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 122 total aligned reads displayed.

AAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCA  >  NC_000913/1505462‑1505658
                                                                                                   |                                                                                                 
aaaaaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAatcat                                                                                                   >  1:2407357/1‑100 (MQ=255)
  aaaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATcg                                                                                                 <  1:1451553/100‑1 (MQ=255)
  aaaTCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACcg                                                                                                 <  2:2890599/100‑1 (MQ=255)
    gatctgggTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACcgcg                                                                                               <  1:2751613/95‑1 (MQ=255)
     tCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                                                                              <  2:3844/100‑1 (MQ=255)
     tCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                                                                              <  2:3944463/100‑1 (MQ=255)
     tCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCgg                                                                                              <  2:1381529/100‑1 (MQ=255)
     tCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCgg                                                                                              <  2:860062/100‑1 (MQ=255)
      cAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGc                                                                                             <  1:797557/100‑1 (MQ=255)
        ggggTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCaa                                                                                           >  1:1673906/1‑100 (MQ=255)
        ggggTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCaa                                                                                           >  1:1701153/1‑100 (MQ=255)
        ggggTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCaa                                                                                           >  1:3727757/1‑100 (MQ=255)
         gggTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACcgcgagatc                                                                                          >  2:2751597/1‑95 (MQ=255)
              atttaAGTGGGGAGAACGGGTTTAAAGACACCATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCCCCGCGGCAATTACCg                                                                                     <  2:3623249/97‑1 (MQ=255)
              atatAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCg                                                                                     <  1:1489103/100‑1 (MQ=255)
               tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGGCGATCCCGATCCCAATAATCACCGCGGCAATTACCGc                                                                                    >  2:2786991/1‑100 (MQ=255)
               tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGc                                                                                    >  1:3725946/1‑100 (MQ=255)
               tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGc                                                                                    >  1:3726081/1‑100 (MQ=255)
               tataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGc                                                                                    <  1:571585/100‑1 (MQ=255)
                ataAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCa                                                                                   >  1:3226306/1‑100 (MQ=255)
                       gggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGcgcg                                                                            <  2:907195/100‑1 (MQ=255)
                       gggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGcgcg                                                                            >  1:3328968/1‑100 (MQ=255)
                        ggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGcgcg                                                                            >  2:1850143/1‑99 (MQ=255)
                        ggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcgg                                                                           >  2:748530/1‑100 (MQ=255)
                        ggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCgcgg                                                                           <  1:2595796/100‑1 (MQ=255)
                        ggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCgcgg                                                                           >  1:766407/1‑100 (MQ=255)
                        ggAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCgcgg                                                                           >  1:424414/1‑100 (MQ=255)
                          agaACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCgcggcg                                                                         <  1:248413/100‑1 (MQ=255)
                          agaACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCgcggcg                                                                         <  2:4115995/100‑1 (MQ=255)
                             aCGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCaa                                                                      >  1:1822027/1‑100 (MQ=255)
                              cGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAAc                                                                     >  2:3704651/1‑100 (MQ=255)
                              cGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAAc                                                                     >  1:3266660/1‑100 (MQ=255)
                              cGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAAc                                                                     >  1:135269/1‑100 (MQ=255)
                               gggTTTAGTGACAACAACAGTGGGGACAACGTCAGCTTGCGCGACGAAGATAACGAGCCCAATTACCACCGCGGAAAGTACCGCATAGCgcggcg                                                                         <  2:537842/95‑1 (MQ=255)
                                ggTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACgg                                                                   >  1:3482731/1‑100 (MQ=255)
                                ggTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACAGCGGAAATTACCGCATAGCGCGGCGCAACgg                                                                   >  1:2460403/1‑100 (MQ=255)
                                  tttAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGAAGATCACGATCCCAATAATCACCGCGGCTATTACCGCATAGCGCGGCGCAACGGc                                                                  >  2:1015947/1‑99 (MQ=255)
                                   ttAAAGACAACGTCAGTTGTGACAACGTCACCTTGCGCGATGACGATAACGATCCCAATACTCCCCGCGGCAATTACCGCATAACGCCGCGACACGGCCt                                                                >  2:2162719/1‑100 (MQ=255)
                                     aaaGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTg                                                              >  1:610208/1‑100 (MQ=255)
                                      aaGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTgg                                                             <  1:2044057/100‑1 (MQ=255)
                                      aaGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTgg                                                             <  1:748504/100‑1 (MQ=255)
                                       aGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGt                                                            <  2:1919047/100‑1 (MQ=255)
                                        gACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGtt                                                           >  1:3751759/1‑100 (MQ=255)
                                         acaacaTCAGTTGTGACAACGTCACCTTGCGCTATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTg                                                          <  2:287319/100‑1 (MQ=255)
                                         acaacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTg                                                          <  2:3725922/100‑1 (MQ=255)
                                         acaacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGtt                                                           <  1:3106097/99‑1 (MQ=255)
                                          caacaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgc                                                         >  1:2422400/1‑100 (MQ=255)
                                          caacaTCAGTTGTGACAACGTCACCGTGCGCGATGACGATCACGATCCCAATGATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgc                                                         <  2:3383747/100‑1 (MQ=255)
                                            acaTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgcca                                                       >  1:2914996/1‑100 (MQ=255)
                                             caTCAGTTGTGACAACGTCACCTTGCGTGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccag                                                      <  2:2698225/100‑1 (MQ=255)
                                             caTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccag                                                      <  2:219562/100‑1 (MQ=255)
                                             caTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccag                                                      <  2:367418/100‑1 (MQ=255)
                                             caTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccag                                                      <  2:219156/100‑1 (MQ=255)
                                             caTCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccag                                                      <  2:522235/100‑1 (MQ=255)
                                               tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccagcc                                                    <  2:3157049/100‑1 (MQ=255)
                                               tCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccagcc                                                    <  2:3344903/100‑1 (MQ=255)
                                                 aGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTgccagcc                                                    <  2:573386/98‑1 (MQ=255)
                                                 aGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCAt                                                  >  1:2132023/1‑100 (MQ=255)
                                                 aGTTGTGACAACGTCACCGTGCGCGATTACGATCACGAGCCCAAGAAGCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCAt                                                  <  2:136532/100‑1 (MQ=255)
                                                  gTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATa                                                 >  2:3196696/1‑100 (MQ=255)
                                                  gTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATa                                                 >  2:1232210/1‑100 (MQ=255)
                                                   ttGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATAc                                                <  1:4043945/100‑1 (MQ=255)
                                                    tgtgACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATAcc                                               <  2:3872901/100‑1 (MQ=255)
                                                    tgtgACAACGTCACCTTGCGCGATGACGATCACGATACCAATAATCATCGGGGCAATTACCGCATAAAGCGGCGCAACGTCCTTGATGGCCAGCCACAcc                                               <  2:13410/100‑1 (MQ=255)
                                                      tgACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAg                                             <  1:1290036/100‑1 (MQ=255)
                                                      tgACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGATTGCCAGCCATACCAg                                             >  1:684052/1‑100 (MQ=255)
                                                        aCAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCa                                           >  2:2321599/1‑100 (MQ=255)
                                                         caccGTCACCCTACGCGATGACGATCACGATCCTAATAATCAGCGCGGCAATCACCGCAAAGCGCGCCGCAACGGCCTTGGTTGCGAGCCATACCAg                                             <  2:2913390/94‑1 (MQ=255)
                                                         cAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCaa                                          >  2:3046798/1‑100 (MQ=255)
                                                         cAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCaa                                          >  2:2519969/1‑100 (MQ=255)
                                                           agggcACCTTGCGCGATGACGATCACGATCCCCATCATCCCCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCGGCAACa                                        <  2:564428/96‑1 (MQ=255)
                                                           aCGTCACCTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACa                                        <  2:135268/100‑1 (MQ=255)
                                                             gTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACAt                                       >  2:1929180/1‑99 (MQ=255)
                                                                 ccTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTc                                  >  2:197185/1‑100 (MQ=255)
                                                                  cTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTcc                                 <  2:4136386/100‑1 (MQ=255)
                                                                  cTTGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTcc                                 >  1:1158204/1‑100 (MQ=255)
                                                                    tGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCcac                               <  1:2675905/100‑1 (MQ=255)
                                                                    tGCGCGATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACATCcac                               >  2:1033201/1‑100 (MQ=255)
                                                                     gcgtGATGACGGTCACGCTCCCAATAATCAGCGCGGCAATTAACGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTcc                                 <  2:1217386/97‑1 (MQ=255)
                                                                       gcgATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCcacaca                            <  1:3428946/100‑1 (MQ=255)
                                                                       gcgATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCcacaca                            >  1:3427164/1‑100 (MQ=255)
                                                                       gcgATGACGATCACGATCCCAATAATCATCACGGCAATTCCCGCATAGCGCGGCGAAACGGCATTGGTTGCCAACCATACCAGCCACATACTTCcacaca                            >  2:4004786/1‑100 (MQ=255)
                                                                       gcgATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCcacaca                            <  1:933110/100‑1 (MQ=255)
                                                                        cgATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACaa                           >  2:733157/1‑100 (MQ=255)
                                                                        cgATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACaa                           <  2:3901058/100‑1 (MQ=255)
                                                                         gATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAAc                          <  1:2993830/100‑1 (MQ=255)
                                                                         gATGACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAAc                          >  2:2349296/1‑100 (MQ=255)
                                                                            gACGATCACGATCCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTa                       <  1:109511/100‑1 (MQ=255)
                                                                               gatcacgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAAt                    >  1:2338735/1‑100 (MQ=255)
                                                                                atcacgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAAtt                   <  1:1962688/100‑1 (MQ=255)
                                                                                 tcacgatcCCAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTg                  <  2:3266640/100‑1 (MQ=255)
                                                                                   acgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGAc                <  2:2775711/100‑1 (MQ=255)
                                                                                       tcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATCGACCGTc            <  2:519447/100‑1 (MQ=255)
                                                                                        cccAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTcc           <  2:3972393/100‑1 (MQ=255)
                                                                                         ccAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCa          >  1:2097458/1‑100 (MQ=255)
                                                                                         ccAATAAGCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCa          <  2:4092500/100‑1 (MQ=255)
                                                                                          cAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAg         >  1:2600428/1‑100 (MQ=255)
                                                                                          cAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAg         >  1:1160551/1‑100 (MQ=255)
                                                                                          cAATAATCACCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAg         <  1:79314/100‑1 (MQ=255)
                                                                                                 natcGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCa  >  1:3109312/2‑100 (MQ=255)
                                                                                                   |                                                                                                 
AAAAATCAGGGGTAATATAAGTGGGGAGAACGGGTTTAAAGACAACATCAGTTGTGACAACGTCACCTTGCGCGATGACGATCACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCA  >  NC_000913/1505462‑1505658

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: