Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 93,288 T→C 5.1% Q102Q (CAA→CAG trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294793,2880TC 5.1% 625.2 / 12.1 193Q102Q (CAA→CAGtrpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (94/89);  new base C (2/8);  total (96/97)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

GCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGAAAGTAAGC  >  NC_002947/93143‑93437
                                                                                                                                                 |                                                                                                                                                     
gctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGc                                                                                                                                                   <  1:3986700/150‑1 (MQ=255)
 ctgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGcc                                                                                                                                                  <  1:3383445/150‑1 (MQ=255)
  tgGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCg                                                                                                                                                 >  1:548976/1‑150 (MQ=255)
    gcgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGAt                                                                                                                                               >  2:4247756/1‑150 (MQ=255)
     cgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATg                                                                                                                                              <  1:601413/150‑1 (MQ=255)
     cgcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATg                                                                                                                                              <  2:4074267/150‑1 (MQ=255)
      gcTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGc                                                                                                                                             >  2:4829034/1‑150 (MQ=255)
       cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                                                                            <  1:2427726/150‑1 (MQ=255)
       cTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCa                                                                                                                                            >  2:5218795/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)
                 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)
                         cATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgc                                                                                                                          >  1:2767397/1‑150 (MQ=255)
                          aTGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTgcgcg                                                                                                                         >  2:1735634/1‑150 (MQ=255)
                             tAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCg                                                                                                                      <  2:4448358/150‑1 (MQ=255)
                                  cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                                                                                                        <  2:3342904/127‑1 (MQ=255)
                                  cACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAgttg                                                                                                                                        >  1:3342989/1‑127 (MQ=255)
                                   aCGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATgg                                                                                                                >  1:270534/1‑150 (MQ=255)
                                     gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                                                                              <  1:3311439/150‑1 (MQ=255)
                                     gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                                                                              >  2:1254068/1‑150 (MQ=255)
                                     gCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGtt                                                                                                              >  1:4276336/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)
                                           caggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGctc                                                                                                        >  2:1522086/1‑150 (MQ=255)
                                             ggcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCtt                                                                                                      <  1:1202865/150‑1 (MQ=255)
                                              gcCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTc                                                                                                     >  1:111192/1‑150 (MQ=255)
                                               ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCa                                                                                                    <  2:3602752/150‑1 (MQ=255)
                                               ccGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCa                                                                                                    <  2:668709/150‑1 (MQ=255)
                                                    gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                                                                               <  2:2672009/150‑1 (MQ=255)
                                                    gcgcgcGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCtt                                                                                               >  2:1655942/1‑150 (MQ=255)
                                                      gcgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTga                                                                                             <  1:3313713/150‑1 (MQ=255)
                                                       cgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaa                                                                                            <  2:3906596/150‑1 (MQ=255)
                                                       cgcgGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTgaa                                                                                            <  2:1020967/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:3957457/150‑1 (MQ=255)
                                                                  tggcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCg                                                                                 >  2:5346256/1‑150 (MQ=255)
                                                                    gcggtggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCt                                                                               >  2:4221908/1‑150 (MQ=255)
                                                                         ggcgACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTccgcc                                                                          >  2:3855112/1‑150 (MQ=255)
                                                                            gACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCa                                                                       >  2:1932201/1‑150 (MQ=255)
                                                                             acacCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCaa                                                                      <  2:113248/150‑1 (MQ=255)
                                                                                   tGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCg                                                                <  1:489149/150‑1 (MQ=255)
                                                                                      aTGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTc                                                             >  1:2607933/1‑150 (MQ=255)
                                                                                        gCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAg                                                           <  1:1683629/150‑1 (MQ=255)
                                                                                         ccGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAgg                                                          <  2:3303982/150‑1 (MQ=255)
                                                                                         ccGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAgg                                                          <  1:1652206/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:2027258/1‑150 (MQ=255)
                                                                                                cGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCgg                                                   <  1:5307344/150‑1 (MQ=255)
                                                                                                cGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCgg                                                   >  1:3494516/1‑150 (MQ=255)
                                                                                                cGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCgg                                                   <  2:1196605/150‑1 (MQ=255)
                                                                                                 ggTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGc                                                  <  1:4461640/150‑1 (MQ=255)
                                                                                                    ttCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTGGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGaa                                               <  1:1747238/150‑1 (MQ=255)
                                                                                                    ttCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGaa                                               <  2:3825468/150‑1 (MQ=255)
                                                                                                       ggCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTa                                            >  2:2209065/1‑150 (MQ=255)
                                                                                                       ggCGATCAGGCGCTTGTTGACAATGCGCTTGGCCAGCAGCACTTGGCCGATGAAGGTGTTCACCTTGTGCTCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGAGCTCGGCGAAGTa                                            >  2:109252/1‑150 (MQ=255)
                                                                                                       cgCGATGATGCGGGTTTTGCCCATGCGCTTGGCAAGCAGGATCTGGCCGAGGCAATTGTTAATCTTGTGGCTGCCGGTATGGTTCAGGTCCTCGCGCTTAAAGTAAATCTTTGc                                                                                <  2:4057982/113‑1 (MQ=255)
                                                                                                        gCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTAc                                           >  1:1043038/1‑150 (MQ=255)
                                                                                                             cAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCgg                                      <  1:1778483/150‑1 (MQ=255)
                                                                                                              aGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCggg                                     <  2:4255684/150‑1 (MQ=255)
                                                                                                              aGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCGTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCggg                                     <  2:1446453/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)
                                                                                                                       ttGCCCATGCGCTTGGCCAGCAGCACTGGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGggcgg                            <  1:298255/150‑1 (MQ=255)
                                                                                                                           ccATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGa                        <  2:4170010/150‑1 (MQ=255)
                                                                                                                           ccATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGa                        >  2:699985/1‑150 (MQ=255)
                                                                                                                           ccATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGa                        <  1:1332140/150‑1 (MQ=255)
                                                                                                                            cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                       <  1:4012083/150‑1 (MQ=255)
                                                                                                                            cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                       <  1:2192097/150‑1 (MQ=255)
                                                                                                                            cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                       >  1:4391909/1‑150 (MQ=255)
                                                                                                                            cATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGCTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAt                       <  1:2912915/150‑1 (MQ=255)
                                                                                                                             aTGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATg                      <  2:805464/150‑1 (MQ=255)
                                                                                                                              tGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGt                     >  2:2846915/1‑150 (MQ=255)
                                                                                                                              tGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGt                     >  2:1259274/1‑150 (MQ=255)
                                                                                                                               gcgcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTa                    <  1:2558994/150‑1 (MQ=255)
                                                                                                                                 gcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGt                  >  2:2392869/1‑150 (MQ=255)
                                                                                                                                 gcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGt                  >  2:5298417/1‑150 (MQ=255)
                                                                                                                                 gcTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGt                  >  1:2517279/1‑150 (MQ=255)
                                                                                                                                  cTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTc                 <  2:3660859/150‑1 (MQ=255)
                                                                                                                                   ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcg                >  2:321712/1‑150 (MQ=255)
                                                                                                                                   ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGAGGTAGTcg                >  2:4878797/1‑150 (MQ=255)
                                                                                                                                   ttGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGCTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcg                <  2:1372309/150‑1 (MQ=255)
                                                                                                                                    tGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTcgc               <  2:548976/150‑1 (MQ=255)
                                                                                                                                           cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa        <  1:1259274/150‑1 (MQ=255)
                                                                                                                                           cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa        <  2:2354733/150‑1 (MQ=255)
                                                                                                                                           cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa        <  2:2571897/150‑1 (MQ=255)
                                                                                                                                           cagcaCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGaaa        <  2:960318/150‑1 (MQ=255)
                                                                                                                                               aCTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGAAAGTaa    >  2:4296385/1‑150 (MQ=255)
                                                                                                                                                 ttGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGAAAGTAAGc  <  2:4682546/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCTGGCGCTCGATGTCGGTGGCACCCATGTAGATCACGCAAGGCAGGCCGAAGCGCGCGGCGACAGTGGCGGTGGCGACACCGTGCATGCCGGCGCCGGTTTCGGCGATCAGGCGCTTTTTGCCCATGCGCTTGGCCAGCAGCACTTGGCCGATGCAGTTGTTCACCTTGTGCGCGCCGGTATGGTTGAGCTCTTCACGCTTGAAGAAAATCTTCGCTCCGCCGCAATGCTCGGTCAGGCGCTCGGCGAAGTACAGCGGGTTGGGGCGGCCGATGTAGTCGCGCTGAAAGTAAGC  >  NC_002947/93143‑93437

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: