Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,302,596 A→G intergenic (+183/+532) eco → / ← yojH ecotin/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,302,5960AG92.0% 177.9 / 10.5 87intergenic (+183/+532)eco/yojHecotin/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/7);  new base G (33/47);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 4.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

GCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAC  >  NC_000913/2302509‑2302684
                                                                                       |                                                                                        
gcgTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCgg                                                                                     <  2:1692037/93‑1 (MQ=255)
 cgTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGt                                                                                    >  1:1039329/1‑93 (MQ=255)
 cgTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGt                                                                                    <  2:1557057/93‑1 (MQ=255)
  gTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGtt                                                                                   >  1:1449879/1‑93 (MQ=255)
  gTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGtt                                                                                   <  1:2137901/93‑1 (MQ=255)
  gTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGtt                                                                                   <  1:1936265/93‑1 (MQ=255)
     ttATCAGGCCTACACCGCTGTGAAGTGCGGGACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGGTGcc                                                                                >  1:1779520/1‑93 (MQ=255)
      tATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCt                                                                               >  2:1511017/1‑93 (MQ=255)
      tATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCt                                                                               <  1:41309/93‑1 (MQ=255)
        tCATGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGa                                                                             <  2:901785/93‑1 (MQ=255)
        tCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCtca                                                                             <  2:3733461/93‑3 (MQ=255)
        tCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGa                                                                             <  2:851271/93‑1 (MQ=255)
        tCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGa                                                                             >  2:3890817/1‑93 (MQ=255)
         cAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGAt                                                                            <  1:4035965/93‑1 (MQ=255)
         cAGCCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGAt                                                                            <  1:2044494/93‑1 (MQ=255)
          aGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATg                                                                           <  2:1165116/93‑1 (MQ=255)
           ggCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGc                                                                          >  2:3562427/1‑93 (MQ=255)
             ccTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAGGGCAGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGa                                                                        <  1:2511710/93‑1 (MQ=255)
             ccTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGa                                                                        <  1:2306907/93‑1 (MQ=255)
              cTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGAc                                                                       >  2:580121/1‑93 (MQ=255)
              cTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGAc                                                                       >  1:3499258/1‑93 (MQ=255)
               tACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACg                                                                      <  2:3203715/93‑1 (MQ=255)
                acacCGCCGTGAAGTGCGTCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGc                                                                     <  1:3812618/93‑1 (MQ=11)
                 cacCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCt                                                                    <  2:2730678/93‑1 (MQ=255)
                  acCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTg                                                                   <  2:3864386/93‑1 (MQ=255)
                     gCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGgcg                                                                >  1:273515/1‑93 (MQ=18)
                      cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGgcgc                                                               >  1:1186583/1‑93 (MQ=17)
                       tgtgAAGTGCGGCACCCCGTAGGTCAGATAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGcgcg                                                              >  2:4534336/1‑93 (MQ=9)
                       tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGcgcg                                                              <  1:2895702/93‑1 (MQ=14)
                       tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGcgcg                                                              <  2:1776667/93‑1 (MQ=14)
                       tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGcgcg                                                              <  2:3984917/93‑1 (MQ=14)
                       tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGcgcg                                                              >  2:3146199/1‑93 (MQ=14)
                       tgtgAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGcgcg                                                              >  1:309458/1‑93 (MQ=14)
                                                      ggCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGGGCGGCTTATCAAGCCTACACCGCTGGGAAGTGCt                               >  2:1227735/1‑93 (MQ=11)
                                                       gCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTc                              <  1:3609586/93‑1 (MQ=17)
                                                       gCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTc                              <  1:1449392/93‑1 (MQ=17)
                                                       gCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTc                              <  2:1678525/93‑1 (MQ=17)
                                                         ggTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCa                            <  2:3219334/93‑1 (MQ=21)
                                                          gTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAc                           >  2:1610753/1‑93 (MQ=25)
                                                          gTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAc                           >  2:387952/1‑93 (MQ=25)
                                                           tCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTc                              >  2:1044321/1‑89 (MQ=17)
                                                            cACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAccc                         >  1:4032983/1‑93 (MQ=25)
                                                             ccgccgcATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCcccccc                        <  2:3227396/92‑1 (MQ=255)
                                                             aCGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                        <  1:4339695/93‑1 (MQ=33)
                                                              cgccgcATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCg                       >  1:74011/1‑93 (MQ=18)
                                                                ccgcATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                     <  1:2128592/93‑1 (MQ=18)
                                                                ccgcATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                     <  2:3109266/93‑1 (MQ=18)
                                                                  gcATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                   >  1:2609031/1‑93 (MQ=21)
                                                                  gcATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                   >  1:2615791/1‑93 (MQ=18)
                                                                  gcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                   <  2:3499258/93‑1 (MQ=33)
                                                                   cATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                  >  2:2202011/1‑93 (MQ=21)
                                                                    aTCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                 <  2:3701229/93‑1 (MQ=18)
                                                                    aTCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                 <  1:3948700/93‑1 (MQ=18)
                                                                    aTCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                 >  1:3904406/1‑93 (MQ=18)
                                                                    aTCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                 <  2:1785461/93‑1 (MQ=18)
                                                                      ccGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg               <  2:4219586/93‑1 (MQ=18)
                                                                       cGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa              >  2:2239031/1‑93 (MQ=18)
                                                                        gACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt             >  2:1299512/1‑93 (MQ=18)
                                                                        gACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt             <  2:3031153/93‑1 (MQ=18)
                                                                         aCATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa            <  1:633678/93‑1 (MQ=18)
                                                                         aCATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa            <  1:4605679/93‑1 (MQ=18)
                                                                         aCATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa            <  1:4038471/93‑1 (MQ=18)
                                                                         aCATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa            >  2:183723/1‑93 (MQ=21)
                                                                          cATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa           <  2:3172672/93‑1 (MQ=18)
                                                                          cATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa           >  1:104391/1‑93 (MQ=21)
                                                                          cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa           <  1:1287736/93‑1 (MQ=33)
                                                                          cATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa           <  1:4354656/93‑1 (MQ=33)
                                                                           aTCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCGTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTCGGATAAg          <  1:946097/93‑1 (MQ=255)
                                                                           aTCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAg          <  1:1545900/93‑1 (MQ=18)
                                                                           aTCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAg          >  1:2130477/1‑93 (MQ=18)
                                                                           aTCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAg          <  2:3254001/93‑1 (MQ=18)
                                                                            tCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg         <  1:2971326/93‑1 (MQ=18)
                                                                            tCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg         <  2:1658470/93‑1 (MQ=18)
                                                                            tCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg         >  1:722629/1‑93 (MQ=21)
                                                                            tCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg         >  1:655765/1‑93 (MQ=18)
                                                                              cAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg       <  1:3710758/93‑1 (MQ=18)
                                                                              cAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg       <  2:637946/93‑1 (MQ=18)
                                                                               aaCGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcc                          >  2:4214523/1‑73 (MQ=19)
                                                                               aaCGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg      <  1:183723/93‑1 (MQ=21)
                                                                               aaCGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg      >  2:2101823/1‑93 (MQ=18)
                                                                                aCGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                     >  1:1524623/1‑77 (MQ=19)
                                                                                aCGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGt     >  2:2086779/1‑93 (MQ=9)
                                                                                aCGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt     >  2:119860/1‑93 (MQ=18)
                                                                                aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                     <  2:1524623/77‑1 (MQ=34)
                                                                                  gCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa   <  1:3942987/93‑1 (MQ=17)
                                                                                  gCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa   <  1:3972660/93‑1 (MQ=17)
                                                                                  gCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa   <  2:1780810/93‑1 (MQ=17)
                                                                                  gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa   <  1:809285/93‑1 (MQ=32)
                                                                                  gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa   <  2:2814151/93‑1 (MQ=32)
                                                                                   cccGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc  <  1:2671958/93‑1 (MQ=14)
                                                                                   cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc  <  1:1726781/93‑1 (MQ=31)
                                                                                       |                                                                                        
GCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAC  >  NC_000913/2302509‑2302684

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