Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,105,120 C→G 8.9% A174A (GCC→GCG) ‡ hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,105,1200CG 8.9% 201.1 / 9.1 67A174A (GCC→GCG) ‡hisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base G (4/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

GAACAAGGTCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGCCCGTCT  >  NC_002947/1104974‑1105268
                                                                                                                                                  |                                                                                                                                                     
gAACAAGGTCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCAc                                                                                                                                                     >  8:347904/1‑149 (MQ=255)
    aaGGTCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATc                                                                                                                                                 <  4:147411/149‑1 (MQ=255)
      ggTCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAg                                                                                                                                               <  5:306710/149‑1 (MQ=255)
       gTCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGTGGCGATCAgt                                                                                                                                              <  7:271867/149‑2 (MQ=255)
       gTCGCCAGATCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGc                                                                                                                                              >  4:97714/1‑149 (MQ=255)
        tCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGct                                                                                                                                             >  1:148284/1‑149 (MQ=255)
          gCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGctc                                                                                                                                            <  3:282508/148‑1 (MQ=255)
                tCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCtg                                                                                                                                     <  4:188406/149‑1 (MQ=255)
                   aCATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCtggtg                                                                                                                                  >  5:100247/1‑149 (MQ=255)
                    catcatCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCtggtgg                                                                                                                                 >  6:185563/1‑149 (MQ=255)
                      tcatcaAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTc                                                                                                                               <  3:266971/149‑1 (MQ=255)
                            aGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACaaa                                                                                                                         <  7:150304/149‑1 (MQ=255)
                                gTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCt                                                                                                                     >  6:80475/1‑149 (MQ=255)
                                   tGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCa                                                                                                                  <  6:100247/149‑1 (MQ=255)
                                   tGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCa                                                                                                                  >  5:327711/1‑149 (MQ=255)
                                   tGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCa                                                                                                                  <  8:16170/149‑1 (MQ=255)
                                       tCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGa                                                                                                               >  7:342513/1‑148 (MQ=255)
                                         gATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGaaa                                                                                                             <  1:287407/148‑1 (MQ=255)
                                           tGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATg                                                                                                          <  6:73809/149‑1 (MQ=255)
                                             gAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATg                                                                                                          >  8:105063/1‑147 (MQ=255)
                                                   gcgcCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGc                                                                                                  >  2:11515/1‑149 (MQ=255)
                                                     gccgcTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGTGGCGATCAGttcctg                                                                                                                                         >  1:22596/1‑105 (MQ=255)
                                                     gccgcTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGTGGCGATCAGttcctg                                                                                                                                         <  2:22596/108‑4 (MQ=255)
                                                     gccgcTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGccc                                                                                                <  5:185563/149‑1 (MQ=255)
                                                      ccgcTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCg                                                                                               >  3:210470/1‑149 (MQ=255)
                                                             tCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCag                                                                                        >  8:80657/1‑149 (MQ=255)
                                                               aaCCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCagag                                                                                      >  3:304291/1‑149 (MQ=255)
                                                               aaCCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCagag                                                                                      <  3:97714/149‑1 (MQ=255)
                                                               aaCCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCaga                                                                                       <  7:29381/148‑1 (MQ=255)
                                                                      ccGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGAt                                                                                <  8:342513/148‑1 (MQ=255)
                                                                      ccGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATc                                                                               <  3:164412/149‑1 (MQ=255)
                                                                       cGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATc                                                                               >  2:371678/1‑148 (MQ=255)
                                                                                 atcGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAAAGCAGCACGCGCGTATCCAGAGCCTGATCGACACGCTGCg                                                                    >  6:214987/1‑149 (MQ=255)
                                                                                         ggTCGACACCGGCAACACGCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCGGGTGGTGAACAAAGCCTCCATGAAAATGAAGAACGCCCGTATCCAGAGCCTGATAGACACGCTGCGTGTAGCg‑g                                                           >  6:340474/1‑149 (MQ=255)
                                                                                          gTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑Gt                                                          <  1:371678/149‑1 (MQ=255)
                                                                                           tCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTc                                                         >  7:8705/1‑149 (MQ=255)
                                                                                              acacCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGaa                                                      >  3:77203/1‑149 (MQ=255)
                                                                                               cacCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGAC‑CGCTGCACGCAGCGTGTCGatc                                                     <  7:105063/149‑3 (MQ=255)
                                                                                                acCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAATTGCAGCACGCCCGTATCGAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATc                                                    >  8:347413/1‑149 (MQ=255)
                                                                                                acCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATc                                                    <  6:327711/149‑1 (MQ=255)
                                                                                                  cGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATcgc                                                  <  8:8705/149‑1 (MQ=255)
                                                                                                      aaCACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGaca                                              >  7:117317/1‑149 (MQ=255)
                                                                                                         acCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACcg                                           <  4:210470/149‑1 (MQ=255)
                                                                                                           ccTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑Gt                                                          <  1:120102/132‑1 (MQ=255)
                                                                                                           ccTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑Gt                                                          >  2:120102/1‑132 (MQ=255)
                                                                                                            cTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCgg                                        >  8:132465/1‑149 (MQ=255)
                                                                                                                    cAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTcc                                >  5:327150/1‑149 (MQ=255)
                                                                                                                    cAACGGACTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCACGGTTGACTcc                                >  4:86078/1‑149 (MQ=255)
                                                                                                                              gAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGAc                      <  7:80657/149‑1 (MQ=255)
                                                                                                                                 ccccAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGaca                                                                           <  6:146698/94‑1 (MQ=255)
                                                                                                                                 ccccAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGaca                                                                           >  5:146698/1‑94 (MQ=255)
                                                                                                                                 ccccAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGTGACCtt                   >  2:306783/1‑149 (MQ=255)
                                                                                                                                     aGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCt               >  6:219942/1‑149 (MQ=255)
                                                                                                                                     aGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCt               >  8:205896/1‑149 (MQ=255)
                                                                                                                                      agagcTGATCGTGGCGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTg              >  5:118101/5‑149 (MQ=255)
                                                                                                                                       gAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGt             >  7:218367/1‑149 (MQ=255)
                                                                                                                                        aaCTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTc            <  4:304291/149‑1 (MQ=255)
                                                                                                                                         aCTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTcg           <  1:7609/149‑1 (MQ=255)
                                                                                                                                          cTGATCGTGGCGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATcgc                                                  >  8:324827/1‑109 (MQ=255)
                                                                                                                                          cTGATCGTGGCGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTcgc          >  2:313146/1‑149 (MQ=255)
                                                                                                                                          cTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTcgc          >  7:299135/1‑149 (MQ=255)
                                                                                                                                           tGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTcgcg         >  4:159055/1‑149 (MQ=255)
                                                                                                                                             aTCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGcc       >  2:338525/1‑149 (MQ=255)
                                                                                                                                              tCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGccc      <  6:327150/149‑1 (MQ=255)
                                                                                                                                              tCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAAAGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGccc      <  5:214987/149‑1 (MQ=255)
                                                                                                                                                gCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGCCCg     <  4:77203/148‑1 (MQ=255)
                                                                                                                                                gCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGCCCGt    <  7:347904/149‑1 (MQ=255)
                                                                                                                                                 ccACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCgg                                        <  3:301895/112‑1 (MQ=255)
                                                                                                                                                 ccACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCgg                                        >  4:301895/1‑112 (MQ=255)
                                                                                                                                                  cACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGCCCGTc   >  5:253963/1‑148 (MQ=255)
                                                                                                                                                  cACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGCCCGTCt  <  5:219942/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GAACAAGGTCGCCAGGTCGACATCATCAAGCTGTATGGCTCGATGGAGCTGGCGCCGCTGATCAACCTTGCCGACAAGATCATCGACGTGGTCGACACCGGCAACACCCTGCGTGCCAACGGCCTGGAACCCCAGGAACTGATCGCCACGATCAGCTCTCGCCTGGTGGTCAACAAAGCCTCCATGAAAATGCAGCACGCCCGTATCCAGAGCCTGATCGACACGCTGCGTGCAGCG‑GTCGAATCGCGACACCGCGGTTGACTCCTGCGCGCGACCTTCGCTGTCGCGCCCGTCT  >  NC_002947/1104974‑1105268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑