Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 474,808 G→C 5.4% Y59* (TAC→TAG)  dnaG ← DNA primase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947474,8080GC 5.4% 285.8 / 3.8 87Y59* (TAC→TAG) dnaGDNA primase
Reads supporting (aligned to +/- strand):  ref base G (45/37);  new base C (2/3);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GCGGCTTCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCG  >  NC_002947/474660‑474952
                                                                                                                                                     |                                                                                                                                                
gcgGCTTCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCagtag                                                                                                                                                  >  5:37547/1‑149 (MQ=255)
    cTTCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGtagtag                                                                                                                                               >  2:80230/1‑148 (MQ=255)
     ttCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGaa                                                                                                                                             >  6:281857/1‑149 (MQ=255)
          gTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGc                                                                                                                                         <  7:318283/148‑1 (MQ=255)
            cgcgGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTTCTACTACTGCttc                                                                                                                                      <  3:227802/149‑2 (MQ=255)
             gcggcgGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTTCTACTACTGCttc                                                                                                                                      <  3:260386/148‑2 (MQ=255)
              cggcggCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTTCTACTACTGCttcgg                                                                                                                                    <  6:305621/149‑4 (MQ=255)
                 cggcCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCggg                                                                                                                                 <  4:343816/149‑1 (MQ=255)
                  ggcCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCgggg                                                                                                                                <  6:83030/149‑1 (MQ=255)
                     cTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCAGGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAGCTGGTTGTCGGGGCTg                                                                                                                             >  6:373049/1‑149 (MQ=255)
                          tCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGt                                                                                                                        >  4:162097/1‑149 (MQ=255)
                            gcgGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTg                                                                                                                       <  4:171801/148‑1 (MQ=255)
                            gcgGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGa                                                                                                                      <  2:88425/149‑1 (MQ=255)
                                    cTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGAcc                                                                                                                          >  6:357106/1‑137 (MQ=255)
                                    cTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGAcc                                                                                                                          <  5:357106/137‑1 (MQ=255)
                                     ttCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACgggg                                                                                                             >  3:262912/1‑149 (MQ=255)
                                     ttCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACgggg                                                                                                             >  8:292468/1‑149 (MQ=255)
                                         aTGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCTTGACAAAGCCCAAGGCGTTGCCACCGGCGCCAGAGCCGAATCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGtttt                                                                                                         <  3:148370/149‑1 (MQ=255)
                                         aTGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGtttt                                                                                                         >  7:2703/1‑149 (MQ=255)
                                               gcggcgCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCttt                                                                                                   <  5:281857/149‑1 (MQ=255)
                                                cggcgCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGGTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtg                                                                                                  <  7:313820/149‑1 (MQ=255)
                                                 ggcgCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtgt                                                                                                 <  2:70633/149‑1 (MQ=255)
                                                    gcgGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGCCAAAGCCCAAGGCGTTGCCACCGGCGCCACCGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTTACCGTGAAGGACGGCGTTTTTTCTTTGTgg                                                                                               >  7:276128/1‑148 (MQ=255)
                                                         ccAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGCCGGGGGTTTTTCTTTGTGGAACggg                                                                                         >  6:252831/1‑149 (MQ=255)
                                                           aGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCGTTGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACggg                                                                                         >  1:254365/1‑148 (MQ=255)
                                                             ttCCTCGACGGCCTGGGGGAAGTCCAGTTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGgcagc                                                                                     <  1:80230/149‑1 (MQ=255)
                                                                cTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGcagcag                                                                                   <  1:357046/148‑1 (MQ=255)
                                                                cTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAgg                                                                                  >  1:324513/1‑149 (MQ=255)
                                                                   gACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGcc                                                                                >  5:384401/1‑148 (MQ=255)
                                                                   gACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGcc                                                                                >  5:157755/1‑148 (MQ=255)
                                                                    aCGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGa                                                                              <  3:280531/149‑1 (MQ=255)
                                                                       gCCTGGGTGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagta                                                                           <  7:292468/149‑1 (MQ=255)
                                                                       gCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagta                                                                           <  1:323979/149‑1 (MQ=255)
                                                                        ccTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagtag                                                                          >  4:68754/1‑149 (MQ=255)
                                                                         cTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagtagt                                                                         >  6:60558/1‑149 (MQ=255)
                                                                         cTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagtagt                                                                         >  2:201257/1‑149 (MQ=255)
                                                                              ggAAGTCCAGGTTGTCG‑TGGTCCGTGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTg                                                                    >  4:62016/1‑149 (MQ=255)
                                                                               gAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGc                                                                   >  5:145680/1‑149 (MQ=255)
                                                                                   tCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGt                                                               <  2:254365/149‑1 (MQ=255)
                                                                                    ccAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGtt                                                              >  8:182791/1‑149 (MQ=255)
                                                                                     cAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtgt                                                                                                 <  1:324519/113‑1 (MQ=255)
                                                                                     cAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtgt                                                                                                 >  2:324519/1‑113 (MQ=255)
                                                                                     cAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGttt                                                             >  3:216495/1‑149 (MQ=255)
                                                                                      aGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGtttt                                                            >  7:257060/1‑149 (MQ=255)
                                                                                           gTCG‑TGGTCCCTGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTc                                                        <  3:62016/148‑1 (MQ=255)
                                                                                            tCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAg                                                      >  7:12149/1‑149 (MQ=255)
                                                                                             cG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCCGc                                                     >  5:214903/1‑149 (MQ=255)
                                                                                              g‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGc                                                     <  6:183943/148‑1 (MQ=255)
                                                                                                  gTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTgg                                                  <  2:324513/148‑1 (MQ=255)
                                                                                                  gTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGa                                                 <  6:265859/149‑1 (MQ=255)
                                                                                                  gTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCAACGCCGAAGAAGTAGTAGAAATGCTTGTGGGGGGTGTCCGTGATGGCAGGGGTTTTTTGTTTGTGTAAGGGGCAGAAGGGCGAGTAGTTCTTGCAGGTTTTTTTGAGCTGGa                                                 >  8:161180/1‑149 (MQ=255)
                                                                                                    ccATGACAATGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagt                                                                            <  5:271345/120‑1 (MQ=255)
                                                                                                    ccATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGgc                                                                                        <  4:111967/108‑1 (MQ=255)
                                                                                                    ccATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGgc                                                                                        >  3:111967/1‑108 (MQ=255)
                                                                                                    ccATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagt                                                                            >  6:271345/1‑120 (MQ=255)
                                                                                                     cATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCTAGTAGTTCTAGACGGTTTTTTTCAGCTGGAcgc                                              >  4:99460/1‑149 (MQ=255)
                                                                                                      aTGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcg                                             >  3:261014/1‑149 (MQ=255)
                                                                                                       tGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcgc                                            <  6:157755/149‑1 (MQ=255)
                                                                                                        gACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGa                                                                              <  3:345166/114‑1 (MQ=255)
                                                                                                        gACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGa                                                                              >  4:345166/1‑114 (MQ=255)
                                                                                                         aCAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcgcg                                           <  5:60558/148‑1 (MQ=255)
                                                                                                          cAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGttttttt                                                         <  6:277651/133‑1 (MQ=255)
                                                                                                          cAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGttttttt                                                         >  5:277651/1‑133 (MQ=255)
                                                                                                          cAAAGCCCAAGGCGTTGCCACCGGCGCAACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGAGAAGGACGGCGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcgcca                                          >  8:101760/1‑146 (MQ=255)
                                                                                                           aaaGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAAc                                        >  8:267358/1‑149 (MQ=255)
                                                                                                           aaaGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAACGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAAc                                        >  6:331927/1‑149 (MQ=255)
                                                                                                               cccAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTcac                                    <  8:12149/149‑1 (MQ=255)
                                                                                                                 cAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTcacca                                  >  1:218211/1‑149 (MQ=255)
                                                                                                                 cAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTcacca                                  >  1:333388/1‑149 (MQ=255)
                                                                                                                     gCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTc                              >  7:30248/1‑149 (MQ=255)
                                                                                                                       gTTGCCACCGGCGCCACATCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTcga                            <  7:331586/149‑1 (MQ=255)
                                                                                                                        ttGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTcgac                           <  6:37547/149‑1 (MQ=255)
                                                                                                                           ccACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGACCTCACCACGTCGACGAt                        >  4:273872/1‑149 (MQ=255)
                                                                                                                           ccACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGACGGTTTTTTTCAGCTGGAAGCGCGAACTCACCACGTCGACGAt                        >  1:282094/1‑149 (MQ=255)
                                                                                                                            cACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATg                       >  5:177230/1‑149 (MQ=255)
                                                                                                                             aCCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGt                      >  3:67078/1‑149 (MQ=255)
                                                                                                                                ggCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGag                                                                             <  3:185687/91‑1 (MQ=255)
                                                                                                                                ggCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGag                                                                             >  4:185687/1‑91 (MQ=255)
                                                                                                                                 gcgcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAg                  <  4:262912/149‑1 (MQ=255)
                                                                                                                                 gcgcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAg                  >  6:274253/1‑149 (MQ=255)
                                                                                                                                 gcgcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAg                  >  5:77424/1‑149 (MQ=255)
                                                                                                                                   gcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGc                <  8:257060/149‑1 (MQ=255)
                                                                                                                                    ccACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCATGTCGACGATGTCGAGGCg               <  1:201257/149‑1 (MQ=255)
                                                                                                                                         cgacaaGCAGTTCTACTACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTGGa          >  2:312218/5‑149 (MQ=255)
                                                                                                                                         cgacaaGCAGTTCTACTACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTATCTTTGTGGAACGCGCAGCAGGCCGAGTACTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGCCGATGTCGAGGCGGTTgc          >  5:283401/5‑148 (MQ=255)
                                                                                                                                           cGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATc        <  7:267358/149‑1 (MQ=255)
                                                                                                                                             aaGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAg      <  3:68754/149‑1 (MQ=255)
                                                                                                                                              aGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTCCGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAgg     <  4:76930/149‑1 (MQ=255)
                                                                                                                                                cAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGgtc   >  2:55667/1‑149 (MQ=255)
                                                                                                                                                 agtagtAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGgtcg  <  4:67078/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCGGCTTCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCG  >  NC_002947/474660‑474952

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