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.6% 800.1 / NA 236A266V (GCA→GTA) RPPX_24600integrase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (104/131);  minor base C (1/0);  total (105/131)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

GTCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTGCAGCCTGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGAT  >  CP009974/5513718‑5514009
                                                                                                                                                    |                                                                                                                                               
 tCCTGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAg                                                                                                                                              >  2:5234495/1‑151 (MQ=255)
    tGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGccc                                                                                                                                           >  1:4196836/1‑150 (MQ=255)
    tGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGccc                                                                                                                                           <  1:4294090/151‑2 (MQ=255)
    tGGTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGccc                                                                                                                                           >  1:971945/1‑150 (MQ=255)
     ggTTCGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTCAGCTCCAccccc                                                                                                                                          >  2:3175153/1‑149 (MQ=255)
      gTTCGAGACTTATGGCTCAGGCCGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTCCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACCGTTCTAAGCTACAGcccgc                                                                                                                                         >  2:4448943/1‑148 (MQ=255)
         cGAGACTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGcc                                                                                                                                      <  1:1978829/151‑1 (MQ=255)
          gagaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGccc                                                                                                                                     >  1:2237935/1‑151 (MQ=255)
          gagaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGccc                                                                                                                                     <  1:4716585/151‑1 (MQ=255)
           agaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCt                                                                                                                                    <  2:4478875/151‑1 (MQ=255)
           agaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCt                                                                                                                                    >  1:5187005/1‑151 (MQ=255)
           agaCTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCt                                                                                                                                    <  1:4095575/151‑1 (MQ=255)
             agttATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTCAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGt                                                                                                                                  <  2:5117940/149‑1 (MQ=255)
              cgtATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTa                                                                                                                                 <  2:4522171/149‑1 (MQ=255)
              cTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGTCTTGCGGCCCTTCGCGGGCTCGCCCCGGCTCCCACAGATTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTgca                                                                                                                                 >  2:96482/1‑149 (MQ=37)
              cTTATGGCTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACATCCCGAGCCCTGTa                                                                                                                                 >  1:1001979/1‑151 (MQ=255)
                     cTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGGACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCcaca                                                                                                                          <  1:82592/151‑1 (MQ=255)
                     cTCAGGCAGGAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCcaca                                                                                                                          <  2:633476/151‑1 (MQ=255)
                        aggcaggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAggg                                                                                                                       <  1:4928806/151‑1 (MQ=255)
                           caggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTc                                                                                                                    <  1:2778886/151‑1 (MQ=255)
                            aggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCGAAAGTTGGACAGTTCGAAGCTACACCCCGAGCCCTGTATTCCACAGGGCTcc                                                                                                                   >  2:4612049/1‑150 (MQ=37)
                            aggAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCa                                                                                                                   >  2:860986/1‑151 (MQ=255)
                             ggggTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAg                                                                                                                  <  2:3861245/148‑1 (MQ=255)
                              gAGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAgg                                                                                                                 <  2:4364329/151‑1 (MQ=255)
                               aGTGGCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGc                                                                                                                >  2:1488412/1‑151 (MQ=255)
                                   gCAGGCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCAtt                                                                                                            <  1:3275896/151‑1 (MQ=255)
                                      ggCGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAg                                                                                                         <  2:1319606/151‑1 (MQ=255)
                                        cGGGTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtc                                                                                                       >  1:2551535/1‑151 (MQ=255)
                                         gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtct                                                                                                      <  1:3838529/151‑1 (MQ=255)
                                         gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtct                                                                                                      <  2:2299502/151‑1 (MQ=255)
                                         gggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGACATTGAGtct                                                                                                      >  2:2852113/1‑151 (MQ=255)
                                          ggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtctc                                                                                                     <  1:2852401/151‑1 (MQ=255)
                                          ggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtctc                                                                                                     >  1:3699024/1‑151 (MQ=255)
                                          ggTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGtctc                                                                                                     <  2:2517407/151‑1 (MQ=255)
                                           gTATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTGCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCa                                                                                                    <  2:1008195/151‑1 (MQ=255)
                                            tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                   <  2:2183105/151‑1 (MQ=255)
                                            tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                   >  1:482727/1‑151 (MQ=255)
                                            tATCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAg                                                                                                   <  2:3633794/151‑1 (MQ=255)
                                             aTCAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGc                                                                                                  <  2:273297/151‑1 (MQ=255)
                                               cAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                                >  1:3584330/1‑151 (MQ=255)
                                               cAATACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTACTCCACAGGGCTCAGGCCATTGAGTCTCAGCtt                                                                                                >  1:554152/1‑151 (MQ=255)
                                                 aTACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGa                                                                                              >  1:3574071/1‑151 (MQ=255)
                                                  tACGATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGAt                                                                                             <  1:3513250/151‑1 (MQ=255)
                                                   aCGATTGAATGACCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATg                                                                                            <  1:2390597/151‑1 (MQ=255)
                                                    cGATTGAATGCACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACCATTTGGACTGCCCCCCAAAAGTTGGACCGTTCTCCGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCATGATGc                                                                                           <  2:1694238/151‑1 (MQ=18)
                                                     gATTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCg                                                                                          >  2:322722/1‑151 (MQ=255)
                                                      aTTGAATGAACCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGa                                                                                         <  1:4496326/151‑1 (MQ=255)
                                                        tgaatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATc                                                                                       <  2:2116080/151‑1 (MQ=255)
                                                         gaatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCGCACAATTTGGGCTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTCAGTCTCAGCTTGATGCGATCa                                                                                      >  2:5128801/1‑151 (MQ=18)
                                                          aatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                     <  2:2395472/151‑1 (MQ=255)
                                                          aatgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                     >  2:1804174/1‑151 (MQ=255)
                                                            tgaaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATgg                                                                                   >  1:4325490/1‑151 (MQ=255)
                                                              ggcgggggCCGCCTTGCGGCGCGTGGCGGGCTCGCCCTGGGTCCCACAATTTGGACTGCCCCCCAAGAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGtt                                                                                 <  1:4614110/147‑1 (MQ=17)
                                                              aaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGCGCCCTGTATTCCACAGGGCTCAGGCCCTTGAGTCTCAGCTTGATGCCATCACGGtg                                                                                 >  2:3503171/1‑150 (MQ=18)
                                                              aaCCGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGtt                                                                                 <  2:1083469/151‑1 (MQ=255)
                                                                    ggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCAt                                                                                     <  2:4155599/141‑1 (MQ=255)
                                                                      ccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTagtagt                                                                         >  1:796541/1‑151 (MQ=255)
                                                                      ccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTagtagt                                                                         <  1:9873/151‑1 (MQ=255)
                                                                      ccgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTagtagt                                                                         <  1:2061098/151‑1 (MQ=255)
                                                                       cgccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTg                                                                        >  1:4673938/1‑151 (MQ=255)
                                                                        gccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGGCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGa                                                                       >  1:3259120/1‑151 (MQ=255)
                                                                        gccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGa                                                                       <  2:3261202/151‑1 (MQ=255)
                                                                        gccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGa                                                                       <  2:691801/151‑1 (MQ=255)
                                                                         ccTTGCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGaa                                                                      <  1:4051650/151‑1 (MQ=255)
                                                                             gCGGCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTa                                                                  <  1:2968579/151‑1 (MQ=255)
                                                                                gCCCTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTc                                                               >  1:1536457/1‑151 (MQ=255)
                                                                                 cccTTCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCa                                                              <  1:1058339/151‑1 (MQ=255)
                                                                                 cccTTCGCGCGCTCGCCCTCGCTCCCACAATCTGGACCGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCTCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCa                                                              <  2:4106102/151‑1 (MQ=17)
                                                                                    ttCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATcc                                                           <  1:1306343/151‑1 (MQ=255)
                                                                                    ttCGCGGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTAATAATcc                                                           >  1:2373136/1‑151 (MQ=255)
                                                                                      cgcgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAg                                                         >  1:3053839/1‑151 (MQ=255)
                                                                                        cgGGCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAc                                                       <  1:4513613/151‑1 (MQ=255)
                                                                                         gggCTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGAcc                                                      >  1:3412869/1‑151 (MQ=255)
                                                                                            cTCGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGc                                                   >  2:4018397/1‑151 (MQ=255)
                                                                                              cGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctc                                                 <  1:3066639/151‑1 (MQ=255)
                                                                                              cGCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTCGGACAGTTCTAAGCTACATCCCGAGCCCCGTATCCCACAGGGCTCAGGCCATCGCGTCCCAGCTTGATACGACCATGGCTCTAGCAGTCAATGTACCCCTCCAGATCCGgtc                                                 >  2:3585272/1‑151 (MQ=255)
                                                                                               gCCCTCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGctct                                                >  2:590/1‑151 (MQ=255)
                                                                                                   tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTCGTAGTGAATGTACTCATCCAGACCCGCTCTCAGc                                            >  1:3926658/1‑151 (MQ=255)
                                                                                                   tCGCTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGc                                            >  1:211291/1‑151 (MQ=255)
                                                                                                      cTCCCACAATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCtca                                         >  2:3772177/1‑151 (MQ=255)
                                                                                                           acaATTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAAt                                    <  2:25700/151‑1 (MQ=255)
                                                                                                              aTTTGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGct                                 >  1:1115176/1‑151 (MQ=19)
                                                                                                                 tGGACTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCtctc                              >  1:70582/1‑151 (MQ=12)
                                                                                                                    aCTGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGaa                           >  1:2997395/1‑151 (MQ=12)
                                                                                                                      tGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                         <  2:3460961/151‑1 (MQ=14)
                                                                                                                      tGCCCCCCAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAAc                         >  1:4439711/1‑151 (MQ=12)
                                                                                                                          ccccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCtt                     <  1:1684171/151‑1 (MQ=14)
                                                                                                                          ccccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCtt                     <  1:504825/151‑1 (MQ=14)
                                                                                                                            ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCa                   <  2:3098895/151‑1 (MQ=14)
                                                                                                                            ccAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCa                   <  1:3130155/151‑1 (MQ=12)
                                                                                                                             cAAAAGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAg                  >  2:1360277/1‑151 (MQ=9)
                                                                                                                              aaaaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAgg                 <  2:972354/151‑1 (MQ=12)
                                                                                                                              aaaaGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAgg                 >  2:5044853/1‑151 (MQ=12)
                                                                                                                                 aGTTGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAg              <  2:806465/151‑1 (MQ=12)
                                                                                                                                    tGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaa           >  1:92647/1‑151 (MQ=12)
                                                                                                                                    tGGACAGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGaaa           >  1:4156529/1‑151 (MQ=12)
                                                                                                                                         aGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTc      >  2:2930388/1‑151 (MQ=9)
                                                                                                                                         aGTTCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTc      >  1:2284209/1‑151 (MQ=9)
                                                                                                                                         aGTTCTAACCTATAGCCCGTGCCCTGTAGTCGACGGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATGATGGTTGAAGTAGTGAATGTACTCACCCAGACCCGCTCTCAGCTCATCGATGCTCTCGAAACGCTTCAGGTAGAAAAACTc      <  2:4691171/151‑1 (MQ=11)
                                                                                                                                            tCTAAGCTACAGCCCGAGCCCTGTATTCCACAGGGCTCAGGCCATTGAGTCTCAGCTTGATGCGATCATGGTTGTAGTAGTGAATGTACTCATCCAGACCCGCTCTCAGCTCATCAATGCTCTCGAAACGCTTCAGGTAGAAAAACTCAGa   <  2:4306996/151‑1 (MQ=14)
                                                                                                                                                    |                                                                                                                                               
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: