Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,019,126 C→A 6.7% L120M (CTG→ATG) ‡ PP_1797 → HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,019,1260CA 6.7% 333.6 / 8.2 105L120M (CTG→ATG) ‡PP_1797HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base C (46/52);  new base A (4/3);  total (50/55)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

Only 100 of 108 total aligned reads displayed.

AGGGTGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCGT  >  NC_002947/2018977‑2019275
                                                                                                                                                     |                                                                                                                                                     
aGGGTGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGCGCCGCAGCTCTGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGcc                                                                                                                                                       >  1:541613/1‑150 (MQ=255)
    tGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCg                                                                                                                                                   <  2:2256900/150‑1 (MQ=255)
            tCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcctctt                                                                                                                                           <  1:1488910/150‑5 (MQ=255)
                  ggCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACg                                                                                                                                     <  2:249273/149‑1 (MQ=255)
                    cTGCATGCCCGCGATGGGGACCTTTCTGAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGc                                                                                                                                   >  1:2158012/1‑150 (MQ=255)
                     tGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCt                                                                                                                                  >  2:686364/1‑150 (MQ=255)
                       cATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTg                                                                                                                                <  2:2322783/150‑1 (MQ=255)
                        aTGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTgg                                                                                                                               <  2:801354/150‑1 (MQ=255)
                         tGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGa                                                                                                                              <  1:1150139/150‑1 (MQ=255)
                         tGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGa                                                                                                                              <  1:2257137/150‑1 (MQ=255)
                         tGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcctcttt                                                                                                                                          <  2:809616/138‑6 (MQ=255)
                         tGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcctcttt                                                                                                                                          >  1:809616/1‑133 (MQ=255)
                            ccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTc                                                                                                                           >  2:124187/1‑150 (MQ=255)
                            ccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCACGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTc                                                                                                                           <  1:551656/150‑1 (MQ=255)
                             cgcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcctctt                                                                                                                                           <  2:1974043/133‑5 (MQ=255)
                             cgcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcctctt                                                                                                                                           >  1:1974043/1‑129 (MQ=255)
                              gcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGa                                                                                                                         <  1:1618814/150‑1 (MQ=255)
                              gcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGa                                                                                                                         >  1:411570/1‑150 (MQ=255)
                                gATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATc                                                                                                                       <  1:1500713/150‑1 (MQ=255)
                                    gggACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATc                                                                                                                       <  1:517913/146‑1 (MQ=255)
                                    gggACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATc                                                                                                                       >  2:517913/1‑146 (MQ=255)
                                    gggACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGt                                                                                                                   <  1:1917288/150‑1 (MQ=255)
                                     ggACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTc                                                                                                                  >  1:1402063/1‑150 (MQ=255)
                                      gACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCg                                                                                                                 <  2:404406/150‑1 (MQ=255)
                                      gACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAACGTAAAATTTTGCACTCGAGCCCGTCg                                                                                                                 >  2:1132601/1‑150 (MQ=255)
                                         cTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTg                                                                                                              <  1:1478441/150‑1 (MQ=255)
                                         cTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTg                                                                                                              >  2:2336937/1‑150 (MQ=255)
                                          tGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGa                                                                                                                         <  1:1102572/138‑1 (MQ=255)
                                          tGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGa                                                                                                                         >  2:1102572/1‑138 (MQ=255)
                                              gaagaaGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAAt                                                                                                         <  2:2407738/150‑1 (MQ=255)
                                              gaagaaGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAAt                                                                                                         <  1:2184273/150‑1 (MQ=255)
                                              gaagaaGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAAt                                                                                                         <  2:432191/150‑1 (MQ=255)
                                                 gaagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGTTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGa                                                                                                      <  1:124187/150‑1 (MQ=255)
                                                 gaagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGCGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGa                                                                                                      >  1:373441/1‑150 (MQ=255)
                                                 gaagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGa                                                                                                      >  2:1832196/1‑150 (MQ=255)
                                                  aagGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTCTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                     <  1:1352008/150‑1 (MQ=255)
                                                   agGGCGATCCGCTTATTGTGCTGGATGAGTCGCAGCGCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAg                                                                                                    >  2:440903/1‑150 (MQ=255)
                                                    gggCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCACAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAgg                                                                                                   >  2:38338/1‑150 (MQ=255)
                                                             gCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTc                                                                                          <  1:1804685/150‑1 (MQ=255)
                                                              cTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCt                                                                                         >  1:2070475/1‑150 (MQ=255)
                                                                gATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGa                                                                                       <  1:1832196/150‑1 (MQ=255)
                                                                  ttGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGc                                                                                                                                                    <  2:589031/87‑1 (MQ=255)
                                                                  ttGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGc                                                                                                                                                    >  1:589031/1‑87 (MQ=255)
                                                                  ttGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGAcc                                                                                     <  1:2182090/150‑1 (MQ=255)
                                                                     tgCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACcgcg                                                                                  >  2:1376114/1‑150 (MQ=255)
                                                                      gCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                                 <  1:2452183/102‑1 (MQ=255)
                                                                      gCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                                 >  2:2452183/1‑102 (MQ=255)
                                                                            tGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCaaga                                                                           >  2:489082/1‑150 (MQ=255)
                                                                                       tCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAg                                                                <  1:489082/150‑1 (MQ=255)
                                                                                       tCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAg                                                                >  1:1789669/1‑150 (MQ=255)
                                                                                         aGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                              <  2:2070475/150‑1 (MQ=255)
                                                                                         aGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGAACGCGCGCAAGAAGCGCTGGCAGGt                                                              <  1:1132601/150‑1 (MQ=255)
                                                                                           cgcCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                     >  1:86035/1‑109 (MQ=255)
                                                                                           cgcCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGaa                                                                                                     <  2:86035/109‑1 (MQ=255)
                                                                                           cgcCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            <  1:1376114/150‑1 (MQ=255)
                                                                                            gcCCAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCATTCCCTAGGGTCGCTGGCAGGTGag                                                           <  2:1402063/150‑1 (MQ=255)
                                                                                              cGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagca                                                         >  2:641682/1‑150 (MQ=255)
                                                                                                  tATGAGTGCACCCGCAACCAAATCATTATTGCCCGTTAGCAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTAGATCCCGTCGGTGATTAGTCTGCGCCCCCGGGCTGTCGGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                     <  1:1004540/150‑1 (MQ=255)
                                                                                                  tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            >  1:1071765/1‑143 (MQ=255)
                                                                                                  tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            <  2:1071765/143‑1 (MQ=255)
                                                                                                  tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                     <  1:641682/150‑1 (MQ=255)
                                                                                                  tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                     >  1:377541/1‑150 (MQ=255)
                                                                                                  tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGCCACGTCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                     <  1:2365142/150‑1 (MQ=255)
                                                                                                      aGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGgtgt                                                 <  1:686364/150‑1 (MQ=255)
                                                                                                          cACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            <  2:1617904/135‑1 (MQ=255)
                                                                                                          cACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            >  1:1617904/1‑135 (MQ=255)
                                                                                                             ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTgg                                                                   <  2:2305229/125‑1 (MQ=255)
                                                                                                             ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTgg                                                                   >  1:2305229/1‑125 (MQ=255)
                                                                                                             ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagc                                                          >  2:1437366/1‑134 (MQ=255)
                                                                                                             ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagc                                                          <  1:1437366/134‑1 (MQ=255)
                                                                                                             ccGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGc                                          >  2:1949256/1‑150 (MQ=255)
                                                                                                              cGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTgg                                                                   <  1:364844/124‑1 (MQ=255)
                                                                                                              cGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTgg                                                                   >  2:364844/1‑124 (MQ=255)
                                                                                                               gCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGccc                                        <  1:1949256/150‑1 (MQ=255)
                                                                                                                   ccaacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCa                                    >  2:326690/1‑150 (MQ=255)
                                                                                                                    caacTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATACCGTCTGTGATAATGGAAGTTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCAc                                   >  2:270729/1‑150 (MQ=255)
                                                                                                                           ttgttgCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGt                            >  1:1852657/1‑150 (MQ=255)
                                                                                                                                 cccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGTTCTGCACGTGGTTAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTgg                      <  2:2158012/150‑1 (MQ=255)
                                                                                                                                  ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTAc                                                                                                <  1:75400/75‑1 (MQ=255)
                                                                                                                                  ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTAc                                                                                                >  2:75400/1‑75 (MQ=255)
                                                                                                                                  ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTggg                     >  1:2016418/1‑150 (MQ=255)
                                                                                                                                  ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTggg                     >  1:541268/1‑150 (MQ=255)
                                                                                                                                   cGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCa                                    <  2:1941371/134‑1 (MQ=255)
                                                                                                                                       ataAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGaaa                >  2:2207102/1‑150 (MQ=255)
                                                                                                                                         aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCCAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtc              >  1:1175942/1‑150 (MQ=255)
                                                                                                                                         aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGAGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtc              >  1:1736896/1‑150 (MQ=255)
                                                                                                                                          aaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtct             >  1:2144621/1‑150 (MQ=255)
                                                                                                                                           ctcggCACGCAGGCGTGCCGAGAGTGACGGCTTGGATTCGATCCCGTCGGGGATTATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtctc            >  1:895237/4‑150 (MQ=255)
                                                                                                                                            gagGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg           <  1:330515/150‑1 (MQ=255)
                                                                                                                                            gagGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg           <  1:2207102/150‑1 (MQ=255)
                                                                                                                                             agGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCgg          <  2:2016418/150‑1 (MQ=255)
                                                                                                                                               gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            <  1:1046794/98‑1 (MQ=255)
                                                                                                                                               gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGa                                                            >  2:1046794/1‑98 (MQ=255)
                                                                                                                                               gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTc                              <  2:1073325/128‑1 (MQ=255)
                                                                                                                                               gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTc                              >  1:1073325/1‑128 (MQ=255)
                                                                                                                                               gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTc        >  1:537545/1‑150 (MQ=255)
                                                                                                                                               gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTc        >  1:353177/1‑150 (MQ=255)
                                                                                                                                                   gCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAAcc    >  1:489444/1‑150 (MQ=255)
                                                                                                                                                     cTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCGt  >  2:590193/1‑150 (MQ=255)
                                                                                                                                                     cTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCGt  <  1:1806279/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGTGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACCGT  >  NC_002947/2018977‑2019275

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