Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,479,733 G→T 100% D806Y (GAT→TAT)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,479,7330GT96.4% 688.9 / 14.6 224D806Y (GAT→TAT) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base T (95/121);  total (100/124)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 224 total aligned reads displayed.

GTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCATC  >  CP009273/1479587‑1479877
                                                                                                                                                  |                                                                                                                                                
gTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGa                                                                                                                                                 >  2:1872798/1‑148 (MQ=255)
    cTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTa                                                                                                                                             <  1:1186643/148‑1 (MQ=255)
    cTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTa                                                                                                                                             <  1:78937/148‑1 (MQ=255)
       gCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtc                                                                                                                                          >  1:3014751/1‑148 (MQ=255)
         aGACGCGTCACGCATTCGTCCGTGAAAACCTGAAACGACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGATTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtctc                                                                                                                                        <  2:2587801/148‑1 (MQ=255)
          gACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTcc                                                                                                                                       >  1:1059621/1‑148 (MQ=255)
             gcgTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCt                                                                                                                                    <  2:3282164/148‑1 (MQ=255)
              cgTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTc                                                                                                                                   >  2:2677658/1‑148 (MQ=255)
                tCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGc                                                                                                                                 >  2:2571043/1‑148 (MQ=255)
                   cGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCAc                                                                                                                              >  1:995125/1‑148 (MQ=255)
                   cGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCAc                                                                                                                              >  2:3416921/1‑148 (MQ=255)
                    gCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTAGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACt                                                                                                                             >  1:1182578/1‑148 (MQ=255)
                     cATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACtt                                                                                                                            >  1:3525527/1‑148 (MQ=255)
                      aTTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTc                                                                                                                           >  1:254862/1‑148 (MQ=255)
                       ttcttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCg                                                                                                                          >  1:1613008/1‑148 (MQ=255)
                        tcttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGa                                                                                                                         <  1:3404042/148‑1 (MQ=255)
                        tcttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGa                                                                                                                         <  1:2275198/148‑1 (MQ=255)
                         cttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGAc                                                                                                                        <  1:89531/148‑1 (MQ=255)
                         cttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGAc                                                                                                                        >  2:2000711/1‑148 (MQ=255)
                           tcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAg                                                                                                                      <  1:3148912/148‑1 (MQ=255)
                           tcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAg                                                                                                                      <  1:812470/148‑1 (MQ=255)
                           tcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAg                                                                                                                      <  1:1383131/148‑1 (MQ=255)
                              gTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtg                                                                                                                   <  1:721000/148‑1 (MQ=255)
                                  aaaCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtggtgg                                                                                                               >  2:2253216/1‑148 (MQ=255)
                                      cTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaa                                                                                                           >  2:3120094/1‑148 (MQ=255)
                                       tGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaaa                                                                                                          <  1:2037837/148‑1 (MQ=255)
                                        gAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaaaa                                                                                                         <  1:151556/148‑1 (MQ=255)
                                         aaaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAg                                                                                                        <  1:2663400/148‑1 (MQ=255)
                                          aaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGt                                                                                                       <  2:459425/148‑1 (MQ=255)
                                          aaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGt                                                                                                       <  1:1776167/148‑1 (MQ=255)
                                          aaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGt                                                                                                       <  1:1853762/148‑1 (MQ=255)
                                           aCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTc                                                                                                      >  1:2305783/1‑148 (MQ=255)
                                             tACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAg                                                                                                    <  1:203817/148‑1 (MQ=255)
                                                    ggtaGTAGAAGATCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAACCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaa                                                                                             <  1:21547/145‑1 (MQ=255)
                                                     gAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaaa                                                                                            <  2:426513/148‑1 (MQ=255)
                                                        gTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGc                                                                                         <  1:599461/148‑1 (MQ=255)
                                                          agaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCt                                                                                       >  1:282117/1‑148 (MQ=255)
                                                           gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                      <  2:3509318/148‑1 (MQ=255)
                                                            aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                     >  2:2631204/1‑148 (MQ=255)
                                                            aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                     <  1:1262611/148‑1 (MQ=255)
                                                            aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                     <  2:757594/148‑1 (MQ=255)
                                                            aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                     <  2:1341390/148‑1 (MQ=255)
                                                             agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                                    <  2:3138832/148‑1 (MQ=255)
                                                             agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                                    <  2:2520376/148‑1 (MQ=255)
                                                              gagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAtt                                                                                   <  1:1943770/148‑1 (MQ=255)
                                                              gagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAtt                                                                                   <  1:144644/148‑1 (MQ=255)
                                                                gCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTg                                                                                 <  2:2078812/148‑1 (MQ=255)
                                                                  ttgAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                               <  1:2340371/146‑1 (MQ=255)
                                                                  tGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                               >  2:3406597/1‑148 (MQ=255)
                                                                    gAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCa                                                                             >  2:319798/1‑148 (MQ=255)
                                                                       cacaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACt                                                                          >  2:2800655/1‑148 (MQ=255)
                                                                       cacaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACt                                                                          <  2:1961407/148‑1 (MQ=255)
                                                                       cacaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACt                                                                          >  2:1379740/1‑148 (MQ=255)
                                                                       cacaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGAGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACt                                                                          >  1:1568858/1‑148 (MQ=255)
                                                                            aTCACGTCGCCGCGATATTCTGGTTGATGACGAAGCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGAAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaa                                                                     >  1:3315857/1‑148 (MQ=255)
                                                                             tCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaa                                                                    >  2:3309300/1‑148 (MQ=255)
                                                                                cGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaaaaa                                                                 <  2:799794/148‑1 (MQ=255)
                                                                                  tCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGc                                                               >  2:177090/1‑148 (MQ=255)
                                                                                   cgccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCa                                                              <  2:1745214/148‑1 (MQ=255)
                                                                                    gccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCAt                                                             <  1:2519845/148‑1 (MQ=255)
                                                                                     ccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATg                                                            <  2:68416/148‑1 (MQ=255)
                                                                                      cgcgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGt                                                           <  2:1098251/148‑1 (MQ=255)
                                                                                      cgcgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGt                                                           <  2:1059621/148‑1 (MQ=255)
                                                                                      cgcgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGt                                                           <  2:3098051/148‑1 (MQ=255)
                                                                                        cgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTg                                                         <  2:1614890/148‑1 (MQ=255)
                                                                                          atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                       <  2:578189/148‑1 (MQ=255)
                                                                                          atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                       <  2:2998225/148‑1 (MQ=255)
                                                                                          atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                       >  1:2971851/1‑148 (MQ=255)
                                                                                          atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                       <  2:2418952/148‑1 (MQ=255)
                                                                                           tatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATc                                                      >  2:1581034/1‑148 (MQ=255)
                                                                                              tCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCaaa                                                   <  1:2321150/148‑1 (MQ=255)
                                                                                              tCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCaaa                                                   >  1:2850370/1‑148 (MQ=255)
                                                                                               cTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAag                                                  <  2:2695622/148‑1 (MQ=255)
                                                                                                  gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                               <  2:2283338/148‑1 (MQ=255)
                                                                                                  gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                               <  2:3408606/148‑1 (MQ=255)
                                                                                                   ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                              >  2:1007365/1‑148 (MQ=255)
                                                                                                   ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                              >  2:2549125/1‑148 (MQ=255)
                                                                                                      atgaCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGgcgc                                           <  2:474291/148‑1 (MQ=255)
                                                                                                      atgaCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGgcgc                                           >  1:1051047/1‑148 (MQ=255)
                                                                                                         aCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                        >  2:2772779/1‑148 (MQ=255)
                                                                                                         aCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                        >  1:2412264/1‑148 (MQ=255)
                                                                                                            aaaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaa                                     <  2:3080139/148‑1 (MQ=255)
                                                                                                            aaaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaa                                     <  2:1524447/148‑1 (MQ=255)
                                                                                                             aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                    >  1:1494077/1‑148 (MQ=255)
                                                                                                             aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                    >  1:3340720/1‑148 (MQ=255)
                                                                                                              aCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaaa                                   >  2:1070699/1‑148 (MQ=255)
                                                                                                                  tgttTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAg                               >  1:1903683/1‑148 (MQ=255)
                                                                                                                           tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                      >  1:1282283/1‑148 (MQ=255)
                                                                                                                           tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                      <  1:465813/148‑1 (MQ=255)
                                                                                                                           tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                      <  1:2987407/148‑1 (MQ=255)
                                                                                                                              aCGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTa                   >  1:684252/1‑148 (MQ=255)
                                                                                                                              aCGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTa                   >  1:2087691/1‑148 (MQ=255)
                                                                                                                                  ccAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCg               <  2:3292426/148‑1 (MQ=255)
                                                                                                                                     gcgcATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAAc            <  1:278143/148‑1 (MQ=255)
                                                                                                                                     gcgcATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAAc            <  1:1872798/148‑1 (MQ=255)
                                                                                                                                      cgcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACt           >  2:189517/1‑148 (MQ=255)
                                                                                                                                       gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt          <  1:236317/148‑1 (MQ=255)
                                                                                                                                       gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt          >  1:381633/1‑148 (MQ=255)
                                                                                                                                       gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATAAGCAAGCTGCATTACCCGAAcaa          >  1:1203466/1‑146 (MQ=255)
                                                                                                                                              ccACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCAt   >  2:2101753/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTGACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCATC  >  CP009273/1479587‑1479877

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: