Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,492 A→T 27.9% 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,4920AT27.9% 394.0 / 191.8 233E63V (GAA→GTA) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base A (96/71);  new base T (33/32);  total (129/104)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 235 total aligned reads displayed.

CCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCT  >  CP009273/1325345‑1325639
                                                                                                                                                   |                                                                                                                                                   
  aaaaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAAc                                                                                                                                                   >  2:2670764/1‑148 (MQ=255)
      cGATGAACAATTATCTGGTTAGTGAATACGTGGTGAAATCCAGCGTCTGTCACATCCGCGTTTTGCCGACCAGTGTCTCTGATGCCAAAAATTGTGCCGACTCTAACTCCACCAAGACGGATAAAAAGCCTAAAAAGGATGTACGTgg                                                                                                                                               <  2:1630177/148‑1 (MQ=255)
        aTCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGgcg                                                                                                                                             <  2:2250009/148‑1 (MQ=255)
              aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCg                                                                                                                                       >  1:1538577/1‑148 (MQ=255)
              aaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCg                                                                                                                                       <  2:2888703/148‑1 (MQ=255)
                  aTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCaa                                                                                                                                   >  1:2069074/1‑148 (MQ=255)
                   tCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAAc                                                                                                                                  >  1:2613139/1‑148 (MQ=255)
                      gggTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGt                                                                                                                               <  2:1979466/148‑1 (MQ=255)
                       ggTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTa                                                                                                                              >  2:2229645/1‑148 (MQ=255)
                         tAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATg                                                                                                                            >  1:687109/1‑148 (MQ=255)
                         tAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATg                                                                                                                            >  1:242149/1‑148 (MQ=255)
                          aGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATgg                                                                                                                           >  1:1435648/1‑148 (MQ=255)
                          aGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATgg                                                                                                                           >  2:2809743/1‑148 (MQ=255)
                            tGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATgggg                                                                                                                         >  2:2566069/1‑148 (MQ=255)
                              aCTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGt                                                                                                                       >  2:1000117/1‑148 (MQ=255)
                                tACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTg                                                                                                                     <  2:504910/148‑1 (MQ=255)
                                  cGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGAc                                                                                                                   >  2:457191/1‑148 (MQ=255)
                                  cGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGAc                                                                                                                   >  2:1506844/1‑148 (MQ=255)
                                     ggtgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCCCCAAGACGGCTAAAAAGCCGAAAAAGGATGAACGTGGCGCTCTCGTCCACCGGATGGGGGTTGACCCg                                                                                                                >  1:2221933/1‑148 (MQ=255)
                                     ggtgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCg                                                                                                                <  1:736774/148‑1 (MQ=255)
                                       tgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGCTGAACGTGGCGCTCTCGTCAACCGTATGTGGGTTTACCCGTg                                                                                                              >  2:3021897/1‑148 (MQ=255)
                                          aaTCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGca                                                                                                           >  1:341206/1‑148 (MQ=255)
                                             ccAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAACAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACaa                                                                                                        <  2:1086673/148‑1 (MQ=255)
                                                  gTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTggg                                                                                                   <  2:3356487/148‑1 (MQ=255)
                                                   tCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGa                                                                                                  >  1:2719571/1‑148 (MQ=255)
                                                     ggTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAgg                                                                                                >  2:712878/1‑148 (MQ=255)
                                                       tCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGgcg                                                                                              >  1:32911/1‑148 (MQ=255)
                                                        cacaTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGgcgc                                                                                             >  1:3030247/1‑148 (MQ=255)
                                                          caTCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCAc                                                                                           >  2:1609914/1‑148 (MQ=255)
                                                            tCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTa                                                                                         <  2:2582402/148‑1 (MQ=255)
                                                            tCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTa                                                                                         <  2:544234/148‑1 (MQ=255)
                                                             ccGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTAt                                                                                        >  2:53150/1‑148 (MQ=255)
                                                                   tttGCCGACCAGTGGCTAAGCTGACAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgtg                                                                                  <  2:2455062/148‑1 (MQ=255)
                                                                    ttGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgtgt                                                                                 >  1:3311396/1‑148 (MQ=255)
                                                                    ttGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgtgt                                                                                 >  1:1513360/1‑148 (MQ=255)
                                                                      gCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTg                                                                               <  2:2872115/148‑1 (MQ=255)
                                                                       ccGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGc                                                                              >  1:2004097/1‑148 (MQ=255)
                                                                       ccGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGc                                                                              >  1:2002349/1‑148 (MQ=255)
                                                                       ccGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGc                                                                              <  1:678611/148‑1 (MQ=255)
                                                                        cGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGcc                                                                             <  2:927948/148‑1 (MQ=255)
                                                                         gACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCt                                                                            >  2:3525845/1‑148 (MQ=255)
                                                                          aCCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTg                                                                           >  1:29546/1‑148 (MQ=255)
                                                                           ccAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTgg                                                                          <  2:2162737/148‑1 (MQ=255)
                                                                              gTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTaa                                                                       >  2:46088/1‑148 (MQ=255)
                                                                               tGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTaaa                                                                      >  2:1826016/1‑148 (MQ=255)
                                                                                ggCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAag                                                                     <  1:3227400/148‑1 (MQ=255)
                                                                                   taaGCTGCAAAAAAGAGTGCCGACTCTAACTCCACCAAGTCTGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTTTGGGGGTTGACCCGTGTCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAgaga                                                                  <  1:621073/146‑1 (MQ=255)
                                                                                      gCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAgg                                                               <  1:1673869/148‑1 (MQ=255)
                                                                                      gCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAgg                                                               <  1:28022/148‑1 (MQ=255)
                                                                                          cccaaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAGGAGAAGgtcgt                                                           >  2:3063745/4‑148 (MQ=255)
                                                                                          ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                           <  1:1004115/148‑1 (MQ=255)
                                                                                          ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                           >  2:3063924/1‑148 (MQ=255)
                                                                                          ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                           <  2:1559864/148‑1 (MQ=255)
                                                                                          ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                           <  2:115069/148‑1 (MQ=255)
                                                                                           cAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgtc                                                          >  1:1253629/1‑148 (MQ=255)
                                                                                            acacagagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGtct                                                         <  2:741806/144‑1 (MQ=255)
                                                                                                agagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctga                                                     <  2:2069159/148‑1 (MQ=255)
                                                                                                  agTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaac                                                   >  1:3309247/1‑148 (MQ=255)
                                                                                                   gTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaact                                                  >  2:164807/1‑148 (MQ=255)
                                                                                                       cGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGaaa                                              <  2:2756953/148‑1 (MQ=255)
                                                                                                       cGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGaaa                                              <  1:2237287/148‑1 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                           <  2:1243483/148‑1 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGTAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                           <  1:2986648/148‑1 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                           >  2:1286575/1‑148 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                           <  1:1890258/148‑1 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                           <  1:3270983/148‑1 (MQ=255)
                                                                                                           tctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaac                                          >  1:1255698/1‑148 (MQ=255)
                                                                                                           tctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaac                                          <  2:1401418/148‑1 (MQ=255)
                                                                                                            ctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACt                                         >  2:825666/1‑148 (MQ=255)
                                                                                                             tACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTg                                        >  2:2090367/1‑148 (MQ=255)
                                                                                                              aCCTCCACCAAGACGGCTAAAAAGCCTAATAAGGATGAACGTGGCGCTCTCGTCAACCGTTTGGGGTTTGACCCTTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAATGAGAAGGTCGTCTCTGAACTTAAACAACTgg                                       <  2:2985521/148‑1 (MQ=255)
                                                                                                              aCCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTgg                                       >  1:1369645/1‑148 (MQ=255)
                                                                                                                cTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCt                                     >  1:1252831/1‑148 (MQ=255)
                                                                                                                 tCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTg                                    >  2:1448399/1‑148 (MQ=255)
                                                                                                                 tCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTg                                    >  1:2123131/1‑148 (MQ=255)
                                                                                                                  ccaccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGATGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGa                                   <  1:1847430/148‑1 (MQ=255)
                                                                                                                    accaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaa                                 >  1:2241440/1‑148 (MQ=255)
                                                                                                                     ccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaaa                                >  1:1196920/1‑148 (MQ=255)
                                                                                                                        aGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGc                             >  2:2038386/1‑148 (MQ=255)
                                                                                                                         gACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGcc                            >  1:175473/1‑148 (MQ=255)
                                                                                                                         gACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGcc                            >  2:3262594/1‑148 (MQ=255)
                                                                                                                          aCGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCg                           >  1:455300/1‑148 (MQ=255)
                                                                                                                           cGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGa                          <  2:2373134/148‑1 (MQ=255)
                                                                                                                             gCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGAcc                        <  2:2613139/148‑1 (MQ=255)
                                                                                                                              cTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTCGCTGAAAAAGCCGACca                       >  2:1937782/1‑148 (MQ=255)
                                                                                                                                 aaaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACat                    >  1:953068/1‑148 (MQ=255)
                                                                                                                                  aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                   >  1:2148964/1‑148 (MQ=255)
                                                                                                                                  aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                   <  1:1826016/148‑1 (MQ=255)
                                                                                                                                    aGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatcta                 >  2:1264527/1‑148 (MQ=255)
                                                                                                                                      ccTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatc               <  1:1355347/148‑1 (MQ=255)
                                                                                                                                       cTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatct              <  2:341206/148‑1 (MQ=255)
                                                                                                                                           aaaGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCa          <  2:1538577/148‑1 (MQ=255)
                                                                                                                                           aaaGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCa          <  1:2363906/148‑1 (MQ=255)
                                                                                                                                            aaGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCaa         >  1:3520729/1‑148 (MQ=255)
                                                                                                                                             aGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAAc        >  1:495472/1‑148 (MQ=255)
                                                                                                                                               gATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATGGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCCCATCTATCTCGCACCCg      >  2:1056262/1‑148 (MQ=255)
                                                                                                                                               gATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCg      >  1:578001/1‑148 (MQ=255)
                                                                                                                                                aTGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGa     >  1:1058242/1‑148 (MQ=255)
                                                                                                                                                aTGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGa     <  2:865824/148‑1 (MQ=255)
                                                                                                                                                 tGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGAc    >  1:980491/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGACCT  >  CP009273/1325345‑1325639

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: