Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,436,200 G→A 5.4% L45L (CTG→CTA) ‡ gpI → tail formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,436,2000GA 5.4% 288.9 / 3.7 92L45L (CTG→CTA) ‡gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base G (42/45);  new base A (2/3);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGG  >  NC_002947/3436056‑3436338
                                                                                                                                                  |                                                                                                                                            
aGGTGGTGACATGAGCCTGTTACCGCACAACTCCTCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTTGCTGGCCt                                                                                                                                          <  7:117617/149‑1 (MQ=255)
  gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                        <  4:311260/149‑1 (MQ=255)
  gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                        <  8:306824/149‑1 (MQ=255)
   tggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTggg                                                                                                                                       <  1:280830/149‑1 (MQ=255)
   tggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCAGGCCAGCCAggg                                                                                                                                       <  1:208816/149‑1 (MQ=255)
     gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCa                                                                                                                                     >  6:313906/1‑149 (MQ=255)
          aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑gg                                                                                                                                 >  5:321721/1‑147 (MQ=255)
          aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGaa                                                                                                                               <  4:110262/149‑1 (MQ=255)
            gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAgg                                                                                                                             >  4:321560/1‑149 (MQ=255)
            gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAgg                                                                                                                             >  6:239112/1‑149 (MQ=255)
              gCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATAATCCGGGGCGTGGCCGACTCGTCACGC‑GGCCGGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTc                                                                                                                           <  2:340323/149‑1 (MQ=255)
                   ttACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAgg                                                                                                                       <  1:325420/148‑1 (MQ=255)
                      ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCt                                                                                                                    <  5:56986/148‑1 (MQ=255)
                       cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTgg                                                                                                                  <  5:313906/149‑1 (MQ=255)
                        gCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTggg                                                                                                                 <  2:292429/149‑1 (MQ=255)
                        gCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTgg                                                                                                                  >  1:292429/1‑148 (MQ=255)
                           caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAgg                                                                                                              >  4:115716/1‑149 (MQ=255)
                           caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAgg                                                                                                              >  8:304838/1‑149 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGaa                                                                                                                               <  4:47522/123‑1 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                     <  5:239112/149‑1 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                     >  3:164309/1‑149 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                     >  7:261286/1‑149 (MQ=255)
                                     cAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTacac                                                                                                    <  6:321721/149‑1 (MQ=255)
                                        ttGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACAGTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGATCGTCGAGGGCTGGGAGGCGGCTTACACcga                                                                                                 >  5:317113/1‑149 (MQ=255)
                                                 gcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTc                                                                                                                           <  6:201618/114‑1 (MQ=255)
                                                 gcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTc                                                                                                                           >  5:201618/1‑114 (MQ=255)
                                                  tccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAgc                                                                                         <  3:60995/146‑1 (MQ=255)
                                                  cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAgg                                                                                                                             >  8:279524/1‑111 (MQ=255)
                                                  cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAgg                                                                                                                             <  7:279524/111‑1 (MQ=255)
                                                    cTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACCCCGACGACCAgc                                                                                         >  4:60995/1‑145 (MQ=255)
                                                     tGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgccg‑c                                                                                   >  2:319493/1‑149 (MQ=255)
                                                        aGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcg‑cgcg                                                                                <  8:206884/149‑1 (MQ=255)
                                                          gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgac                                                                                                >  8:213249/1‑132 (MQ=255)
                                                          gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgac                                                                                                <  7:213249/132‑1 (MQ=255)
                                                                  cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCcgc                                                                       >  8:252962/1‑148 (MQ=255)
                                                                   gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGa                                                                     <  7:304838/149‑1 (MQ=255)
                                                                   gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGa                                                                     <  8:261286/149‑1 (MQ=255)
                                                                       ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGCCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGCCGGCTTACACCGACTACCAGCGCCG‑CGCGCTGATCCGCTAGGc                                                                  >  5:207316/1‑148 (MQ=255)
                                                                          gccacgccccGGATAATGTCCGGGTCAAGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCAtt                                                              >  6:112842/8‑149 (MQ=255)
                                                                           ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTc                                                             <  5:87358/149‑1 (MQ=255)
                                                                            cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTcc                                                            >  8:9681/1‑149 (MQ=255)
                                                                             ttGACCCGGACATTATCCGGGGCGTGGCGGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCg                                                           >  1:197314/1‑149 (MQ=255)
                                                                             ttGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCg                                                           >  3:262793/1‑149 (MQ=255)
                                                                               gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCgg                                                          <  2:197314/148‑1 (MQ=255)
                                                                                aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGt                                                         <  6:65868/148‑1 (MQ=255)
                                                                                 cccGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGtt                                                        >  4:117850/1‑148 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCGTGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGtt                                                        >  1:233464/1‑148 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCg                                                           <  4:342033/144‑1 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCg                                                           >  3:342033/1‑144 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTca                                                      >  1:313209/1‑149 (MQ=255)
                                                                                        aTTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGATGCGGCTTACACCGACGACCAGCGCCGCCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGa                                                 <  1:223893/149‑1 (MQ=255)
                                                                                        aTTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                >  4:147945/1‑149 (MQ=255)
                                                                                        aTTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                <  3:117850/149‑1 (MQ=255)
                                                                                            tCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAg                                            <  4:262793/149‑1 (MQ=255)
                                                                                              cGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                           <  3:321560/148‑1 (MQ=255)
                                                                                                  gCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGa                                                 >  3:73276/1‑138 (MQ=255)
                                                                                                  gCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGa                                                 <  4:73276/138‑1 (MQ=255)
                                                                                                    gTGGCCGACTCGTCACGC‑TGCCCGCCGAACTACC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGt                                    >  7:318710/1‑149 (MQ=255)
                                                                                                       gCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTgg                                                                                                                  >  6:295999/1‑69 (MQ=255)
                                                                                                       gCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTgg                                                                                                                  <  5:295999/69‑1 (MQ=255)
                                                                                                        ccGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgc                                <  7:252962/149‑1 (MQ=255)
                                                                                                            cTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGATGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                            <  3:115716/149‑1 (MQ=255)
                                                                                                               gTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCaa                                             <  8:233834/129‑1 (MQ=255)
                                                                                                               gTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCaa                                             >  7:233834/1‑129 (MQ=255)
                                                                                                                tCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCa                                        >  1:242452/1‑134 (MQ=255)
                                                                                                                tCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCa                                        <  2:242452/133‑1 (MQ=255)
                                                                                                                 cACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGAt                                                                           <  6:264219/97‑1 (MQ=255)
                                                                                                                    gC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGTCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                    <  1:77313/149‑1 (MQ=255)
                                                                                                                    gC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                    <  3:147945/149‑1 (MQ=255)
                                                                                                                    gC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                    >  2:276336/1‑149 (MQ=255)
                                                                                                                       tGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGtt                                                        <  6:209928/111‑1 (MQ=255)
                                                                                                                       tGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGtt                                                        >  5:209928/1‑111 (MQ=255)
                                                                                                                          ccgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGCGCACGGTCGGCGCGATTCGGCGGGTGc                >  2:332337/1‑148 (MQ=255)
                                                                                                                           cgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGc                                                                  <  3:157111/97‑1 (MQ=255)
                                                                                                                           cgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGc                                                                  >  4:157111/1‑97 (MQ=255)
                                                                                                                           cgccgAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc              <  1:319493/149‑1 (MQ=255)
                                                                                                                             tcgACCTTCCATGCCCAGGCCAGCCAGGGCAT‑‑GAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTc                                                       >  7:120122/2‑106 (MQ=21)
                                                                                                                             tcgACCTTCCATGCCCAGGCCAGCCAGGGCAT‑‑GAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTc                                                       <  8:120122/105‑1 (MQ=21)
                                                                                                                             ccgAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc              >  4:278080/1‑147 (MQ=255)
                                                                                                                              cgACCTTCCATGCCCAGGCCAGCCAGGGCA‑‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                            <  8:118937/132‑1 (MQ=37)
                                                                                                                              cgACCTTCCATGCCCAGGCCAGCCAGGGCA‑‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                            >  7:118937/1‑132 (MQ=37)
                                                                                                                              cgAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAg           <  1:276336/149‑1 (MQ=255)
                                                                                                                               gAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg          <  8:318710/149‑1 (MQ=255)
                                                                                                                                 aCTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCg        <  5:33727/149‑1 (MQ=255)
                                                                                                                                  cTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCCGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCGCAcggcgg       >  6:26275/1‑149 (MQ=255)
                                                                                                                                  cTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCTCAAGGCgg       >  6:21046/1‑149 (MQ=255)
                                                                                                                                     cc‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcg‑cgcg                                                                                <  2:223279/73‑1 (MQ=255)
                                                                                                                                     cc‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcg‑cgcg                                                                                >  1:223279/1‑73 (MQ=255)
                                                                                                                                      c‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg   >  7:179699/1‑149 (MQ=255)
                                                                                                                                        tGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg  >  4:177577/1‑149 (MQ=255)
                                                                                                                                          cccTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg   >  3:212791/1‑146 (MQ=255)
                                                                                                                                          cccTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg   <  4:212791/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
AGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGG  >  NC_002947/3436056‑3436338

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: ‑