Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,513,872 T→C Q264R (CAG→CGG)  RPPX_24600 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,513,8720TC99.1% 1020.2 / ‑6.1 332Q264R (CAG→CGG) RPPX_24600integrase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (165/164);  total (165/167)
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 332 total aligned reads displayed.

GTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGATTTGAG  >  CP009974/5513724‑5514014
                                                                                                                                                    |                                                                                                                                              
gTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGa                                                                                                                                              <  2:1962508/151‑1 (MQ=255)
  tCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGc                                                                                                                                            <  1:6036630/151‑1 (MQ=255)
   ggagaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGccc                                                                                                                                          <  2:2655279/150‑1 (MQ=255)
       aCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAg                                                                                                                                             <  2:1116792/145‑1 (MQ=255)
          tATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTAtt                                                                                                                                    <  1:2405866/151‑1 (MQ=255)
           aTGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAAATTGGACTGCCCCCCAAAAGCTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTc                                                                                                                                   >  2:1848643/1‑151 (MQ=255)
            tGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTcc                                                                                                                                  >  2:3183449/1‑151 (MQ=255)
            tGGCTAAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCATCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTcc                                                                                                                                  <  2:1485458/151‑1 (MQ=255)
             ggCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCca                                                                                                                                 <  1:6671702/151‑1 (MQ=255)
                 caggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAgg                                                                                                                             >  2:4167957/1‑151 (MQ=255)
                  aggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAcgg                                                                                                                            >  2:2438762/1‑151 (MQ=255)
                     caggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTc                                                                                                                         >  2:7159138/1‑151 (MQ=255)
                     caggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTc                                                                                                                         >  2:5189661/1‑151 (MQ=255)
                       ggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAg                                                                                                                       >  1:5621364/1‑151 (MQ=255)
                            ggCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTGGCGGCCCTTCGCGGGCTCGCGCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCAt                                                                                                                  <  1:5452995/151‑1 (MQ=255)
                               aGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGCTCTAAGCTACAGCCCGAGCCCTGTATTCCACCGGGCTCAGGCCCTTGa                                                                                                               >  2:6581423/1‑151 (MQ=255)
                                 gCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGt                                                                                                             <  1:1762043/151‑1 (MQ=255)
                                   gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtct                                                                                                           <  1:7113971/151‑1 (MQ=255)
                                   gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtct                                                                                                           >  2:3454345/1‑151 (MQ=255)
                                   gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtcc                                                                                                           >  2:5413783/1‑150 (MQ=255)
                                    ggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtctc                                                                                                          <  2:4736629/151‑1 (MQ=255)
                                    ggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtctc                                                                                                          >  2:2230954/1‑151 (MQ=255)
                                     gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                         <  1:4770462/151‑1 (MQ=255)
                                     gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                         <  1:931848/151‑1 (MQ=255)
                                      tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCCCAGGGCTCAGGCCATTGAGTCGCAg                                                                                                        >  1:5309359/1‑151 (MQ=255)
                                         cAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                                     <  1:2719315/151‑1 (MQ=255)
                                         cAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCGCGGCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                                     <  2:5212820/151‑1 (MQ=255)
                                           gtACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGGTGGACCTTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGa                                                                                                   <  2:7214164/150‑1 (MQ=21)
                                             aCGATGGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATg                                                                                                 <  2:6272566/151‑1 (MQ=255)
                                              cGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGc                                                                                                <  1:6121516/151‑1 (MQ=255)
                                              cGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGc                                                                                                >  2:4295213/1‑151 (MQ=255)
                                              cGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGAACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGc                                                                                                <  2:3074607/151‑1 (MQ=255)
                                               gATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCg                                                                                               >  1:6609571/1‑151 (MQ=255)
                                                aTTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGa                                                                                              >  1:55489/1‑151 (MQ=255)
                                                    aatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                          <  1:3735183/151‑1 (MQ=255)
                                                     atgaaCCGGGGCCGCCTTGCGTCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATg                                                                                         <  2:5612014/151‑1 (MQ=255)
                                                     atgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATg                                                                                         <  2:5114730/151‑1 (MQ=255)
                                                     atgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACACCCCGCGCCCTGTATTCCACCCGGCCCACGCCATTGAGTCTCCCCCTGATGCGCTCATg                                                                                         >  1:6336198/1‑151 (MQ=18)
                                                            ggggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGACCATAGTCgtag                                                                                  >  2:4438877/1‑151 (MQ=18)
                                                             gggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtagt                                                                                 <  1:5740776/151‑1 (MQ=255)
                                                              ggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTgtagta                                                                                >  2:2255162/1‑151 (MQ=255)
                                                               gccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGtagtag                                                                               <  2:4899618/151‑1 (MQ=255)
                                                                ccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTagtagt                                                                              <  1:1908024/151‑1 (MQ=255)
                                                                       gCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTa                                                                       <  1:4179467/151‑1 (MQ=255)
                                                                       gCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTa                                                                       <  2:3226612/151‑1 (MQ=255)
                                                                        cGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTAc                                                                      >  1:3935666/1‑151 (MQ=255)
                                                                         ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACCGGGCTCAGGCCATTGAGCCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                     >  2:6782900/1‑151 (MQ=255)
                                                                         ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCCGGCCATTGAGGCTCAGCCTGATGCGCTCATGCTTGTAGGAGGGAAAGTACt                                                                     >  2:6114822/1‑151 (MQ=17)
                                                                         ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                     >  1:206770/1‑151 (MQ=255)
                                                                         ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                     <  2:38264/151‑1 (MQ=255)
                                                                         ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                     >  1:4708897/1‑151 (MQ=255)
                                                                         ggCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACt                                                                     >  2:1494323/1‑151 (MQ=255)
                                                                          gCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTc                                                                    <  1:5375652/151‑1 (MQ=255)
                                                                           cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCCATCATGGTTGTAGTAGTGAATGTACTCa                                                                   >  1:4527487/1‑151 (MQ=255)
                                                                             cTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATc                                                                 >  2:688179/1‑151 (MQ=255)
                                                                               tCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATccc                                                               >  2:4669095/1‑150 (MQ=255)
                                                                                  cgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAc                                                            >  1:5046540/1‑151 (MQ=255)
                                                                                   gggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAcc                                                           <  1:2227806/151‑1 (MQ=255)
                                                                                    ggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAccc                                                          <  1:1505642/151‑1 (MQ=255)
                                                                                         gCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctct                                                     >  1:5581107/1‑151 (MQ=255)
                                                                                            cTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAg                                                  >  1:879511/1‑151 (MQ=255)
                                                                                             tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGc                                                 <  2:1558074/151‑1 (MQ=255)
                                                                                              cGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCt                                                <  1:3002622/151‑1 (MQ=255)
                                                                                              cGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCt                                                >  1:3170203/1‑151 (MQ=255)
                                                                                               gCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtc                                               <  1:2555547/151‑1 (MQ=255)
                                                                                               gCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtc                                               <  1:5238097/151‑1 (MQ=255)
                                                                                                cTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtca                                              >  1:5163396/1‑151 (MQ=255)
                                                                                                 tCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcat                                             <  2:5737307/151‑1 (MQ=255)
                                                                                                 tCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcat                                             >  2:616994/1‑151 (MQ=255)
                                                                                                  cccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcatc                                            >  1:1251005/1‑151 (MQ=255)
                                                                                                  cccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtcatc                                            >  1:1341117/1‑151 (MQ=255)
                                                                                                   ccACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTAATCCAGACCCGCTCTCAGCtcatca                                           >  2:3680164/1‑151 (MQ=255)
                                                                                                    cacaATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCaa                                          <  1:7159175/151‑1 (MQ=255)
                                                                                                      caATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATg                                        <  2:5905753/151‑1 (MQ=255)
                                                                                                       aaTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGc                                       >  1:3228236/1‑151 (MQ=255)
                                                                                                        aTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGct                                      >  1:2280509/1‑151 (MQ=19)
                                                                                                        aTTGGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGct                                      <  2:5144738/151‑1 (MQ=12)
                                                                                                           tGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCtctc                                   <  2:2488795/151‑1 (MQ=14)
                                                                                                           tGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCtctc                                   <  1:2874028/151‑1 (MQ=12)
                                                                                                           tGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCtctc                                   <  1:688179/151‑1 (MQ=12)
                                                                                                             gACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGa                                 <  1:5531558/151‑1 (MQ=14)
                                                                                                              aCTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGaa                                <  1:7029989/151‑1 (MQ=14)
                                                                                                               cTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAACGCTCTCGaaa                               >  2:2549537/1‑151 (MQ=12)
                                                                                                                 gCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACg                             >  1:5494587/1‑151 (MQ=12)
                                                                                                                   cccccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCt                           >  1:3546973/1‑151 (MQ=12)
                                                                                                                     cccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTc                         <  1:1033836/151‑1 (MQ=12)
                                                                                                                      ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCa                        <  2:4569581/151‑1 (MQ=9)
                                                                                                                      ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCa                        <  1:4105770/151‑1 (MQ=12)
                                                                                                                       cAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAg                       <  1:4021515/151‑1 (MQ=17)
                                                                                                                       cAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAg                       >  2:3843915/1‑151 (MQ=12)
                                                                                                                            gTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGa                  >  1:2809906/1‑151 (MQ=12)
                                                                                                                              tGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaa                <  1:435295/151‑1 (MQ=9)
                                                                                                                               ggACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaaa               >  2:4357117/1‑151 (MQ=12)
                                                                                                                                  cAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAAcc            >  2:2183658/1‑150 (MQ=12)
                                                                                                                                    gTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCa          >  1:5898932/1‑151 (MQ=14)
                                                                                                                                      tCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGa        <  2:2797280/151‑1 (MQ=14)
                                                                                                                                       cTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGAt       >  1:1788713/1‑151 (MQ=11)
                                                                                                                                            ctacAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGATTTGAg  <  1:4595426/148‑1 (MQ=1)
                                                                                                                                            ctacAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGATTTGAg  >  2:1722924/4‑151 (MQ=1)
                                                                                                                                            ctacAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGATTTGAg  <  2:880331/148‑1 (MQ=1)
                                                                                                                                                    |                                                                                                                                              
GTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGATTTGAG  >  CP009974/5513724‑5514014

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: