Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,436,195 T→A 5.4% W44R (TGG→AGG)  gpI → tail formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,436,1950TA 5.4% 281.8 / 4.6 92W44R (TGG→AGG) gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base T (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 = 7.26e-01

AAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGG  >  NC_002947/3436055‑3436338
                                                                                                                                              |                                                                                                                                                 
aaGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCcaggccagc                                                                                                                                            >  2:208816/1‑144 (MQ=255)
 aGGTGGTGACATGAGCCTGTTACCGCACAACTCCTCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTTGCTGGCCt                                                                                                                                          <  7:117617/149‑1 (MQ=255)
   gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                        <  8:306824/149‑1 (MQ=255)
   gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                        <  4:311260/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                                                                                                     >  7:261286/1‑149 (MQ=255)
                                     ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                     <  5:239112/149‑1 (MQ=255)
                                     ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                     >  3:164309/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                                                                                                                             <  7:279524/111‑1 (MQ=255)
                                                   cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAgg                                                                                                                             >  8:279524/1‑111 (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                                                                                                <  7:213249/132‑1 (MQ=255)
                                                           gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgac                                                                                                >  8:213249/1‑132 (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                                                <  3:117850/149‑1 (MQ=255)
                                                                                         aTTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                >  4:147945/1‑149 (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                                                                                                                  <  5:295999/69‑1 (MQ=255)
                                                                                                        gCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTgg                                                                                                                  >  6:295999/1‑69 (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                                             >  7:233834/1‑129 (MQ=255)
                                                                                                                gTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCaa                                             <  8:233834/129‑1 (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                            >  7:118937/1‑132 (MQ=37)
                                                                                                                               cgACCTTCCATGCCCAGGCCAGCCAGGGCA‑‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                            <  8:118937/132‑1 (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)
                                                                                                                                              |                                                                                                                                                 
AAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGC‑TGCCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCAT‑GGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCG‑CGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGG  >  NC_002947/3436055‑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: ‑