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.

ACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  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)
            gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg‑c‑g                                                                                     <  1:2544065/100‑1 (MQ=255)
              gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑Gcc                                                                                   >  2:2191099/1‑100 (MQ=255)
               aTCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑Gcc                                                                                   <  2:1522049/99‑1 (MQ=255)
                tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAg                                                                                 <  1:4213100/100‑1 (MQ=255)
                tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAg                                                                                 <  1:1830563/100‑1 (MQ=255)
                tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAg                                                                                 <  1:756480/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGc                                                                                <  2:1041344/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAg                                                                                 <  2:408929/99‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAGc                                                                                <  2:766870/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAGc                                                                                <  2:3245599/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAGc                                                                                <  2:3348482/100‑1 (MQ=255)
                 cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCC‑GCCAGc                                                                                <  2:1695878/100‑1 (MQ=255)
                   aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCgg                                                                              <  2:846545/100‑1 (MQ=255)
                   aTTTACAAGAAGGCATCTTACCGTTTGCCGGTGGGCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCGGGCGGCGAAATTTCG‑C‑GCCAGCgg                                                                              <  1:2466661/100‑1 (MQ=255)
                     ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTc                                                                            >  1:3471136/1‑100 (MQ=255)
                     ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTc                                                                            <  2:4358673/100‑1 (MQ=255)
                       aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAg                                                                          >  1:3653673/1‑100 (MQ=255)
                        cAAGGAGGCAGCTTACCTTTTGCCGGTGGTCGACATACTGCGGATGATGTAGGTAAGCGCGCCGGGAAGGTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGc                                                                         <  1:903642/100‑1 (MQ=255)
                        cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAg                                                                          >  2:1832087/1‑99 (MQ=255)
                        cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGc                                                                         >  2:2009411/1‑100 (MQ=255)
                         tagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTTATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGcc                                                                        <  1:2438853/99‑1 (MQ=255)
                          agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGccc                                                                       >  2:1933753/1‑100 (MQ=255)
                          agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGcca                                                                       >  2:2438864/1‑99 (MQ=255)
                            aagGCATCTTAGCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGccc                                                                       >  1:781048/1‑98 (MQ=255)
                             agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCgtg                                                                    >  2:1175129/1‑100 (MQ=255)
                             agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCgtg                                                                    <  2:884544/100‑1 (MQ=255)
                              ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCgtg                                                                    <  1:2935622/99‑1 (MQ=255)
                               gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGtt                                                                  <  2:4393484/100‑1 (MQ=255)
                               gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGtt                                                                  <  2:1466537/100‑1 (MQ=255)
                                  tCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCt                                                               >  2:2260730/1‑100 (MQ=255)
                                   cTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTg                                                              >  1:1961521/1‑100 (MQ=255)
                                    ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTgg                                                             <  2:3726350/100‑1 (MQ=255)
                                    ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTgg                                                             >  2:2674854/1‑100 (MQ=255)
                                    ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTgg                                                             >  2:1095252/1‑100 (MQ=255)
                                     tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGt                                                            <  2:1994290/100‑1 (MQ=255)
                                     tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGt                                                            <  2:828481/100‑1 (MQ=255)
                                     tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGt                                                            >  2:2819308/1‑100 (MQ=255)
                                      aCCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTg                                                           >  2:1653051/1‑100 (MQ=255)
                                       ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAACCCGTGGTTCTGGTgc                                                          >  1:3615427/1‑100 (MQ=255)
                                        cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                         >  2:584876/1‑100 (MQ=255)
                                        cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                         <  1:2191334/100‑1 (MQ=255)
                                           ttGCCGGTGGTCCACATACTTCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGcgcg                                                       <  1:1652883/99‑1 (MQ=255)
                                            tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     >  1:1894789/1‑100 (MQ=255)
                                            tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     <  2:3695607/100‑1 (MQ=255)
                                            tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                     <  2:4185219/100‑1 (MQ=255)
                                             gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    <  1:3001699/100‑1 (MQ=255)
                                             gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    >  2:3955030/1‑100 (MQ=255)
                                             gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTg                                                    >  2:3865085/1‑100 (MQ=255)
                                                 gtggtCCATATACTGCCGATGAGGTAGTTAATCGCGCCGGGCAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGCTTATGGAGCGCGTTGGCTg                                                <  1:650944/100‑1 (MQ=255)
                                                   ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                              >  2:4114778/1‑100 (MQ=255)
                                                   ggtCCACATACTGCCGATGAGGTAGTTAATCGCACCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                              <  2:4267002/100‑1 (MQ=255)
                                                    gtCCAGATACGGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTGCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             <  1:3752640/100‑1 (MQ=255)
                                                    gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                              <  1:986419/99‑1 (MQ=255)
                                                    gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                             <  2:2397025/100‑1 (MQ=255)
                                                      ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                            >  1:1063038/1‑99 (MQ=255)
                                                      ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                           <  1:1183544/100‑1 (MQ=255)
                                                      ccACATACTGAAAATGAGGAAGTGAATTGCGCATGTAAGCTGGCGGCGAAATTATT‑T‑GACAGCGGTAAGTCCGTGTTTCTGGTGCGCGTTgg                                                   >  1:1253156/1‑92 (MQ=255)
                                                        acaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGcc                                         >  1:1048911/1‑100 (MQ=255)
                                                         caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                        <  2:1430869/100‑1 (MQ=255)
                                                          aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                       <  2:2153765/100‑1 (MQ=255)
                                                          aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                       <  2:3554328/100‑1 (MQ=255)
                                                             cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    >  2:1154403/1‑100 (MQ=255)
                                                             cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                    <  1:2259067/100‑1 (MQ=255)
                                                              tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAc                                   <  1:1095249/100‑1 (MQ=255)
                                                               gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                  <  1:2040527/100‑1 (MQ=255)
                                                                ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 <  1:4404098/100‑1 (MQ=255)
                                                                  gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTg                               <  1:533602/100‑1 (MQ=255)
                                                                  gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTCGTGCGCGTTGGCTGGTCTGCCGATTACGc                                 >  1:3511058/1‑98 (MQ=255)
                                                                    tGAGGTAGTTAATCGCGCCGGGATGCTGGCGGGGAAAGTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATGACGc                                 <  1:2669959/96‑1 (MQ=255)
                                                                           gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      >  1:2497930/1‑100 (MQ=255)
                                                                           gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      <  1:1426932/100‑1 (MQ=255)
                                                                           gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                      >  2:477539/1‑100 (MQ=255)
                                                                             tAATCGCGCCGGGAAGCTTGCGGTGAAATTTCG‑T‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                     <  1:311298/99‑1 (MQ=255)
                                                                               aTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                  >  2:3810692/1‑100 (MQ=255)
                                                                                  gcgcCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGcc               >  2:2883601/1‑100 (MQ=255)
                                                                                   cgcCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCg              <  2:2417541/100‑1 (MQ=255)
                                                                                       gggAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg          >  1:3962198/1‑100 (MQ=255)
                                                                                         gAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt        <  1:2764834/100‑1 (MQ=255)
                                                                                         accgCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCCTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt        >  1:2354835/4‑100 (MQ=38)
                                                                                          ccgCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt        >  1:833385/3‑99 (MQ=38)
                                                                                          ccgCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:3844558/3‑100 (MQ=39)
                                                                                          ccgCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg       >  2:2390300/3‑100 (MQ=39)
                                                                                           cgCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc      >  2:2546175/2‑100 (MQ=39)
                                                                                           cgCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc      >  2:4240833/2‑100 (MQ=39)
                                                                                           cgCTGGC‑GCGAAATTTCG‑CTGCCAGCGGTCAGCCCGTGTTACAGGTGCGCGTTGGATGGTCTGCCGCTTACGCTGAAGCTTTAAACCAGCTGAGCGATg       >  1:2434324/2‑99 (MQ=255)
                                                                                            gCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg     <  2:1790706/100‑1 (MQ=255)
                                                                                            gCTGGC‑GCGAAATTTCGCC‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg     >  1:4352362/1‑100 (MQ=255)
                                                                                              tGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc   >  2:3283864/1‑100 (MQ=255)
                                                                                              tGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc   <  2:2464244/100‑1 (MQ=255)
                                                                                              tGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc   <  2:1346172/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
ACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑C‑GCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCC  >  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: ‑