Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,906,096 T→A 6.5% 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 6.5% 723.9 / 30.9 211K246N (AAA→AAThtpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (91/105);  new base A (6/8);  total (97/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 215 total aligned reads displayed.

CATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCG  >  CP009273/1905948‑1906244
                                                                                                                                                    |                                                                                                                                                    
cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGttt                                                                                                                                                     >  1:3971434/1‑150 (MQ=255)
   atTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCa                                                                                                                                                  >  1:2068823/1‑150 (MQ=255)
     tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAgg                                                                                                                                                <  1:2761255/150‑1 (MQ=255)
     tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAgg                                                                                                                                                >  2:3983391/1‑150 (MQ=255)
         tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCt                                                                                                                                            <  1:1316033/150‑1 (MQ=255)
         tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCt                                                                                                                                            >  2:2236293/1‑150 (MQ=255)
            ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCa                                                                                                                                         >  1:920497/1‑150 (MQ=255)
                 tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGcgcg                                                                                                                                    >  1:4135282/1‑150 (MQ=255)
                  cACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCgg                                                                                                                                   <  2:1092377/150‑1 (MQ=255)
                  cACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCgg                                                                                                                                   >  1:2179971/1‑150 (MQ=255)
                      cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAAt                                                                                                                               >  2:639189/1‑150 (MQ=255)
                      cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAAt                                                                                                                               >  2:4934323/1‑150 (MQ=255)
                       gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATc                                                                                                                              >  2:26051/1‑150 (MQ=255)
                        tACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCa                                                                                                                             <  2:4271969/150‑1 (MQ=255)
                        tACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCa                                                                                                                             <  1:2971955/150‑1 (MQ=255)
                         aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAt                                                                                                                            <  2:3912062/150‑1 (MQ=255)
                         aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAt                                                                                                                            >  1:4845925/1‑150 (MQ=255)
                          cGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAtt                                                                                                                           >  2:1371730/1‑150 (MQ=255)
                           gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCAttt                                                                                                                          >  1:5144916/1‑150 (MQ=255)
                               agCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTtctc                                                                                                                      >  2:3548510/1‑150 (MQ=255)
                                gCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTcg                                                                                                                     >  2:2568856/1‑150 (MQ=255)
                                          gTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGaccaacca                                                                                                           >  2:4608254/1‑150 (MQ=255)
                                              aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                       <  2:1432519/150‑1 (MQ=255)
                                              aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGttt                                                                                                       <  2:895091/150‑1 (MQ=255)
                                               tCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGtttt                                                                                                      <  1:4808863/150‑1 (MQ=255)
                                                ccAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTg                                                                                                     <  2:2104183/150‑1 (MQ=255)
                                                ccAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTg                                                                                                     <  2:2962696/150‑1 (MQ=255)
                                                  aGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGcc                                                                                                   >  1:5148457/1‑150 (MQ=255)
                                                            gggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCAt                                                                                         >  1:2414350/1‑150 (MQ=255)
                                                             ggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATc                                                                                        <  2:1295100/150‑1 (MQ=255)
                                                              gtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCa                                                                                       <  1:1190465/150‑1 (MQ=255)
                                                                cATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGc                                                                                     <  1:5024773/150‑1 (MQ=255)
                                                                 aTGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCa                                                                                    >  1:2255275/1‑150 (MQ=255)
                                                                  tGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCAt                                                                                   <  2:153807/150‑1 (MQ=255)
                                                                   gAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATg                                                                                  <  2:5119307/150‑1 (MQ=255)
                                                                   gAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATg                                                                                  >  2:4573713/1‑150 (MQ=255)
                                                                   gAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATg                                                                                  >  2:560609/1‑150 (MQ=255)
                                                                    aacaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATgg                                                                                 >  2:2421970/1‑150 (MQ=255)
                                                                    aacaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATgg                                                                                 >  2:2144219/1‑150 (MQ=255)
                                                                     acaacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGa                                                                                <  1:3073991/150‑1 (MQ=255)
                                                                      caacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGaa                                                                               >  2:3099634/1‑150 (MQ=255)
                                                                       aacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAAt                                                                              <  1:3312653/150‑1 (MQ=255)
                                                                       aacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAAt                                                                              <  1:3304860/150‑1 (MQ=255)
                                                                       aacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAAt                                                                              <  2:5102300/150‑1 (MQ=255)
                                                                       aacTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAAt                                                                              <  1:699573/150‑1 (MQ=255)
                                                                          tCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCa                                                                           <  2:64974/150‑1 (MQ=255)
                                                                           cACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCAc                                                                          <  1:2177742/150‑1 (MQ=255)
                                                                           cACTGAGCGAGTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCAc                                                                          <  2:3727842/150‑1 (MQ=255)
                                                                            aCTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACg                                                                         <  1:4416357/150‑1 (MQ=255)
                                                                                  cGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCCATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGAcg                                                                   >  1:4381416/1‑150 (MQ=255)
                                                                                  cGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGAcg                                                                   >  1:2805271/1‑150 (MQ=255)
                                                                                         gACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAc                                                            >  2:2749733/1‑150 (MQ=255)
                                                                                         gACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAc                                                            >  1:4654084/1‑150 (MQ=255)
                                                                                         gACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAc                                                            >  1:3140430/1‑150 (MQ=255)
                                                                                         gACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAc                                                            >  1:3968091/1‑150 (MQ=255)
                                                                                          aCTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAAcc                                                           >  1:4614453/1‑150 (MQ=255)
                                                                                           cTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACca                                                          >  2:878067/1‑150 (MQ=255)
                                                                                            ttACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACcac                                                         <  2:352666/150‑1 (MQ=255)
                                                                                             tACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACcaca                                                        >  2:4035138/1‑150 (MQ=255)
                                                                                              aCCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACAt                                                       >  1:4914910/1‑150 (MQ=255)
                                                                                                  ttAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTg                                                   <  2:5151300/150‑1 (MQ=255)
                                                                                                   tAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGa                                                  <  1:4484633/150‑1 (MQ=255)
                                                                                                    aaTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGAt                                                 <  2:2475968/150‑1 (MQ=255)
                                                                                                     aTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATa                                                >  1:2538349/1‑150 (MQ=255)
                                                                                                     aTGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATa                                                >  2:4891332/1‑150 (MQ=255)
                                                                                                       gCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAAt                                              <  1:4158563/150‑1 (MQ=255)
                                                                                                       gCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAAt                                              <  1:3479889/150‑1 (MQ=255)
                                                                                                         agagagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGc                                            >  1:3676449/1‑150 (MQ=255)
                                                                                                         agagagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGc                                            >  1:4742650/1‑150 (MQ=255)
                                                                                                           agagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTc                                          <  1:2236293/150‑1 (MQ=255)
                                                                                                            gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCg                                         <  1:986176/150‑1 (MQ=255)
                                                                                                            gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCg                                         <  2:66842/150‑1 (MQ=255)
                                                                                                            gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCg                                         >  2:2833397/1‑150 (MQ=255)
                                                                                                            gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCg                                         >  2:3058323/1‑150 (MQ=255)
                                                                                                            gagCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCg                                         <  1:985566/150‑1 (MQ=255)
                                                                                                               ccATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCa                                      >  2:1831877/1‑150 (MQ=255)
                                                                                                                   catGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                  <  2:5050900/150‑1 (MQ=255)
                                                                                                                   catGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                  >  1:3769882/1‑150 (MQ=255)
                                                                                                                   catGCTGGTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAAt                                  <  1:3461958/150‑1 (MQ=255)
                                                                                                                      gCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATAcc                               <  1:4030857/150‑1 (MQ=255)
                                                                                                                       cTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCa                              <  2:4212494/150‑1 (MQ=255)
                                                                                                                         ggTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCaaa                            >  1:603736/1‑150 (MQ=255)
                                                                                                                          gTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAac                           >  2:485052/1‑150 (MQ=255)
                                                                                                                          gTTGCTTCTTGCGGTTCATAGCTGGTATTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAac                           >  2:4831858/1‑150 (MQ=255)
                                                                                                                              cttcttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCa                       <  1:3857751/150‑1 (MQ=255)
                                                                                                                                 cttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGtt                    <  1:1812333/150‑1 (MQ=255)
                                                                                                                                  ttGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTc                   <  1:952050/150‑1 (MQ=255)
                                                                                                                                   tGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTcc                  <  1:1019599/150‑1 (MQ=255)
                                                                                                                                     cGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAg                <  1:1089526/150‑1 (MQ=255)
                                                                                                                                     cGGTTCATAGCTGGTGTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAg                <  2:5008248/150‑1 (MQ=255)
                                                                                                                                      ggTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGa               >  1:1566706/1‑150 (MQ=255)
                                                                                                                                        ttCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAAc             <  1:5103067/150‑1 (MQ=255)
                                                                                                                                        ttCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAAc             >  1:1485684/1‑150 (MQ=255)
                                                                                                                                         tCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAAcc            <  2:5144916/150‑1 (MQ=255)
                                                                                                                                           aTAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGt          >  2:780418/1‑150 (MQ=255)
                                                                                                                                            tAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGtt         <  2:4811746/150‑1 (MQ=255)
                                                                                                                                            tAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGtt         <  2:4654023/150‑1 (MQ=255)
                                                                                                                                                 ggTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATACTGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAAc    >  2:2307372/1‑150 (MQ=255)
                                                                                                                                                  gTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAAcc   >  1:2962695/1‑150 (MQ=255)
                                                                                                                                                  gTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAAcc   <  2:2068823/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTGGTTTTCAGGCGCTGCAGCGCGGCAATCATTTTCTCGCGACCAACCAGTTTTGCCGAACCGGCATCAGCATGGAATTCACGATGACGCGAGAACCACATGGTGATAATGCTCGCCAGAATACCAAACACCAGTTCCAGAACCGTTGCAACCG  >  CP009273/1905948‑1906244

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