Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,955 Δ1 bp 16.7% 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.16.7% 380.8 / 77.1 132coding (123/567 nt)yecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base G (51/59);  new base . (16/6);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 3.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 143 total aligned reads displayed.

CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCA  >  NC_000913/1950855‑1951053
                                                                                                   |                                                                                                    
cACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcg                                                                                                      <  1:477706/100‑1 (MQ=255)
 aCTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcgg                                                                                                     <  1:771265/100‑1 (MQ=255)
 aCTGCGCGGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcgg                                                                                                     <  2:327419/100‑1 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcggc                                                                                                    >  1:2927646/1‑100 (MQ=255)
  cTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTAGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTggcagc                                                                                                    >  2:285554/1‑100 (MQ=255)
   tGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGgcggcg                                                                                                   <  2:47873/100‑1 (MQ=255)
   tGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGgcggcg                                                                                                   <  2:2146299/100‑1 (MQ=255)
   tGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGgcggcg                                                                                                   <  1:2185193/100‑1 (MQ=255)
    gcgcTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGa                                                                                                  >  1:695270/1‑100 (MQ=255)
     cgcTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGaa                                                                                                 <  1:677702/100‑1 (MQ=255)
        tggtggTGATCGATTTACAAGTAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAATTAATCGCGCCGGGAAGCTGGCGGCGAAAtt                                                                                              <  2:1220899/100‑1 (MQ=255)
         ggtggtGAACGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAAtgt                                                                                             >  2:3409114/1‑98 (MQ=255)
           tggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg                                                                                           <  1:2168936/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgc                                                                                         <  1:488175/100‑1 (MQ=255)
            ggtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgc                                                                                         <  1:3013351/100‑1 (MQ=255)
             gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgcg                                                                                        <  1:51878/100‑1 (MQ=255)
             gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgcg                                                                                        >  1:3213414/1‑100 (MQ=255)
             gtgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcgcg                                                                                        <  2:2692594/100‑1 (MQ=255)
              tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCgcgc                                                                                       <  2:2604132/100‑1 (MQ=255)
              tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCgcgc                                                                                       <  1:2053738/100‑1 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTcg                                                                                           <  2:3142088/96‑1 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                      >  1:2532849/1‑100 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                      >  1:2688042/1‑100 (MQ=255)
               gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                      <  2:1388084/100‑1 (MQ=255)
                aTCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTcgccgcca                                                                                    <  2:1370575/100‑6 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAg                                                                                    <  1:3639091/100‑1 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAg                                                                                    <  2:2742287/100‑1 (MQ=255)
                 tCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAg                                                                                    <  1:2157209/100‑1 (MQ=255)
                  cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAg                                                                                    >  2:3085433/1‑99 (MQ=255)
                    aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCgg                                                                                 <  1:1265457/100‑1 (MQ=255)
                     tttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGcttc                                                                                <  1:2215176/100‑4 (MQ=38)
                      ttACAAGAAGGCGTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTc                                                                               >  1:2396829/1‑100 (MQ=255)
                       tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCa                                                                              >  2:294003/1‑100 (MQ=255)
                       tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCAACAAAATAACGAAGAGGTAGTTAATCGCGACGGGAAGCTGGCGGCGAAATTTCGCGACAGAGGt                                                                                >  2:240340/1‑98 (MQ=255)
                        aCAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAg                                                                             >  2:789372/1‑100 (MQ=255)
                              agGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                       >  2:145327/1‑100 (MQ=255)
                               ggCCTCTTTCCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                       <  2:2329471/99‑1 (MQ=255)
                               gacATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtgt                                                                      >  1:2452983/3‑100 (MQ=255)
                                 cATCTTAGCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtg                                                                       <  2:2331769/97‑1 (MQ=255)
                                 cATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGttt                                                                    <  1:2820426/100‑1 (MQ=255)
                                 cATCTTACCTTTTGCCGGTGGTCCACATACTGCAGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCAGTGttt                                                                    >  2:203775/1‑100 (MQ=255)
                                    cTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTg                                                                 >  1:3521595/1‑100 (MQ=255)
                                    cTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTg                                                                 >  1:1227871/1‑100 (MQ=255)
                                     ttACCTTTTGCCGGTGGTCCTCATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTgg                                                                <  2:2093882/100‑1 (MQ=255)
                                     ttACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTgg                                                                <  1:789377/100‑1 (MQ=255)
                                      tACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGt                                                               >  2:3120785/1‑100 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACATACTGCTGATGAGGTCGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCCCCAGCGGTCAGCCCGTGTATCTGGTgc                                                             >  2:741394/1‑100 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                             <  1:1000491/100‑1 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                             >  2:446368/1‑100 (MQ=255)
                                        ccTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                             <  1:799075/100‑1 (MQ=255)
                                         cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcg                                                            >  2:3309254/1‑100 (MQ=255)
                                                cGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGc                                                     <  1:107646/100‑1 (MQ=255)
                                                 ggtggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCt                                                    <  1:3363132/100‑1 (MQ=255)
                                                 ggtggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCt                                                    <  1:3431741/100‑1 (MQ=255)
                                                  gtggtCGACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCATGTTTCTGGTGCGCGTTGGCTg                                                   <  2:951671/100‑1 (MQ=255)
                                                   tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                  <  1:2914538/100‑1 (MQ=255)
                                                   tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                  <  2:176139/100‑1 (MQ=255)
                                                    ggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGt                                                 >  2:1242140/1‑100 (MQ=255)
                                                     gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                >  1:2343708/1‑100 (MQ=255)
                                                     gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                >  1:3343662/1‑100 (MQ=255)
                                                      tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                               >  1:132640/1‑100 (MQ=255)
                                                        cacaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGc                                             >  2:3454790/1‑100 (MQ=255)
                                                          caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                           <  1:3018208/100‑1 (MQ=255)
                                                           aTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                          <  2:1005171/100‑1 (MQ=255)
                                                               tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAc                                      >  2:276295/1‑100 (MQ=255)
                                                                 ccGATGAGGTAGTTAATCGCTCCGGGAAGCTGTGGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAtt                                        <  2:3211091/96‑1 (MQ=255)
                                                                 ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                    <  1:209738/100‑1 (MQ=255)
                                                                 ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                    >  1:3230571/1‑100 (MQ=255)
                                                                     tGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGaa                                >  2:2059445/1‑100 (MQ=255)
                                                                      gAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAg                               <  1:1123137/100‑1 (MQ=255)
                                                                       aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTTTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                              <  2:1956136/100‑1 (MQ=255)
                                                                       aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                              <  1:1517925/100‑1 (MQ=255)
                                                                        ggTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCa                             >  1:699906/1‑100 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGTCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                          <  2:462691/100‑1 (MQ=255)
                                                                           aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCCGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                          >  2:545443/1‑100 (MQ=255)
                                                                            gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaa                         >  2:1910024/1‑100 (MQ=255)
                                                                             ttAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaa                        <  1:2102332/100‑1 (MQ=255)
                                                                              tAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGGTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTaaaa                       >  2:1209970/1‑100 (MQ=255)
                                                                                   gcgcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGcc                  >  2:1573385/1‑100 (MQ=255)
                                                                                      gctgaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGt               >  1:2734730/8‑100 (MQ=37)
                                                                                       ctgaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc              >  1:1094364/7‑100 (MQ=38)
                                                                                         gaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGa            >  2:511399/5‑100 (MQ=38)
                                                                                          gAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt           >  1:504204/1‑100 (MQ=255)
                                                                                           ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          >  1:1620190/3‑100 (MQ=39)
                                                                                           aaGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          <  1:2710726/100‑1 (MQ=255)
                                                                                           aaGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg          <  1:3235732/100‑1 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  1:2756589/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  2:126057/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  1:1622442/1‑100 (MQ=255)
                                                                                             gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg        >  1:287904/1‑100 (MQ=255)
                                                                                              cTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc       >  2:922761/1‑100 (MQ=255)
                                                                                              cTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc       >  2:1353773/1‑100 (MQ=255)
                                                                                               tGGCGGTGAAATTTCGCGCCAGTGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcgc       <  1:259663/99‑1 (MQ=255)
                                                                                               tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc      >  1:281513/1‑100 (MQ=255)
                                                                                               tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc      <  1:956347/100‑1 (MQ=255)
                                                                                               tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc      >  1:367900/1‑100 (MQ=255)
                                                                                               tGGCGGCGAAATTCCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc      <  1:1910034/100‑1 (MQ=255)
                                                                                                 gcggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCt    >  1:3630925/1‑100 (MQ=255)
                                                                                                   ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa  >  1:2576192/1‑100 (MQ=255)
                                                                                                   ggcgAAATTGCGCGCCAGCGGTCAGCTCTTGTTTCTGGTGCTCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa  <  2:2047336/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
CACTGCGCTGGTGGTGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCA  >  NC_000913/1950855‑1951053

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: