Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,337,753 A→T 5.1% D305V (GAC→GTC)  PP_2934 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,337,7530AT 5.1% 317.0 / 4.2 100D305V (GAC→GTC) PP_2934hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (47/47);  new base T (3/2);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 102 total aligned reads displayed.

CCGTAGACCGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAG  >  NC_002947/3337606‑3337902
                                                                                                                                                   |                                                                                                                                                     
ccGTAGACCGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACg                                                                                                                                                     >  2:2367901/1‑150 (MQ=255)
  gTAGACCGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGc                                                                                                                                                   >  1:2248535/1‑150 (MQ=255)
    aGACCGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGt                                                                                                                                                 <  1:2323128/150‑1 (MQ=255)
       ccGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGTCGCTGTGCAAGGACGCCGACGGCGTCGa                                                                                                                                              <  1:873281/150‑1 (MQ=255)
       ccGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGa                                                                                                                                              <  2:1511537/150‑1 (MQ=255)
       ccGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGa                                                                                                                                              <  1:759603/150‑1 (MQ=255)
       ccGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCAGGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGTCGGCGTcct                                                                                                                                              <  2:1180495/150‑3 (MQ=255)
             tGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGAccgccg                                                                                                                                        <  2:1567571/150‑1 (MQ=255)
                         ggCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCcatca                                                                                                                            >  1:410280/1‑150 (MQ=255)
                            ccGCAAACTGCCGGTCAGACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTc                                                                                                                         >  2:1236846/1‑150 (MQ=255)
                            ccGCAAACTGACGGTCACACCGGACTACCGCCTCTTGGCAGCTGCAAAACGTGGCGGATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGCACGCCGACTGCGTCGACGGGCGGTATATCcatcatgc                                                                                                                         >  2:339122/1‑148 (MQ=255)
                            aaccAAACTGCCGTTCGCCCCGGACTAGCGCCTGTTGGCCGTTTGCGGGCACGTCGTATTCATGGCCCACTGCTATGCCGAGGAGTACGTGCGAATGCCGGCGCTGTGCGACGACGCCGACTTCGTCGACCGCCGGTATATCCAGCATTc                                                                                                                         <  2:665845/147‑1 (MQ=255)
                             cGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCg                                                                                                                        >  1:1568309/1‑150 (MQ=255)
                              gCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGGCGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGc                                                                                                                       <  2:1983905/150‑1 (MQ=255)
                                  aCTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCa                                                                                                                   >  1:3414/1‑150 (MQ=255)
                                   cTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGAATTTCTTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAg                                                                                                                  <  1:262946/150‑1 (MQ=255)
                                      ccGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGGATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGGGCAAGGACGCCGACGGCGGCGACCGCCGGTATATCCATCATTCGCTGCAGCga                                                                                                               <  1:322612/150‑1 (MQ=255)
                                          tCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGaaa                                                                                                           >  1:1511829/1‑150 (MQ=255)
                                              acCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGc                                                                                                       >  1:801186/1‑150 (MQ=255)
                                              acCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGc                                                                                                       >  2:359834/1‑150 (MQ=255)
                                              acCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGc                                                                                                       <  1:2367901/150‑1 (MQ=255)
                                               ccGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCa                                                                                                      <  1:1236846/150‑1 (MQ=255)
                                                cGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGCAACGGCAg                                                                                                     <  2:2212050/150‑1 (MQ=255)
                                                     tACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGttc                                                                                                <  2:733874/150‑1 (MQ=255)
                                                     tACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGttc                                                                                                <  1:380485/150‑1 (MQ=255)
                                                       ccGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGttctt                                                                                              <  2:801186/150‑1 (MQ=255)
                                                        cGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGttcttc                                                                                             <  2:1749403/150‑1 (MQ=255)
                                                        cGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCACTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCCGCGttcttc                                                                                             >  2:660959/1‑150 (MQ=255)
                                                         gCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGt                                                                                                                                                 >  2:466138/1‑97 (MQ=255)
                                                         gCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGt                                                                                                                                                 <  1:466138/97‑1 (MQ=255)
                                                         gCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCgg                                                                                                                                       <  1:1343227/107‑1 (MQ=255)
                                                         gCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCgg                                                                                                                                       >  2:1343227/1‑107 (MQ=255)
                                                          ccTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAg                                                                                           >  1:1905953/1‑150 (MQ=255)
                                                          ccTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAg                                                                                           >  1:836347/1‑150 (MQ=255)
                                                           cTGTTGGCCGGTGCCGGGTATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGc                                                                                          <  1:2015151/150‑1 (MQ=255)
                                                           cTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCAAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGc                                                                                          <  1:1122952/150‑1 (MQ=255)
                                                            tGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGcc                                                                                         <  2:2230591/150‑1 (MQ=255)
                                                              ttGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGc                                                                                                                                                   <  1:603518/90‑1 (MQ=255)
                                                              ttGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGc                                                                                                                                                   >  2:603518/1‑90 (MQ=255)
                                                                 gCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCt                                                                                                                      <  1:52283/116‑1 (MQ=255)
                                                                 gCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCt                                                                                                                      >  2:52283/1‑116 (MQ=255)
                                                                   cGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCt                                                                                  >  2:2357066/1‑150 (MQ=255)
                                                                     gTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTgg                                                                                <  1:2037614/150‑1 (MQ=255)
                                                                          gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCg                                                                           >  2:160559/1‑150 (MQ=255)
                                                                           ggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGc                                                                          <  2:3414/150‑1 (MQ=255)
                                                                            gCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGAAGCCAGACACGACAGCGTTCTTCAGCCAGACGCTGGGGTCGcc                                                                         >  2:1123934/1‑150 (MQ=255)
                                                                             cATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCa                                                                        <  1:359834/150‑1 (MQ=255)
                                                                               tttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgc                                                                      >  1:72298/1‑150 (MQ=255)
                                                                                 tCGAATGCATTGCCCATTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgcgc                                                                    <  2:262484/150‑1 (MQ=255)
                                                                                  cGTATTCATGGCCCAAGGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGNAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGCGACCCCAGCGcc                                                                   <  1:1442420/150‑1 (MQ=255)
                                                                                   gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCC‑GGCGCTGGCGACCCCAGCGcctg                                                                 <  1:534673/150‑3 (MQ=255)
                                                                                    tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCgg                                                                 <  2:564013/150‑1 (MQ=255)
                                                                                    tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCgg                                                                 >  1:1521198/1‑150 (MQ=255)
                                                                                      ttCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCt                                                               >  1:2029577/1‑150 (MQ=255)
                                                                                        cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGa                                                             <  2:3089/150‑1 (MQ=255)
                                                                                          tGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGc                                                                                                                       >  2:432335/1‑90 (MQ=255)
                                                                                          tGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGc                                                                                                                       <  1:432335/90‑1 (MQ=255)
                                                                                            gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGt                                                         >  2:228855/1‑150 (MQ=255)
                                                                                            gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGt                                                         >  2:1318422/1‑150 (MQ=255)
                                                                                              ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAATGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTc                                                       >  2:1069523/1‑150 (MQ=255)
                                                                                               cACTGTCATGCCGAGCAGTACGCGCGAATGCCGGCGCTGGGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGTGCGTTCa                                                      <  2:1344389/150‑1 (MQ=255)
                                                                                               cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCAtt                                                                                                                          >  2:2256922/1‑82 (MQ=255)
                                                                                               cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCAtt                                                                                                                          <  1:2256922/82‑1 (MQ=255)
                                                                                               cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCa                                                      <  1:743448/150‑1 (MQ=255)
                                                                                                aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCg                                                                  >  1:1599488/1‑137 (MQ=255)
                                                                                                aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCg                                                                  <  2:1599488/137‑1 (MQ=255)
                                                                                                aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAg                                                     <  1:1123934/150‑1 (MQ=255)
                                                                                                 cTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGc                                                    <  2:2036110/150‑1 (MQ=255)
                                                                                                   gCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTg                                                  <  1:228855/150‑1 (MQ=255)
                                                                                                       tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGAc                                              >  2:1817228/1‑150 (MQ=255)
                                                                                                           gAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGgc                                          >  1:1089931/1‑150 (MQ=255)
                                                                                                                 tatgcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGGCGGCTGAGCGTTCAGCTGCGACGGGcgcgccg                                    >  1:953352/4‑150 (MQ=255)
                                                                                                                  aCGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGcgcgccgc                                   <  2:1089931/150‑1 (MQ=255)
                                                                                                                    gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGc                                                                                                                       <  2:1651244/64‑1 (MQ=255)
                                                                                                                    gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGc                                                                                                                       >  1:1651244/1‑64 (MQ=255)
                                                                                                                    gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCa                                                      <  1:43016/129‑1 (MQ=255)
                                                                                                                    gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCa                                                      >  2:43016/1‑129 (MQ=255)
                                                                                                                    gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                 >  2:899390/1‑150 (MQ=255)
                                                                                                                    gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                 >  1:1799052/1‑150 (MQ=255)
                                                                                                                        gAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                             >  1:649909/1‑150 (MQ=255)
                                                                                                                           tGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                          >  2:222648/1‑150 (MQ=255)
                                                                                                                             ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCTATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                        <  1:2068151/150‑1 (MQ=255)
                                                                                                                             ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                        >  2:1503000/1‑150 (MQ=255)
                                                                                                                              cGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                       >  2:1630088/1‑150 (MQ=255)
                                                                                                                               ggCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACg                      >  1:915163/1‑150 (MQ=255)
                                                                                                                                gcgcTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTggcgac                                                                            >  1:629865/1‑91 (MQ=255)
                                                                                                                                gcgcTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTggcgac                                                                            <  2:629865/95‑5 (MQ=255)
                                                                                                                                    tgtgCAAGGACGCCGGCGGAGACCAGCCGCTGTATATCCATCATTAGCTACAGCGCGAACCCACAGCTGGCTGCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                 <  1:44472/150‑1 (MQ=255)
                                                                                                                                    tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                 <  1:1318422/150‑1 (MQ=255)
                                                                                                                                      tgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGcc               <  2:72298/150‑1 (MQ=255)
                                                                                                                                          tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCt                                                                                  >  1:2258005/3‑79 (MQ=255)
                                                                                                                                          tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCt                                                                                  <  2:2258005/77‑1 (MQ=255)
                                                                                                                                          tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAcc           >  2:2325104/3‑150 (MQ=255)
                                                                                                                                          tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGACGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACGACGCTAGGCCGAcc           >  2:1126018/3‑150 (MQ=255)
                                                                                                                                           ggACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAccc          >  2:1362265/1‑150 (MQ=255)
                                                                                                                                                 cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAAAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACGACGTTAGGCCGACCCCGGCCa    >  2:1239473/1‑150 (MQ=255)
                                                                                                                                                 cGACGGCGTCGACCGCCGGTACATACACGAATCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa    >  1:65048/1‑150 (MQ=255)
                                                                                                                                                  gACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCaa   >  1:2032291/1‑150 (MQ=255)
                                                                                                                                                   aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCggactac  <  1:1069523/150‑6 (MQ=255)
                                                                                                                                                   aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg  <  2:654210/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGTAGACCGTAATGCCGAGGCCCTGGCCCGCAAACTGCCGGTCACACCGGACTACCGCCTGTTGGCCGGTGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCCGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAG  >  NC_002947/3337606‑3337902

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