Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,285 C→G 5.4% V103V (GTG→GTC) ‡ trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,2850CG 5.4% 614.0 / 18.4 184V103V (GTG→GTC) ‡trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (91/83);  new base G (2/8);  total (93/91)
Fisher's exact test for biased strand distribution p-value = 5.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 191 total aligned reads displayed.

CGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGAAA  >  NC_002947/93138‑93431
                                                                                                                                                   |                                                                                                                                                  
cGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCAc                                                                                                                                                  <  1:4184727/150‑1 (MQ=255)
      ctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGcc                                                                                                                                            <  1:3383445/150‑1 (MQ=255)
       tgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCg                                                                                                                                           >  1:548976/1‑150 (MQ=255)
         gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                                                                                         >  2:1079777/1‑150 (MQ=255)
         gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                                                                                         >  2:4247756/1‑150 (MQ=255)
          cgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATg                                                                                                                                        <  2:4074267/150‑1 (MQ=255)
            cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                                                                      >  2:5218795/1‑150 (MQ=255)
            cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                                                                      >  2:3913644/1‑150 (MQ=255)
             tCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAg                                                                                                                                     <  1:4645392/150‑1 (MQ=255)
               aacgTCAGTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAatt                                                                                                                                   <  1:1296038/147‑1 (MQ=255)
                   tcagTGGCACCCATCGAGACCCCGCCAGGCAGGCCGAACCGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgtta                                                                                                                               <  2:3024270/147‑2 (MQ=255)
                   tcagTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgtta                                                                                                                               <  2:3024314/147‑2 (MQ=255)
                    cGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCa                                                                                                                              >  1:2055586/1‑150 (MQ=255)
                      gTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAcc                                                                                                                            <  2:1686963/150‑1 (MQ=255)
                      gTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCAcc                                                                                                                            >  1:5101163/1‑150 (MQ=255)
                      gTGGCACCCATATAGACCACGCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAattgttaatc                                                                                                                            >  1:3459595/1‑146 (MQ=255)
                       tGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCt                                                                                                                           >  1:1312396/1‑150 (MQ=255)
                         gCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTtg                                                                                                                         >  2:3549636/1‑150 (MQ=255)
                           aCCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTtgtg                                                                                                                       >  2:4890271/1‑150 (MQ=255)
                            cccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgc                                                                                                                      <  2:1126703/150‑1 (MQ=255)
                            cccATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgc                                                                                                                      <  1:1712426/150‑1 (MQ=255)
                              cATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgc                                                                                                                    >  1:2767397/1‑150 (MQ=255)
                               aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                                                                                                   >  2:1735634/1‑150 (MQ=255)
                               aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                                                                                                   >  2:4769517/1‑150 (MQ=255)
                                  tAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCg                                                                                                                <  2:4448358/150‑1 (MQ=255)
                                      tCACGCAAGGCAGGCCGAAGCGCGCGGTGACAGTGGAGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCGTGTGCGCGCCGGTAt                                                                                                            >  2:4875624/1‑150 (MQ=255)
                                       cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                                                                                                  >  1:3342989/1‑127 (MQ=255)
                                       cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                                                                                                  <  2:3342904/127‑1 (MQ=255)
                                        aCGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATgg                                                                                                          >  1:270534/1‑150 (MQ=255)
                                          gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                                                                        <  1:3311439/150‑1 (MQ=255)
                                          gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                                                                        >  1:4276336/1‑150 (MQ=255)
                                          gCAAGGCAGGCCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATGGtt                                                                                                        >  1:2948899/1‑150 (MQ=255)
                                            aaGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGa                                                                                                      >  2:1520583/1‑150 (MQ=255)
                                               gcaggcCGAAACGGGCTGCAACGGTTGCTGATGCAACGCCATGCTGACCAGCGCCGGTTTCCGCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATGGTTCAggt                                                                                                   <  1:5000131/150‑3 (MQ=255)
                                                caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                                                                                                  >  2:3167655/1‑150 (MQ=255)
                                                caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                                                                                                  <  2:3480785/150‑1 (MQ=255)
                                                  ggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCtt                                                                                                <  1:1202865/150‑1 (MQ=255)
                                                   gcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTc                                                                                               >  1:111192/1‑150 (MQ=255)
                                                    ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCa                                                                                              <  2:668709/150‑1 (MQ=255)
                                                    ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCa                                                                                              <  2:3602752/150‑1 (MQ=255)
                                                    ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCa                                                                                              <  2:3901461/150‑1 (MQ=255)
                                                     cGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCAc                                                                                             >  2:1629273/1‑150 (MQ=255)
                                                         gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                                                                         <  2:2672009/150‑1 (MQ=255)
                                                           gcgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTga                                                                                       <  1:3313713/150‑1 (MQ=255)
                                                             gcggcgACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaag                                                                                     >  2:2189276/1‑150 (MQ=255)
                                                               ggcgACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCTGTATGGTTGAGCTCTTCACGCTTgaagaa                                                                                   <  1:1165575/150‑1 (MQ=255)
                                                                     aGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCtt                                                                             >  1:4525777/1‑150 (MQ=255)
                                                                      gtggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTc                                                                            <  1:3957467/150‑1 (MQ=255)
                                                                       tggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCg                                                                           >  2:5346256/1‑150 (MQ=255)
                                                                         gcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCt                                                                         >  2:4221908/1‑150 (MQ=255)
                                                                         gcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCt                                                                         >  1:594735/1‑150 (MQ=255)
                                                                              ggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTccgcc                                                                    >  2:3855112/1‑150 (MQ=255)
                                                                                 gACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCa                                                                 >  2:1932201/1‑150 (MQ=255)
                                                                                        tGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCg                                                          <  1:489149/150‑1 (MQ=255)
                                                                                           aTGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTc                                                       >  1:2607933/1‑150 (MQ=255)
                                                                                             gCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAg                                                     <  1:1683629/150‑1 (MQ=255)
                                                                                              ccGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAgg                                                    <  1:1652206/150‑1 (MQ=255)
                                                                                              ccGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAgg                                                    <  2:3303982/150‑1 (MQ=255)
                                                                                                 gcgcCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGgcgc                                                 >  1:170855/1‑150 (MQ=255)
                                                                                                  cgcCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCt                                                <  2:64405/150‑1 (MQ=255)
                                                                                                    ccGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCg                                              <  1:3679574/150‑1 (MQ=255)
                                                                                                     cGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCgg                                             <  1:5307344/150‑1 (MQ=255)
                                                                                                     cGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCgg                                             >  1:3494516/1‑150 (MQ=255)
                                                                                                      ggTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGTCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGc                                            <  1:70004/150‑1 (MQ=255)
                                                                                                      ggTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGc                                            <  1:4461640/150‑1 (MQ=255)
                                                                                                        tttCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGa                                          >  2:1134590/1‑150 (MQ=255)
                                                                                                         ttCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGaa                                         <  2:3825468/150‑1 (MQ=255)
                                                                                                            ggCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTa                                      >  2:2209065/1‑150 (MQ=255)
                                                                                                            cgCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATGGTTCAGGTCCTCGCGCTTAAAGTAAATCTTTGc                                                                          <  2:4057982/113‑1 (MQ=255)
                                                                                                             gCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTAc                                     >  1:1043038/1‑150 (MQ=255)
                                                                                                                aTCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGc                                  <  1:198458/150‑1 (MQ=255)
                                                                                                                  cAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCgg                                <  1:1778446/150‑1 (MQ=255)
                                                                                                                     gcgcTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                                                                         <  2:2527087/90‑1 (MQ=255)
                                                                                                                         tttttGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTgggg                         >  1:1374894/1‑150 (MQ=255)
                                                                                                                            ttGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGggcgg                      <  1:3855112/150‑1 (MQ=255)
                                                                                                                             tGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGggcggc                     <  2:4285219/150‑1 (MQ=255)
                                                                                                                                ccATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGa                  >  2:699985/1‑150 (MQ=255)
                                                                                                                                ccATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGa                  <  1:1332140/150‑1 (MQ=255)
                                                                                                                                ccATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGa                  <  2:4170010/150‑1 (MQ=255)
                                                                                                                                 cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                 >  1:4391909/1‑150 (MQ=255)
                                                                                                                                 cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                 <  1:4012083/150‑1 (MQ=255)
                                                                                                                                 cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGCTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                 <  1:2912915/150‑1 (MQ=255)
                                                                                                                                  aTGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATg                <  2:805464/150‑1 (MQ=255)
                                                                                                                                   tGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGt               >  2:1259274/1‑150 (MQ=255)
                                                                                                                                   tGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGt               >  2:2846915/1‑150 (MQ=255)
                                                                                                                                    gcgcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTa              <  1:2558994/150‑1 (MQ=255)
                                                                                                                                      gcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGt            >  1:2517279/1‑150 (MQ=255)
                                                                                                                                      gcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGt            >  2:2392869/1‑150 (MQ=255)
                                                                                                                                      gcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGt            >  2:5298417/1‑150 (MQ=255)
                                                                                                                                       cTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTc           <  2:3660859/150‑1 (MQ=255)
                                                                                                                                        ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcg          >  2:321712/1‑150 (MQ=255)
                                                                                                                                        ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcg          <  1:2209065/150‑1 (MQ=255)
                                                                                                                                        ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAGGTAGTcg          >  2:4878797/1‑150 (MQ=255)
                                                                                                                                        ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGCTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcg          <  2:1372309/150‑1 (MQ=255)
                                                                                                                                         tGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcgc         <  2:548976/150‑1 (MQ=255)
                                                                                                                                          ggCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcgcg        >  2:1241328/1‑150 (MQ=255)
                                                                                                                                                cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa  <  2:2571897/150‑1 (MQ=255)
                                                                                                                                                cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa  <  2:2354733/150‑1 (MQ=255)
                                                                                                                                                cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa  <  1:1259274/150‑1 (MQ=255)
                                                                                                                                                cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa  <  2:960318/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCCTGCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGAAA  >  NC_002947/93138‑93431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: