Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,492 A→T 18.0% E63V (GAA→GTA)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,325,4920AT18.0% 550.2 / 108.8 216E63V (GAA→GTA) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base A (92/85);  new base T (22/17);  total (114/102)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 223 total aligned reads displayed.

CAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCT  >  CP009273/1325346‑1325639
                                                                                                                                                  |                                                                                                                                                   
cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGaa                                                                                                                                                    <  1:3177670/148‑1 (MQ=255)
cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGaa                                                                                                                                                    <  2:1851012/148‑1 (MQ=255)
cAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCAAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGaa                                                                                                                                                    <  1:722428/148‑1 (MQ=255)
     cGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTTCCGACTCTACCTCCACCAAGACGGCTAAAAAGACTAAAAAGGATGAACGTgg                                                                                                                                               >  2:563887/1‑148 (MQ=255)
     cGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTgg                                                                                                                                               >  2:2937080/1‑148 (MQ=255)
       aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGgcg                                                                                                                                             >  1:3325236/1‑148 (MQ=255)
         caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGct                                                                                                                                           <  2:169058/148‑1 (MQ=255)
          aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGctc                                                                                                                                          <  2:2069629/148‑1 (MQ=255)
               gTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTc                                                                                                                                     >  2:3228423/1‑148 (MQ=255)
                   cTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAAcc                                                                                                                                 >  1:528156/1‑148 (MQ=255)
                        tAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATg                                                                                                                            <  2:1623473/148‑1 (MQ=255)
                         aGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATgg                                                                                                                           >  1:2465209/1‑148 (MQ=255)
                            gACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATggggg                                                                                                                        <  1:2817408/148‑1 (MQ=255)
                             actactTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGt                                                                                                                       <  2:1215345/148‑1 (MQ=255)
                             aCTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGt                                                                                                                       <  2:590027/148‑1 (MQ=255)
                                 cttggtgAAATCCTGCGTCGGTCACATCCGCTTTTTGCCTACCAGTGGCTCAGCTGCAAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGAc                                                                                                                   <  2:821800/146‑1 (MQ=255)
                                      tgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTg                                                                                                              >  1:1156132/1‑148 (MQ=255)
                                         aaTCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGca                                                                                                           <  1:71196/148‑1 (MQ=255)
                                           tCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGcaca                                                                                                         >  1:785506/1‑148 (MQ=255)
                                            ccAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACaa                                                                                                        <  1:1623415/148‑1 (MQ=255)
                                             cAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAAt                                                                                                       <  1:606680/148‑1 (MQ=255)
                                             cAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAAt                                                                                                       >  2:3240050/1‑148 (MQ=255)
                                              aGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCCCAAtt                                                                                                      >  1:1476358/1‑148 (MQ=255)
                                               gCGTCGGTCACATCCGCGATTTGCCGACCCGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTg                                                                                                     <  2:314999/148‑1 (MQ=255)
                                               gCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTg                                                                                                     <  2:3379184/148‑1 (MQ=255)
                                                   cGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAg                                                                                                 <  2:2921049/148‑1 (MQ=255)
                                                      tCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGgcg                                                                                              >  2:522245/1‑148 (MQ=255)
                                                        acaTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCa                                                                                            <  1:474274/148‑1 (MQ=255)
                                                           tCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTa                                                                                         <  2:1653851/148‑1 (MQ=255)
                                                           tCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTa                                                                                         <  2:3015027/148‑1 (MQ=255)
                                                           tCCGCGATTTGCCGACCAGGGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTa                                                                                         <  2:2423642/148‑1 (MQ=255)
                                                            ccGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTAt                                                                                        <  1:2716702/148‑1 (MQ=255)
                                                              gcgATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGa                                                                                      <  1:2309722/148‑1 (MQ=255)
                                                                 aTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgt                                                                                   <  2:2852111/148‑1 (MQ=255)
                                                                 aTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgt                                                                                   >  1:3383676/1‑148 (MQ=255)
                                                                  tttGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgtg                                                                                  >  2:1063705/1‑148 (MQ=255)
                                                                    tGCCGCCCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGtt                                                                                <  2:2597644/148‑1 (MQ=255)
                                                                     gCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTg                                                                               <  1:1838064/148‑1 (MQ=255)
                                                                         aCCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTg                                                                           >  1:1004182/1‑148 (MQ=255)
                                                                          ccAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTgg                                                                          >  1:345791/1‑148 (MQ=255)
                                                                           cAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGt                                                                         <  1:138866/148‑1 (MQ=255)
                                                                                gCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAaga                                                                    >  2:2715308/1‑148 (MQ=255)
                                                                                gCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAaga                                                                    <  2:2855677/148‑1 (MQ=255)
                                                                                 cTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAagag                                                                   <  2:3034779/148‑1 (MQ=255)
                                                                                  tCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAgaga                                                                  <  1:3493136/148‑1 (MQ=255)
                                                                                     gCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAgg                                                               <  2:3265100/148‑1 (MQ=255)
                                                                                      cTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgt                                                              >  2:2307482/1‑148 (MQ=255)
                                                                                        gCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcg                                                            >  1:2956228/1‑148 (MQ=255)
                                                                                         ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                           <  2:2471337/148‑1 (MQ=255)
                                                                                           aaaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGtct                                                         <  2:1327449/148‑1 (MQ=255)
                                                                                           aaaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGtct                                                         <  2:2847078/148‑1 (MQ=255)
                                                                                            aaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGtctc                                                        >  1:3540462/1‑148 (MQ=255)
                                                                                             aaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTctct                                                       <  2:1986472/148‑1 (MQ=255)
                                                                                              aaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctg                                                      >  2:1252756/1‑148 (MQ=255)
                                                                                              aaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctg                                                      >  2:2468216/1‑148 (MQ=255)
                                                                                              aaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctg                                                      >  2:2521031/1‑148 (MQ=255)
                                                                                                gagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaa                                                    <  2:1439386/148‑1 (MQ=255)
                                                                                                gagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaa                                                    >  1:3658144/1‑148 (MQ=255)
                                                                                                gagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaa                                                    >  2:1492165/1‑148 (MQ=255)
                                                                                                gagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaa                                                    >  2:555116/1‑148 (MQ=255)
                                                                                                 agTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaac                                                   <  2:3598282/148‑1 (MQ=255)
                                                                                                 agTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaac                                                   >  1:384472/1‑148 (MQ=255)
                                                                                                 agTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaac                                                   <  1:2152710/148‑1 (MQ=255)
                                                                                                 agTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaac                                                   <  1:990803/148‑1 (MQ=255)
                                                                                                  gTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaact                                                  >  1:3638544/1‑148 (MQ=255)
                                                                                                  gTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaact                                                  >  1:1745048/1‑148 (MQ=255)
                                                                                                  gTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaact                                                  >  1:1770574/1‑148 (MQ=255)
                                                                                                   tGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTCTGGGGGTTGACCCGTGGCACAAGTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactg                                                 >  1:3581254/1‑148 (MQ=255)
                                                                                                        aCTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaaca                                            <  1:563786/148‑1 (MQ=255)
                                                                                                        aCTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaaca                                            <  1:3266662/148‑1 (MQ=255)
                                                                                                          tctACCTCCACCAAGCCGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaac                                          <  1:2250323/148‑1 (MQ=255)
                                                                                                           ctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACt                                         <  2:1156132/148‑1 (MQ=255)
                                                                                                            tACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTg                                        <  1:3793906/148‑1 (MQ=255)
                                                                                                             aCCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTgg                                       <  2:2159215/148‑1 (MQ=255)
                                                                                                             aCCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTgg                                       <  2:754362/148‑1 (MQ=255)
                                                                                                             aCCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTgg                                       >  1:811812/1‑148 (MQ=255)
                                                                                                                 ccaccaAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGa                                   >  2:1417387/1‑148 (MQ=255)
                                                                                                                    ccaAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaaa                                >  2:73293/1‑148 (MQ=255)
                                                                                                                    ccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaaa                                >  1:2324566/1‑148 (MQ=255)
                                                                                                                    ccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaaa                                <  1:1201462/148‑1 (MQ=255)
                                                                                                                       aGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGc                             >  2:137348/1‑148 (MQ=255)
                                                                                                                        gACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGcc                            >  2:1084589/1‑148 (MQ=255)
                                                                                                                         aCGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCg                           >  1:2119759/1‑148 (MQ=255)
                                                                                                                         aCGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCg                           <  2:384463/148‑1 (MQ=255)
                                                                                                                             cTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACca                       >  2:846932/1‑148 (MQ=255)
                                                                                                                               aaaaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACcaca                     >  1:1578932/1‑148 (MQ=255)
                                                                                                                                aaaaGCCTAAAAAGGATTAACGTGGCGCTCTCGTCAACCGTAAGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACAGGCTGAAAAAGCCGACCACat                    >  1:1304302/1‑148 (MQ=255)
                                                                                                                                 aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                   >  1:92907/1‑148 (MQ=255)
                                                                                                                                 aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                   >  1:3294724/1‑148 (MQ=255)
                                                                                                                                  aaGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatct                  <  2:3456407/148‑1 (MQ=255)
                                                                                                                                  aaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGGGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGACCTGAAACAACTGGCTGAAAAAGCCGACCACatct                  <  1:3240050/148‑1 (MQ=255)
                                                                                                                                   aGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatcta                 >  1:1381994/1‑148 (MQ=255)
                                                                                                                                   aGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatcta                 >  1:2352587/1‑148 (MQ=255)
                                                                                                                                    gCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctat                >  2:2097622/1‑148 (MQ=255)
                                                                                                                                       tAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTAtctc             >  1:3044328/1‑148 (MQ=255)
                                                                                                                                          aaaGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCa          <  1:1067152/148‑1 (MQ=255)
                                                                                                                                            aGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAAc        >  2:1456494/1‑148 (MQ=255)
                                                                                                                                              gATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCg      >  2:360681/1‑148 (MQ=255)
                                                                                                                                                 gAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGAcc   <  1:1402426/148‑1 (MQ=255)
                                                                                                                                                 gAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATATATCTCGCAACCGAcc   <  1:2268594/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCT  >  CP009273/1325346‑1325639

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: