Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 970,240 T→G 6.2% F29C (TTC→TGC) ‡ tgt → queuine tRNA‑ribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947970,2400TG 6.2% 276.7 / 9.1 98F29C (TTC→TGC) ‡tgtqueuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (46/45);  new base G (3/3);  total (50/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 = 4.87e-01

TGTTCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCA  >  NC_002947/970091‑970384
                                                                                                                                                     |                                                                                                                                                   
tGTTCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCCGAACGCTTTGCACCGTGGAAACCCCGGCAtt                                                                                                                                                     >  1:941531/1‑150 (MQ=255)
  ttCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCa                                                                                                                                                   >  2:915319/1‑150 (MQ=255)
  ttCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCa                                                                                                                                                   >  1:536009/1‑150 (MQ=255)
     atcaCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAGCCCCGGCATTCATGc                                                                                                                                                >  1:572510/1‑150 (MQ=255)
      tcaCCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGcc                                                                                                                                              >  1:674365/1‑150 (MQ=255)
          ccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                         <  2:536009/150‑1 (MQ=255)
          ccGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTg                                                                                                                                         >  1:1550037/1‑150 (MQ=255)
           cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        <  1:350859/150‑1 (MQ=255)
           cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        <  1:281700/150‑1 (MQ=255)
           cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        >  1:1919195/1‑150 (MQ=255)
           cGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGGGTAGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTgg                                                                                                                                        >  1:961236/1‑150 (MQ=255)
            gCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTggg                                                                                                                                       >  2:1797345/1‑150 (MQ=255)
               aTCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGAc                                                                                                                                    <  2:145724/150‑1 (MQ=255)
                   ggCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGCCGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                <  2:641074/150‑1 (MQ=255)
                   ggCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                <  2:573457/150‑1 (MQ=255)
                   ggCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                <  1:1540247/150‑1 (MQ=255)
                     cgcCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATgg                                                                                                                              <  2:1584403/150‑1 (MQ=255)
                     cgcCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATgg                                                                                                                              <  2:1550037/150‑1 (MQ=255)
                           gcgGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCTGCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGt                                                                                                                       <  2:914150/150‑1 (MQ=255)
                               gCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAg                                                                                                                   >  2:427865/1‑150 (MQ=255)
                               gCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAg                                                                                                                   <  1:1297262/150‑1 (MQ=255)
                                 ccGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                <  2:1351860/136‑1 (MQ=255)
                                 ccGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                >  1:1351860/1‑136 (MQ=255)
                                    aGGATCAACCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCAt                                                                                                              <  1:32515/150‑1 (MQ=255)
                                      gATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCAtgc                                                                                                            >  2:1745979/1‑150 (MQ=255)
                                      gATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCAtgc                                                                                                            >  2:33056/1‑150 (MQ=255)
                                      gATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGAACTATGGCACGGTCAAGGGCAtgc                                                                                                            <  2:1682613/150‑1 (MQ=255)
                                      gATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCAACGTGGAAACCCCGCCA‑‑CACGCTGGGTGGGGGACCTATGGCACGGTCAAGGGCAtgc                                                                                                            >  2:152257/1‑150 (MQ=255)
                                           aGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCa                                                                                                       >  1:1661577/1‑150 (MQ=255)
                                              aTGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgc                                                                                                    >  2:1619731/1‑150 (MQ=255)
                                              aTGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgc                                                                                                    <  2:1345090/150‑1 (MQ=255)
                                                gAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGa                                                                                                  >  2:1917250/1‑150 (MQ=255)
                                                 aGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGAc                                                                                                 <  1:1166755/150‑1 (MQ=255)
                                                   tCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAgg                                                                                                                  >  2:48801/1‑131 (MQ=255)
                                                   tCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAgg                                                                                                                  <  1:48801/131‑1 (MQ=255)
                                                   tCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACAt                                                                                               >  1:549436/1‑150 (MQ=255)
                                                    cGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATc                                                                                              <  2:1689701/150‑1 (MQ=255)
                                                       aCCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGac                                                                                           >  2:433230/1‑149 (MQ=255)
                                                       aCCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAg                                                                                           <  2:1074149/150‑1 (MQ=255)
                                                           gTCGAATGTCCTTCGAACTGCTGGCCAGCGACGGCAAGGCCAGACGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCa                                                                                       >  2:21651/1‑150 (MQ=255)
                                                            tCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgc                                                                                                    >  2:804459/1‑136 (MQ=255)
                                                            tCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCAcgc                                                                                                    <  1:804459/136‑1 (MQ=255)
                                                             cGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                <  2:1836306/108‑1 (MQ=255)
                                                             cGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTAt                                                                                                                                >  1:1836306/1‑108 (MQ=255)
                                                               aaTGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCgg                                                                                   >  2:1681084/1‑150 (MQ=255)
                                                                  gTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGgcgc                                                                                <  2:1628793/150‑1 (MQ=255)
                                                                  gTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGgcgc                                                                                <  2:1776216/150‑1 (MQ=255)
                                                                     cTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACACCGGCATTc                                                                                                                                                    <  2:1532402/82‑1 (MQ=255)
                                                                     cTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACACCGGCATTc                                                                                                                                                    >  1:1532402/1‑82 (MQ=255)
                                                                      ttCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGAGAATCGCGCAGCAACGAATGGCAGGCGGCCACGCGACATCGAGGCCATCGGCGCCgag                                                                            >  2:988544/1‑150 (MQ=255)
                                                                       tCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCgaga                                                                           <  1:915319/150‑1 (MQ=255)
                                                                             tgctgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCa                                                                                                                                                   <  2:1410496/75‑1 (MQ=255)
                                                                             tgctgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCa                                                                                                                                                   >  1:1410496/1‑75 (MQ=255)
                                                                             tgctgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATcc                                                                     >  1:768415/1‑150 (MQ=255)
                                                                                tgGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCg                                                                  >  2:847333/1‑150 (MQ=255)
                                                                                 ggCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCgg                                                                 <  1:1619731/150‑1 (MQ=255)
                                                                                  gCCACCGACGGCAAGGCCCGTCGGTGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGc                                                                <  2:1266987/150‑1 (MQ=255)
                                                                                  gCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGAAACCTATGGCACGGTCAAGGGCATTTTGCCACGCGGCATCGAGGCCATCGGCGCCGAGATGATCCTCGGc                                                                <  1:515749/150‑1 (MQ=255)
                                                                                    cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGGCAGCTGGCATCGAGGCCATCGGCGCCGAGATGATCCTCGGCaa                                                              >  1:959504/1‑150 (MQ=255)
                                                                                    cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCaa                                                              >  1:242126/1‑150 (MQ=255)
                                                                                    cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCggg                                                                                                                                           <  2:345096/75‑2 (MQ=255)
                                                                                    cACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATGAATGCCggg                                                                                                                                           >  1:345096/1‑74 (MQ=255)
                                                                                    cACCGACGGCAAGGCCCGTCGGGGCCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCaa                                                              >  2:1285819/1‑150 (MQ=255)
                                                                                      ccGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAaca                                                            <  1:33056/150‑1 (MQ=255)
                                                                                         aCGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCt                                                         >  2:514433/1‑150 (MQ=255)
                                                                                           ggCAAGCCCCGTCGGGGTCGCATCACCGCCCCCCGTGGCACCGTGGAATCCCCGGCATTCGTGCCGGTGGGGACCGAGGGCACGGTCAAGGGCATGCGGCCACGCGACACCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTc                                                       <  2:61272/150‑1 (MQ=255)
                                                                                            gCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTcc                                                      >  2:1611211/1‑150 (MQ=255)
                                                                                              aaGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCAc                                                    <  1:1466204/150‑1 (MQ=255)
                                                                                                   ccGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCCGGATCAAGGGCATGCTGCCACGCGACATCGCGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgtg                                               >  1:1068659/1‑150 (MQ=255)
                                                                                                    cGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATg                                                                                                                               <  2:913147/70‑1 (MQ=255)
                                                                                                    cGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATg                                                                                                                               >  1:913147/1‑70 (MQ=255)
                                                                                                           gTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCTGGGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTc                                       >  1:901705/1‑150 (MQ=255)
                                                                                                              gCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCgg                                    >  1:773725/1‑150 (MQ=255)
                                                                                                              gCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCgg                                    <  1:1611211/150‑1 (MQ=255)
                                                                                                                 tCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCa                                 <  2:202877/150‑1 (MQ=255)
                                                                                                                 tCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCa                                 <  1:427865/150‑1 (MQ=255)
                                                                                                                    ccTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCg                              >  2:1372652/1‑150 (MQ=255)
                                                                                                                     cTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCTTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGa                             >  2:174242/1‑150 (MQ=255)
                                                                                                                       tCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAgg                           <  1:679067/150‑1 (MQ=255)
                                                                                                                          cACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGa                        <  2:768415/150‑1 (MQ=255)
                                                                                                                           aCGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTgg                                              >  1:22807/1‑127 (MQ=255)
                                                                                                                           aCGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTgg                                              <  2:22807/127‑1 (MQ=255)
                                                                                                                           aCGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGAt                       <  1:99937/150‑1 (MQ=255)
                                                                                                                           aCGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGAt                       >  2:684081/1‑150 (MQ=255)
                                                                                                                                 cACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCaagaa                 >  2:290198/1‑150 (MQ=255)
                                                                                                                                  aCCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCaagaag                >  1:556022/1‑150 (MQ=255)
                                                                                                                                    cGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGca              >  1:67060/1‑150 (MQ=255)
                                                                                                                                     gTGGAAACCACGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGcac             >  2:1404891/1‑150 (MQ=255)
                                                                                                                                       ggAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACaa           <  2:67060/150‑1 (MQ=255)
                                                                                                                                       ggAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACaa           <  1:1745979/150‑1 (MQ=255)
                                                                                                                                           ccaccGGCATGAATGCCGGGGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAacac                                                           >  2:1742379/4‑98 (MQ=255)
                                                                                                                                           ccaccGGCATGAATGCCGGGGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgtg                                               <  2:1340521/107‑1 (MQ=255)
                                                                                                                                           ccaccGGCATGAATGCCGGGGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCtgtg                                               >  1:1340521/4‑110 (MQ=255)
                                                                                                                                           ccaccGGCATGAATGCCGGGGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACACCGAGGCCATCGGCGCCGAGATGATCCTCGGCAacac                                                           <  1:1742379/95‑1 (MQ=255)
                                                                                                                                           aCCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGc       <  1:1911520/150‑1 (MQ=255)
                                                                                                                                           aCCCCGGCATTCATGCCGGTGAGGTCCTCTGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCCACACCTTCAACCTGCGGCTGCATGCGCGCGCCGTGGTGAGCAAGAAGCACAACGGc       <  2:961236/150‑1 (MQ=255)
                                                                                                                                               cGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGc   >  2:629255/1‑150 (MQ=255)
                                                                                                                                               cGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGc   >  2:295572/1‑150 (MQ=255)
                                                                                                                                                ggCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCa  <  2:1919195/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTTCATCACCCGCAATCCGGCGCCACGCGGGCCCGAGGATCAAGCATGAGTCGCACCTGTCGAATGTCCTTCGAACTGCTGGCCACCGACGGCAAGGCCCGTCGGGGTCGCATCACCTTCCCACGTGGCACCGTGGAAACCCCGGCATTCATGCCGGTGGGGACCTATGGCACGGTCAAGGGCATGCTGCCACGCGACATCGAGGCCATCGGCGCCGAGATGATCCTCGGCAACACCTTCCACCTGTGGCTGCGTCCGGGCACCGAGGTGATCAAGAAGCACAACGGCCTGCA  >  NC_002947/970091‑970384

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