Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,955 Δ1 bp 17.3% coding (124/567 nt) yecD → isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,9540G.17.3% 536.5 / 113.7 185coding (123/567 nt)yecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base G (68/85);  new base . (20/12);  total (88/97)
Fisher's exact test for biased strand distribution p-value = 7.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 194 total aligned reads displayed.

ACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  NC_000913/1950856‑1951050
                                                                                                  |                                                                                                  
aCTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcgg                                                                                                   >  2:2092236/1‑100 (MQ=255)
 cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                  >  1:3161406/1‑100 (MQ=255)
 cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                  <  2:2628009/100‑1 (MQ=255)
     gttggtggTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGaaa                                                                                              <  1:2100816/98‑1 (MQ=255)
     gcTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCACCGGGAAGCTGGCGGCGaaa                                                                                              >  1:4266970/1‑100 (MQ=255)
      cTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAt                                                                                             <  1:1513074/100‑1 (MQ=255)
        ggtggtGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAttt                                                                                           >  2:2710516/1‑100 (MQ=255)
         gtggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTc                                                                                          <  1:3878486/100‑1 (MQ=255)
         ggtgtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAtt                                                                                            <  1:510277/95‑1 (MQ=255)
            gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCGGGGAAGCTGGCGGCGAAAtttgg                                                                                         >  1:4304262/1‑96 (MQ=255)
            gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgcg                                                                                     <  1:2544065/100‑1 (MQ=255)
              gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                   >  2:2191099/1‑100 (MQ=255)
               aTCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                   <  2:1522049/99‑1 (MQ=255)
                tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAg                                                                                 <  1:4213100/100‑1 (MQ=255)
                tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAg                                                                                 <  1:1830563/100‑1 (MQ=255)
                tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAg                                                                                 <  1:756480/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGc                                                                                <  2:1041344/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAg                                                                                 <  2:408929/99‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:766870/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:3245599/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:3348482/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                <  2:1695878/100‑1 (MQ=255)
                   aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCgg                                                                              <  2:846545/100‑1 (MQ=255)
                   aTTTACAAGAAGGCATCTTACCGTTTGCCGGTGGGCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCGGGCGGCGAAATTTCGCGCCAGCgg                                                                              <  1:2466661/100‑1 (MQ=255)
                     ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTc                                                                            >  1:3471136/1‑100 (MQ=255)
                     ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTc                                                                            <  2:4358673/100‑1 (MQ=255)
                       aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAg                                                                          >  1:3653673/1‑100 (MQ=255)
                        cAAGGAGGCAGCTTACCTTTTGCCGGTGGTCGACATACTGCGGATGATGTAGGTAAGCGCGCCGGGAAGGTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                         <  1:903642/100‑1 (MQ=255)
                        cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAg                                                                          >  2:1832087/1‑99 (MQ=255)
                        cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                         >  2:2009411/1‑100 (MQ=255)
                         tagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTTATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                        <  1:2438853/99‑1 (MQ=255)
                          agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGccc                                                                       >  2:1933753/1‑100 (MQ=255)
                          agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcca                                                                       >  2:2438864/1‑99 (MQ=255)
                            aagGCATCTTAGCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGccc                                                                       >  1:781048/1‑98 (MQ=255)
                             agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                    >  2:1175129/1‑100 (MQ=255)
                             agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                    <  2:884544/100‑1 (MQ=255)
                              ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                    <  1:2935622/99‑1 (MQ=255)
                               gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGtt                                                                  <  2:4393484/100‑1 (MQ=255)
                               gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGtt                                                                  <  2:1466537/100‑1 (MQ=255)
                                  tCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCt                                                               >  2:2260730/1‑100 (MQ=255)
                                   cTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTg                                                              >  1:1961521/1‑100 (MQ=255)
                                    ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTgg                                                             <  2:3726350/100‑1 (MQ=255)
                                    ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTgg                                                             >  2:2674854/1‑100 (MQ=255)
                                    ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTgg                                                             >  2:1095252/1‑100 (MQ=255)
                                     tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGt                                                            <  2:1994290/100‑1 (MQ=255)
                                     tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGt                                                            <  2:828481/100‑1 (MQ=255)
                                     tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGt                                                            >  2:2819308/1‑100 (MQ=255)
                                      aCCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTg                                                           >  2:1653051/1‑100 (MQ=255)
                                       ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAACCCGTGGTTCTGGTgc                                                          >  1:3615427/1‑100 (MQ=255)
                                        cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                         >  2:584876/1‑100 (MQ=255)
                                        cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                         <  1:2191334/100‑1 (MQ=255)
                                           ttGCCGGTGGTCCACATACTTCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGcgcg                                                       <  1:1652883/99‑1 (MQ=255)
                                            tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     >  1:1894789/1‑100 (MQ=255)
                                            tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     <  2:3695607/100‑1 (MQ=255)
                                            tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     <  2:4185219/100‑1 (MQ=255)
                                             gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    <  1:3001699/100‑1 (MQ=255)
                                             gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    >  2:3955030/1‑100 (MQ=255)
                                             gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    >  2:3865085/1‑100 (MQ=255)
                                                 gtggtCCATATACTGCCGATGAGGTAGTTAATCGCGCCGGGCAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGCTTATGGAGCGCGTTGGCTg                                                <  1:650944/100‑1 (MQ=255)
                                                   ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                              >  2:4114778/1‑100 (MQ=255)
                                                   ggtCCACATACTGCCGATGAGGTAGTTAATCGCACCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                              <  2:4267002/100‑1 (MQ=255)
                                                    gtCCAGATACGGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTGCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             <  1:3752640/100‑1 (MQ=255)
                                                    gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                              <  1:986419/99‑1 (MQ=255)
                                                    gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             <  2:2397025/100‑1 (MQ=255)
                                                      ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                            >  1:1063038/1‑99 (MQ=255)
                                                      ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                           <  1:1183544/100‑1 (MQ=255)
                                                      ccACATACTGAAAATGAGGAAGTGAATTGCGCATGTAAGCTGGCGGCGAAATTATTTGACAGCGGTAAGTCCGTGTTTCTGGTGCGCGTTgg                                                   >  1:1253156/1‑92 (MQ=255)
                                                        acaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGcc                                         >  1:1048911/1‑100 (MQ=255)
                                                         caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                        <  2:1430869/100‑1 (MQ=255)
                                                          aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                       <  2:2153765/100‑1 (MQ=255)
                                                          aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                       <  2:3554328/100‑1 (MQ=255)
                                                             cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    >  2:1154403/1‑100 (MQ=255)
                                                             cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    <  1:2259067/100‑1 (MQ=255)
                                                              tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAc                                   <  1:1095249/100‑1 (MQ=255)
                                                               gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  <  1:2040527/100‑1 (MQ=255)
                                                                ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 <  1:4404098/100‑1 (MQ=255)
                                                                  gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTg                               <  1:533602/100‑1 (MQ=255)
                                                                  gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTCGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 >  1:3511058/1‑98 (MQ=255)
                                                                    tGAGGTAGTTAATCGCGCCGGGATGCTGGCGGGGAAAGTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATGACGc                                 <  1:2669959/96‑1 (MQ=255)
                                                                           gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      >  1:2497930/1‑100 (MQ=255)
                                                                           gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      <  1:1426932/100‑1 (MQ=255)
                                                                           gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      >  2:477539/1‑100 (MQ=255)
                                                                             tAATCGCGCCGGGAAGCTTGCGGTGAAATTTCGTGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  1:311298/99‑1 (MQ=255)
                                                                               aTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                  >  2:3810692/1‑100 (MQ=255)
                                                                                  gcgcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGcc               >  2:2883601/1‑100 (MQ=255)
                                                                                   cgcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCg              <  2:2417541/100‑1 (MQ=255)
                                                                                       gggAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg          >  1:3962198/1‑100 (MQ=255)
                                                                                         gAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt        <  1:2764834/100‑1 (MQ=255)
                                                                                         accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCCTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt        >  1:2354835/4‑100 (MQ=38)
                                                                                          ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt        >  1:833385/3‑99 (MQ=38)
                                                                                          ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:3844558/3‑100 (MQ=39)
                                                                                          ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2390300/3‑100 (MQ=39)
                                                                                           cgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc      >  2:2546175/2‑100 (MQ=39)
                                                                                           cgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc      >  2:4240833/2‑100 (MQ=39)
                                                                                           cgCTGGCGCGAAATTTCGCTGCCAGCGGTCAGCCCGTGTTACAGGTGCGCGTTGGATGGTCTGCCGCTTACGCTGAAGCTTTAAACCAGCTGAGCGATg       >  1:2434324/2‑99 (MQ=255)
                                                                                            gCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg     <  2:1790706/100‑1 (MQ=255)
                                                                                            gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg     >  1:4352362/1‑100 (MQ=255)
                                                                                              tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc   >  2:3283864/1‑100 (MQ=255)
                                                                                              tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc   <  2:2464244/100‑1 (MQ=255)
                                                                                              tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc   <  2:1346172/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
ACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  NC_000913/1950856‑1951050

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: