Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,513,866 G→A A266V (GCA→GTA)  RPPX_24600 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,513,8660GA99.1% 1159.9 / NA 346A266V (GCA→GTA) RPPX_24600integrase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (172/171);  total (172/174)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 346 total aligned reads displayed.

GTCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGAT  >  CP009974/5513718‑5514009
                                                                                                                                                    |                                                                                                                                               
 tCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAg                                                                                                                                              <  2:4819306/151‑1 (MQ=255)
 tCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAg                                                                                                                                              <  2:2908547/151‑1 (MQ=255)
 tCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAg                                                                                                                                              <  1:3108932/151‑1 (MQ=255)
  ccTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGc                                                                                                                                             <  2:924933/151‑1 (MQ=255)
  ccTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGc                                                                                                                                             >  1:1643317/1‑151 (MQ=255)
  catGGTTCGAGTCATGTGGCGGAGGCAGGCGTGGCCGGCGGGTATCAATTCGATTGAGCGCTTCGGGGCTGCCTTGCGGCGCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGc                                                                                                                                             <  2:3357904/149‑1 (MQ=25)
   cTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACCAGTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACCGcc                                                                                                                                            >  2:5245089/1‑151 (MQ=255)
    tGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGccc                                                                                                                                           <  2:3945243/151‑2 (MQ=255)
    tGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGccc                                                                                                                                           >  1:6781389/1‑150 (MQ=255)
      ggtCGAGACTTATGGCGCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGa                                                                                                                                         <  2:5162519/149‑1 (MQ=255)
       ttCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACCGCCCGAg                                                                                                                                        >  2:7084284/1‑151 (MQ=255)
       ttCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAg                                                                                                                                        <  1:717185/151‑1 (MQ=255)
           agaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCt                                                                                                                                    >  1:880331/1‑151 (MQ=255)
             aCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGt                                                                                                                                  <  2:2584591/151‑1 (MQ=255)
                  tGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTcc                                                                                                                             >  2:3183847/1‑151 (MQ=255)
                  tGGCTAAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCATCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTcc                                                                                                                             <  2:1485458/151‑1 (MQ=255)
                      tCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAg                                                                                                                         <  1:2266540/151‑1 (MQ=255)
                       caggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAgg                                                                                                                        >  2:4167957/1‑151 (MQ=255)
                        aggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAggg                                                                                                                       >  2:3488859/1‑151 (MQ=255)
                        aggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAggg                                                                                                                       >  1:4132554/1‑151 (MQ=255)
                         ggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGc                                                                                                                      <  2:3138947/151‑1 (MQ=255)
                           caggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGTGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCGACAATTTGGACTGCCCCCCAAAAGATGGTCAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTc                                                                                                                    >  2:5191089/1‑151 (MQ=37)
                             ggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTg                                                                                                                                   >  2:54776/1‑134 (MQ=255)
                             ggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAg                                                                                                                  >  1:5621364/1‑151 (MQ=255)
                               aGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGc                                                                                                                <  2:4593421/151‑1 (MQ=255)
                                gTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGcc                                                                                                               <  1:3995495/151‑1 (MQ=255)
                                gTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGcc                                                                                                               <  2:4528351/151‑1 (MQ=255)
                                 tGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGccc                                                                                                              >  2:225446/1‑150 (MQ=255)
                                  ggCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTGGCGGCCCTTCGCGGGCTCGCGCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCAt                                                                                                             <  1:5452995/151‑1 (MQ=255)
                                    cAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTg                                                                                                           <  1:1866396/151‑1 (MQ=255)
                                     aGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGa                                                                                                          >  2:2716306/1‑151 (MQ=255)
                                        cGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACCGGGCTCAGGCCATTGAGtc                                                                                                       >  1:2554803/1‑151 (MQ=255)
                                           gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                    <  1:931848/151‑1 (MQ=255)
                                           gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                    <  1:4770462/151‑1 (MQ=255)
                                            tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCCCAGGGCTCAGGCCATTGAGTCGCAg                                                                                                   >  1:5309359/1‑151 (MQ=255)
                                            tATCAATACGATTGAATGAACCGGGGACGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                   <  2:1518463/151‑1 (MQ=255)
                                            tATCAATACGATGGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                   <  1:3176289/151‑1 (MQ=255)
                                               aaaTACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                                <  2:2120737/150‑1 (MQ=255)
                                                 gtACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGGTGGACCTTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGa                                                                                              <  2:7214164/150‑1 (MQ=21)
                                                   aCGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATg                                                                                            >  1:4317176/1‑151 (MQ=255)
                                                   aCGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATg                                                                                            <  1:4094344/151‑1 (MQ=255)
                                                    cGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGc                                                                                           >  2:4295213/1‑151 (MQ=255)
                                                     gATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCg                                                                                          >  1:6609571/1‑151 (MQ=255)
                                                       ttGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGAt                                                                                        >  1:2334295/1‑151 (MQ=255)
                                                           atgacCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATg                                                                                    <  2:5115945/151‑1 (MQ=255)
                                                           atgaaCCGGGGCCGCCTTGCGTCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATg                                                                                    <  2:5612014/151‑1 (MQ=255)
                                                           atgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATg                                                                                    >  1:6167106/1‑151 (MQ=255)
                                                                ccGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgt                                                                               <  2:4363788/151‑1 (MQ=255)
                                                                  ggggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATACCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtcg                                                                             >  2:6960152/1‑149 (MQ=255)
                                                                     gccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGtagtag                                                                          <  2:4675542/151‑1 (MQ=255)
                                                                     gccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGtagtag                                                                          <  2:4899618/151‑1 (MQ=255)
                                                                       cgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTg                                                                        >  1:639051/1‑151 (MQ=255)
                                                                       cgccTGGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGGCTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTg                                                                        <  1:6045125/151‑1 (MQ=255)
                                                                          cTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAAt                                                                     <  1:2735253/151‑1 (MQ=255)
                                                                           ttGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATg                                                                    <  2:231314/151‑1 (MQ=255)
                                                                             gCGGCCATTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTCGTAGTGAATGTa                                                                  <  2:638458/151‑1 (MQ=255)
                                                                              cGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTAc                                                                 >  1:4389632/1‑151 (MQ=255)
                                                                               ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCCGGCCATTGAGGCTCAGCCTGATGCGCTCATGCTTGTAGGAGGGAAAGTACt                                                                >  2:6114822/1‑151 (MQ=17)
                                                                               ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                >  1:2865630/1‑151 (MQ=255)
                                                                               ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                >  1:4708897/1‑151 (MQ=255)
                                                                                gCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTc                                                               <  1:5375652/151‑1 (MQ=255)
                                                                                 cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCa                                                              >  2:6591067/1‑151 (MQ=255)
                                                                                 cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCCATCATGGTTGTAGTAGTGAATGTACTCa                                                              >  1:4527487/1‑151 (MQ=255)
                                                                                  ccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCAt                                                             <  1:5824814/151‑1 (MQ=255)
                                                                                       gcgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGa                                                        >  1:2997840/1‑151 (MQ=255)
                                                                                        cgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAc                                                       >  1:5046540/1‑151 (MQ=255)
                                                                                         gggCTCGCCGTCGCGCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAcc                                                      <  2:6181056/151‑1 (MQ=255)
                                                                                          ggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAccc                                                     <  1:1506114/151‑1 (MQ=255)
                                                                                           gCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCg                                                    >  1:6756262/1‑151 (MQ=255)
                                                                                            cTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGc                                                   >  2:605471/1‑151 (MQ=255)
                                                                                              cGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctc                                                 <  1:3852080/151‑1 (MQ=255)
                                                                                               gCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctct                                                >  1:5581107/1‑151 (MQ=255)
                                                                                                cccTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCtctc                                               <  1:2037772/151‑1 (MQ=255)
                                                                                                 ccTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCa                                              <  1:7142106/151‑1 (MQ=255)
                                                                                                  cTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAg                                             >  2:6791484/1‑151 (MQ=255)
                                                                                                   tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGc                                            >  2:2446213/1‑151 (MQ=255)
                                                                                                    cGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGcc                                           >  1:3548882/1‑150 (MQ=255)
                                                                                                      cTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtca                                         >  1:5163396/1‑151 (MQ=255)
                                                                                                      cTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtca                                         >  2:3768002/1‑151 (MQ=255)
                                                                                                       tCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcat                                        <  2:5737307/151‑1 (MQ=255)
                                                                                                        cccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcatc                                       >  1:1341117/1‑151 (MQ=255)
                                                                                                        cccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcatc                                       >  1:1251005/1‑151 (MQ=255)
                                                                                                             aaTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGc                                  <  1:2697962/151‑1 (MQ=255)
                                                                                                                   gACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGa                            <  1:5531558/151‑1 (MQ=14)
                                                                                                                    aCTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGaa                           <  1:7029989/151‑1 (MQ=14)
                                                                                                                       gCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACg                        >  1:5494587/1‑151 (MQ=12)
                                                                                                                         cccccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCGCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCCCATCAATGCCCTCGAAACGCt                      >  2:1729080/1‑151 (MQ=11)
                                                                                                                           cccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTc                    <  1:1033836/151‑1 (MQ=12)
                                                                                                                            ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGCAACGCTTcg                   >  1:5246612/1‑150 (MQ=9)
                                                                                                                            ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCttac                   >  2:1894123/1‑149 (MQ=12)
                                                                                                                             cAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAg                  <  1:5554132/151‑1 (MQ=17)
                                                                                                                                 aGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAg              >  1:3422113/1‑151 (MQ=12)
                                                                                                                                  gTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGa             <  1:2133045/151‑1 (MQ=12)
                                                                                                                                   gtagACAGTTCTAACCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaa            <  2:3832939/148‑1 (MQ=9)
                                                                                                                                    gggACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaa           <  2:2540361/150‑1 (MQ=9)
                                                                                                                                       aCAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAAc        >  2:1264353/1‑151 (MQ=12)
                                                                                                                                        cAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGCAGAAAAACt       >  2:1560063/1‑151 (MQ=12)
                                                                                                                                          gTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCa     >  1:5898932/1‑151 (MQ=14)
                                                                                                                                             cTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGAt  <  2:2203331/151‑1 (MQ=11)
                                                                                                                                             cTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGAt  >  1:1788713/1‑151 (MQ=11)
                                                                                                                                                    |                                                                                                                                               
GTCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGAT  >  CP009974/5513718‑5514009

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