Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,436,204 T→A 5.3% W47R (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,2040TA 5.3% 290.3 / 2.8 93W47R (TGG→AGG) gpItail formation protein
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base A (2/3);  total (45/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 = 1.00e+00

AGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATT  >  NC_002947/3436056‑3436351
                                                                                                                                                      |                                                                                                                                                     
aGGTGGTGACATGAGCCTGTTACCGCACAACTCCTCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTTGCTGGCCt                                                                                                                                                       <  7:117617/149‑1 (MQ=255)
  gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                     <  8:306824/149‑1 (MQ=255)
  gtggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTgg                                                                                                                                                     <  4:311260/149‑1 (MQ=255)
   tggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTggg                                                                                                                                                    <  1:280830/149‑1 (MQ=255)
   tggtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCAGGCCAGCCAggg                                                                                                                                                    <  1:208816/149‑1 (MQ=255)
     gtgACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCa                                                                                                                                                  >  6:313906/1‑149 (MQ=255)
          aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATgg                                                                                                                                              >  5:321721/1‑147 (MQ=255)
          aTGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGaa                                                                                                                                            <  4:110262/149‑1 (MQ=255)
            gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAgg                                                                                                                                          >  4:321560/1‑149 (MQ=255)
            gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAgg                                                                                                                                          >  6:239112/1‑149 (MQ=255)
              gCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATAATCCGGGGCGTGGCCGACTCGTCACGCGGCCGGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTc                                                                                                                                        <  2:340323/149‑1 (MQ=255)
                   ttACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAgg                                                                                                                                    <  1:325420/148‑1 (MQ=255)
                      ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                                 <  5:56986/148‑1 (MQ=255)
                       cGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                               <  5:313906/149‑1 (MQ=255)
                        gCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTggg                                                                                                                              <  2:292429/149‑1 (MQ=255)
                        gCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                               >  1:292429/1‑148 (MQ=255)
                           caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAgg                                                                                                                           >  8:304838/1‑149 (MQ=255)
                           caACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAgg                                                                                                                           >  4:115716/1‑149 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGaa                                                                                                                                            <  4:47522/123‑1 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                                  <  5:239112/149‑1 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                                  >  3:164309/1‑149 (MQ=255)
                                    ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTaca                                                                                                                  >  7:261286/1‑149 (MQ=255)
                                     cAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTacac                                                                                                                 <  6:321721/149‑1 (MQ=255)
                                        ttGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACAGTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGATCGTCGAGGGCTGGGAGGCGGCTTACACcga                                                                                                              >  5:317113/1‑149 (MQ=255)
                                                 gcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTc                                                                                                                                        <  6:201618/114‑1 (MQ=255)
                                                 gcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTc                                                                                                                                        >  5:201618/1‑114 (MQ=255)
                                                  tccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAgc                                                                                                      <  3:60995/146‑1 (MQ=255)
                                                  cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAgg                                                                                                                                          >  8:279524/1‑111 (MQ=255)
                                                  cccTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAgg                                                                                                                                          <  7:279524/111‑1 (MQ=255)
                                                    cTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACCCCGACGACCAgc                                                                                                      >  4:60995/1‑145 (MQ=255)
                                                     tGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGcgccgc                                                                                                >  2:319493/1‑149 (MQ=255)
                                                        aGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcgcg                                                                                             <  8:206884/149‑1 (MQ=255)
                                                          gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgac                                                                                                             >  8:213249/1‑132 (MQ=255)
                                                          gcggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcgac                                                                                                             <  7:213249/132‑1 (MQ=255)
                                                                  cgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                                    >  8:252962/1‑148 (MQ=255)
                                                                   gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                                  <  8:261286/149‑1 (MQ=255)
                                                                   gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGa                                                                                  <  7:304838/149‑1 (MQ=255)
                                                                       ttGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGCCCTGGGCATGGAAGGTCGAGGGCTGGGAGCCGGCTTACACCGACTACCAGCGCCGCGCGCTGATCCGCTAGGc                                                                               >  5:207316/1‑148 (MQ=255)
                                                                          gccacgccccGGATAATGTCCGGGTCAAGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCAtt                                                                           >  6:112842/8‑149 (MQ=255)
                                                                           ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTc                                                                          <  5:87358/149‑1 (MQ=255)
                                                                            cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTcc                                                                         >  8:9681/1‑149 (MQ=255)
                                                                             ttGACCCGGACATTATCCGGGGCGTGGCGGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                                        >  1:197314/1‑149 (MQ=255)
                                                                             ttGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                                        >  3:262793/1‑149 (MQ=255)
                                                                               gACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCgg                                                                       <  2:197314/148‑1 (MQ=255)
                                                                                aCCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGt                                                                      <  6:65868/148‑1 (MQ=255)
                                                                                 cccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                     >  4:117850/1‑148 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCGTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                     >  1:233464/1‑148 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                                        >  3:342033/1‑144 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCg                                                                        <  4:342033/144‑1 (MQ=255)
                                                                                  ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                                   >  1:313209/1‑149 (MQ=255)
                                                                                        aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGATGCGGCTTACACCGACGACCAGCGCCGCCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGa                                                              <  1:223893/149‑1 (MQ=255)
                                                                                        aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                             >  4:147945/1‑149 (MQ=255)
                                                                                        aTTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGAc                                                             <  3:117850/149‑1 (MQ=255)
                                                                                            tCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAg                                                         <  4:262793/149‑1 (MQ=255)
                                                                                              cGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                                        <  3:321560/148‑1 (MQ=255)
                                                                                                  gCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGa                                                              <  4:73276/138‑1 (MQ=255)
                                                                                                  gCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGa                                                              >  3:73276/1‑138 (MQ=255)
                                                                                                    gTGGCCGACTCGTCACGCTGCCCGCCGAACTACCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGt                                                 >  7:318710/1‑149 (MQ=255)
                                                                                                       gCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                               <  5:295999/69‑1 (MQ=255)
                                                                                                       gCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTgg                                                                                                                               >  6:295999/1‑69 (MQ=255)
                                                                                                        ccGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgc                                             <  7:252962/149‑1 (MQ=255)
                                                                                                            cTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGATGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                         <  3:115716/149‑1 (MQ=255)
                                                                                                               gTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCaa                                                          <  8:233834/129‑1 (MQ=255)
                                                                                                               gTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCaa                                                          >  7:233834/1‑129 (MQ=255)
                                                                                                                tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCa                                                     >  1:242452/1‑134 (MQ=255)
                                                                                                                tCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCa                                                     <  2:242452/133‑1 (MQ=255)
                                                                                                                 cACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGAt                                                                                        <  6:264219/97‑1 (MQ=255)
                                                                                                                    gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGTCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                 <  1:77313/149‑1 (MQ=255)
                                                                                                                    gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                 >  2:276336/1‑149 (MQ=255)
                                                                                                                    gCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                                 <  3:147945/149‑1 (MQ=255)
                                                                                                                       tGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                     >  5:209928/1‑111 (MQ=255)
                                                                                                                       tGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                                     <  6:209928/111‑1 (MQ=255)
                                                                                                                          ccgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGCGCACGGTCGGCGCGATTCGGCGGGTGc                             >  2:332337/1‑148 (MQ=255)
                                                                                                                           cgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGc                                                                               <  3:157111/97‑1 (MQ=255)
                                                                                                                           cgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGc                                                                               >  4:157111/1‑97 (MQ=255)
                                                                                                                           cgccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                           <  1:319493/149‑1 (MQ=255)
                                                                                                                             tcgACCTTCCATGCCCAGGCCAGCCAGGGCATGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTc                                                                    >  7:120122/2‑106 (MQ=21)
                                                                                                                             tcgACCTTCCATGCCCAGGCCAGCCAGGGCATGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTc                                                                    <  8:120122/105‑1 (MQ=21)
                                                                                                                             ccgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTc                           >  4:278080/1‑147 (MQ=255)
                                                                                                                              cgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                         >  7:118937/1‑132 (MQ=37)
                                                                                                                              cgACCTTCCATGCCCAGGCCAGCCAGGGCAGGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGa                                         <  8:118937/132‑1 (MQ=37)
                                                                                                                              cgAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAg                        <  1:276336/149‑1 (MQ=255)
                                                                                                                               gAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg                       <  8:318710/149‑1 (MQ=255)
                                                                                                                                 aCTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCg                     <  5:33727/149‑1 (MQ=255)
                                                                                                                                  cTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCCGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCGCAcggcgg                    >  6:26275/1‑149 (MQ=255)
                                                                                                                                  cTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCAAGGTCGGCGCGATTCGGCGGGTGCTCAAGGCgg                    >  6:21046/1‑149 (MQ=255)
                                                                                                                                     ccTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcgcg                                                                                             <  2:223279/73‑1 (MQ=255)
                                                                                                                                     ccTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCcgcgcg                                                                                             >  1:223279/1‑73 (MQ=255)
                                                                                                                                      cTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                >  7:179699/1‑149 (MQ=255)
                                                                                                                                        tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg               >  4:177577/1‑149 (MQ=255)
                                                                                                                                          cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                >  3:212791/1‑146 (MQ=255)
                                                                                                                                          cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCg                <  4:212791/146‑1 (MQ=255)
                                                                                                                                                     cTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGAtt  <  2:313209/149‑1 (MQ=255)
                                                                                                                                                     cTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGAt   >  2:29842/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGCCCGCCGAACTTCCTGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGGGTCAACGCGGATT  >  NC_002947/3436056‑3436351

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: