Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,906,096 T→A 5.3% K246N (AAA→AAT)  htpX ← putative endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,906,0960TA 5.3% 1070.2 / 34.3 311K246N (AAA→AAT) htpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (161/132);  new base A (13/4);  total (174/137)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 322 total aligned reads displayed.

TCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGC  >  CP009273/1905947‑1906245
                                                                                                                                                     |                                                                                                                                                     
 cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGttt                                                                                                                                                      >  2:2981061/1‑150 (MQ=255)
 cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGttt                                                                                                                                                      >  1:4694307/1‑150 (MQ=255)
  atatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGtttt                                                                                                                                                     >  1:2395887/1‑150 (MQ=255)
   tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTc                                                                                                                                                    <  1:3424125/150‑1 (MQ=255)
   tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTc                                                                                                                                                    <  1:1514512/150‑1 (MQ=255)
    atTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa                                                                                                                                                   >  1:3989772/1‑150 (MQ=255)
      tactacAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGCTTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGAGTTTCATAGCTGGTTTTCAgg                                                                                                                                                 >  1:1494336/1‑150 (MQ=255)
      tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAgg                                                                                                                                                 >  2:4082940/1‑150 (MQ=255)
        cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGgcg                                                                                                                                               <  2:1611924/150‑1 (MQ=255)
         ttCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGgcgc                                                                                                                                              >  1:4119128/1‑150 (MQ=255)
          tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCt                                                                                                                                             >  1:3384851/1‑150 (MQ=255)
            aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGc                                                                                                                                           >  2:1098441/1‑150 (MQ=255)
                atacACCCGTACGCAGAGCTTCAATTCGTTTATCCAGATGCGGGTGGTACATGAACAACTCACTGAGCGATTTCGACTAACCGATAATGCAGAGAGCCATCATGATGGTTGCTACTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAgag                                                                                                                                       >  1:2596631/4‑148 (MQ=255)
                   cACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCgg                                                                                                                                    >  2:2758337/1‑150 (MQ=255)
                    aCCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGc                                                                                                                                   >  1:1069204/1‑150 (MQ=255)
                    aCCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGc                                                                                                                                   >  1:2166936/1‑150 (MQ=255)
                    aCCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTGACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGc                                                                                                                                   <  2:1971985/150‑1 (MQ=255)
                       cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAAt                                                                                                                                >  2:881373/1‑150 (MQ=255)
                        gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATc                                                                                                                               >  1:2503521/1‑150 (MQ=255)
                        gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATc                                                                                                                               >  1:943480/1‑150 (MQ=255)
                        gTACGCAGAGCTTCAATACGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATc                                                                                                                               >  2:3413354/1‑150 (MQ=255)
                         tACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATCGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCa                                                                                                                              >  2:5187838/1‑150 (MQ=255)
                         tACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCa                                                                                                                              <  1:4451675/150‑1 (MQ=255)
                          aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACGGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAt                                                                                                                             <  2:278476/150‑1 (MQ=255)
                            gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAttt                                                                                                                           >  1:3626779/1‑150 (MQ=255)
                              agagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtc                                                                                                                         >  1:299483/1‑150 (MQ=255)
                              agagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtc                                                                                                                         >  1:797507/1‑150 (MQ=255)
                               gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCTGGCGCTGCAGCGCGGCAATCAttttat                                                                                                                        >  2:2468708/1‑148 (MQ=255)
                                  cTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcgc                                                                                                                     <  2:2651211/150‑1 (MQ=255)
                                   ttCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcgcg                                                                                                                    >  1:4189523/1‑150 (MQ=255)
                                   ttCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcgcg                                                                                                                    >  2:3551102/1‑150 (MQ=255)
                                       aTTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGacca                                                                                                                >  2:5071811/1‑150 (MQ=255)
                                       aTTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGGTTTCAGGCGCTGCAGCGCGGCAATCAGTTTCTCGCGaccc                                                                                                                >  1:2825499/1‑149 (MQ=255)
                                        gtCGTTTATCGAGCGGCGGGTGGGTCATGTTAAACTCACTGAGCGATTTCGTCTTACCGTTTTTGAAGTGAGCGATCATGCTGGTTGCTTCTTGCGGTTGATTGCTGGTTTTCAGGCGCTGCTGGGCGGCATTCTTTTTGTCgagaagta                                                                                                               <  2:4666602/149‑5 (MQ=255)
                                           gTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCACTCATTTTCTCGCGaccaacca                                                                                                            >  2:5112751/1‑150 (MQ=255)
                                              tATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGtt                                                                                                         >  1:279663/1‑150 (MQ=255)
                                               aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                        <  2:305074/150‑1 (MQ=255)
                                               aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                        >  2:1531844/1‑150 (MQ=255)
                                               aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                        <  2:1129215/150‑1 (MQ=255)
                                               aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                        <  1:5001849/150‑1 (MQ=255)
                                                  cAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGc                                                                                                     >  2:665541/1‑150 (MQ=255)
                                                     cggcggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGa                                                                                                  <  2:393361/150‑1 (MQ=255)
                                                        cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAAcc                                                                                               >  1:1259597/1‑150 (MQ=255)
                                                            tgggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCa                                                                                           >  1:1239437/1‑150 (MQ=255)
                                                             gggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTGTTGCCGAACCGGCAt                                                                                          >  1:1634457/1‑150 (MQ=255)
                                                              ggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATc                                                                                         <  1:3715302/150‑1 (MQ=255)
                                                                  ctGAACCACTCCCTGAGCGATTTCGACTGACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCa                                                                                     <  2:413522/149‑1 (MQ=255)
                                                                  aTGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCa                                                                                     <  2:2886423/150‑1 (MQ=255)
                                                                  aTGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCa                                                                                     >  2:2148000/1‑150 (MQ=255)
                                                                  aTGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCa                                                                                     >  1:2343333/1‑150 (MQ=255)
                                                                   tGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCAt                                                                                    <  1:71144/150‑1 (MQ=255)
                                                                   tGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCAt                                                                                    >  1:2925816/1‑150 (MQ=255)
                                                                    gAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATg                                                                                   >  2:5011570/1‑150 (MQ=255)
                                                                      acaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGa                                                                                 <  1:358312/150‑1 (MQ=255)
                                                                      acaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGa                                                                                 <  2:1145118/150‑1 (MQ=255)
                                                                        aacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAAt                                                                               <  1:3846491/150‑1 (MQ=255)
                                                                           tcccTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCa                                                                            <  2:4334630/147‑1 (MQ=255)
                                                                             aCTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACg                                                                          <  1:2539824/150‑1 (MQ=255)
                                                                              cTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGa                                                                         >  2:3608009/1‑150 (MQ=255)
                                                                              cTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGa                                                                         >  1:1143766/1‑150 (MQ=255)
                                                                              cTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGa                                                                         >  2:3330368/1‑150 (MQ=255)
                                                                               tGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATCGCTGGATTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGAt                                                                        <  1:2077905/150‑1 (MQ=255)
                                                                                gAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATg                                                                       <  1:4264689/150‑1 (MQ=255)
                                                                                gAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATg                                                                       >  1:1736616/1‑150 (MQ=255)
                                                                                    gATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGAcgc                                                                   >  1:755815/1‑150 (MQ=255)
                                                                                    gATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGATTTGCCGAACCGGCATCAGCATGGAATACACGATTAcgc                                                                   >  2:4613399/1‑150 (MQ=255)
                                                                                     aTTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGAcgcg                                                                  <  1:4567461/150‑1 (MQ=255)
                                                                                      tttCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCga                                                                 >  2:389796/1‑150 (MQ=255)
                                                                                        tCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCgaga                                                               >  2:4118846/1‑150 (MQ=255)
                                                                                          gACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAc                                                             >  2:4018984/1‑150 (MQ=255)
                                                                                           aCTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAcc                                                            <  2:666398/150‑1 (MQ=255)
                                                                                           aCTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAcc                                                            >  1:4772883/1‑150 (MQ=255)
                                                                                              tACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACcaca                                                         <  1:5067682/150‑1 (MQ=255)
                                                                                               aCCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACAt                                                        <  2:3268709/150‑1 (MQ=255)
                                                                                               aCCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACAt                                                        >  1:5061149/1‑150 (MQ=255)
                                                                                                 cGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATgg                                                      <  1:1366153/150‑1 (MQ=255)
                                                                                                  gTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGt                                                     <  2:2231581/150‑1 (MQ=255)
                                                                                                       tGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATaa                                                <  1:575282/150‑1 (MQ=255)
                                                                                                        tcAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAAt                                               >  1:4907909/2‑150 (MQ=255)
                                                                                                             gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACAAGGTGATAATGCTCg                                          >  2:2208297/1‑150 (MQ=255)
                                                                                                               gCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGcc                                        <  1:1977372/150‑1 (MQ=255)
                                                                                                               gCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGcc                                        <  1:4582288/150‑1 (MQ=255)
                                                                                                                 catcatGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAg                                      >  2:2013251/1‑150 (MQ=255)
                                                                                                                   tcatGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGaa                                    <  2:3593014/150‑1 (MQ=255)
                                                                                                                    catGCTGGTTGCTTCTTGCGGTTCATAGCTTTTTTTCGGGAGTTTCAGAGCGGCTATCTTTTTCTCGCGTCCATCCATTTTTGCAGAAACGGCATCTGCATTTAATTCACTATGACGCGAGAACCACTTGGTGATAATGCTCGCCAGTAt                                   <  2:2137544/150‑1 (MQ=255)
                                                                                                                    catGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                   >  2:2401842/1‑150 (MQ=255)
                                                                                                                        cTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCa                               >  2:4700457/1‑150 (MQ=255)
                                                                                                                         tGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCaa                              >  2:4279351/1‑150 (MQ=255)
                                                                                                                          ggTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCaaa                             >  2:3074891/1‑150 (MQ=255)
                                                                                                                            tagCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAaca                           <  1:3042626/148‑1 (MQ=255)
                                                                                                                              gCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACAcc                         <  1:1098441/150‑1 (MQ=255)
                                                                                                                              gCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAAACCAAACAcc                         >  1:3487233/1‑150 (MQ=255)
                                                                                                                               cttcttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCa                        >  1:4384563/1‑150 (MQ=255)
                                                                                                                                   ttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTc                    <  1:3150773/150‑1 (MQ=255)
                                                                                                                                      cGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAg                 >  2:3693823/1‑150 (MQ=255)
                                                                                                                                           cATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCg            <  2:1750058/150‑1 (MQ=255)
                                                                                                                                             tAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCCGAATACCAAACACCAGTTCCAGAACCGtt          >  2:1420019/1‑150 (MQ=255)
                                                                                                                                                 tGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGATTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCATAACCGTTGCaa      >  1:3247014/1‑150 (MQ=255)
                                                                                                                                                   gTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAAcc    <  2:4872079/150‑1 (MQ=255)
                                                                                                                                                   gTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCGCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAAcc    <  2:4820004/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCGC  >  CP009273/1905947‑1906245

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