Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,479,733 G→T 91.7% 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,7330GT91.7% 390.7 / 28.1 145D806Y (GAT→TAT) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base T (52/80);  total (59/86)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

ACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTG  >  CP009273/1479590‑1479872
                                                                                                                                               |                                                                                                                                           
aCTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGt                                                                                                                                         <  1:2217062/148‑1 (MQ=255)
    gCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtc                                                                                                                                     <  2:2797282/148‑1 (MQ=255)
     cAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtct                                                                                                                                    <  2:1286805/148‑1 (MQ=255)
     cAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtct                                                                                                                                    <  2:1287184/148‑1 (MQ=255)
     cAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtct                                                                                                                                    <  1:2154829/148‑1 (MQ=255)
     cAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtct                                                                                                                                    <  2:3074884/148‑1 (MQ=255)
        aCGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCg                                                                                                                                 <  2:2850525/148‑1 (MQ=255)
             tCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGc                                                                                                                            <  2:841593/148‑1 (MQ=255)
              cACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGcc                                                                                                                           >  2:2585738/1‑148 (MQ=255)
              cACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGcc                                                                                                                           >  2:1588456/1‑148 (MQ=255)
               aCGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCa                                                                                                                          <  2:1829687/148‑1 (MQ=255)
                cGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCAc                                                                                                                         <  1:1696782/148‑1 (MQ=255)
                cGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCAc                                                                                                                         <  1:1699213/148‑1 (MQ=255)
                  cATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACtt                                                                                                                       <  2:801846/148‑1 (MQ=255)
                  cATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACtt                                                                                                                       <  1:1772571/148‑1 (MQ=255)
                  cATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACtt                                                                                                                       <  1:1772881/148‑1 (MQ=255)
                   aTTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTc                                                                                                                      >  2:2428249/1‑148 (MQ=255)
                    ttcttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCg                                                                                                                     <  1:2810865/148‑1 (MQ=255)
                      cttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGAc                                                                                                                   <  2:1140601/148‑1 (MQ=255)
                       ttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACa                                                                                                                  <  1:1561653/148‑1 (MQ=255)
                       ttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCCCTCGCCACTTCGACa                                                                                                                  <  1:1561635/148‑1 (MQ=255)
                            tGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtgg                                                                                                             <  2:2881519/148‑1 (MQ=255)
                            tGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtgg                                                                                                             >  2:1933800/1‑148 (MQ=255)
                             gAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtggt                                                                                                            <  1:1686448/148‑1 (MQ=255)
                                 aCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaa                                                                                                        <  2:362731/148‑1 (MQ=255)
                                  ccTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaa                                                                                                       <  2:3216644/148‑1 (MQ=255)
                                   cTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaa                                                                                                      <  2:2699329/148‑1 (MQ=255)
                                    tGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGCCAGCTGGTGGaaaaa                                                                                                     >  2:370143/1‑148 (MQ=255)
                                     gAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaaaa                                                                                                    <  1:1265290/148‑1 (MQ=255)
                                      aaaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAg                                                                                                   <  1:1429997/148‑1 (MQ=255)
                                       aaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGt                                                                                                  <  1:496499/148‑1 (MQ=255)
                                        aCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTc                                                                                                 >  2:408328/1‑148 (MQ=255)
                                           aCGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGc                                                                                              >  2:1250973/1‑148 (MQ=255)
                                                 ggAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaa                                                                                        <  1:1065804/148‑1 (MQ=255)
                                                  gAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaaa                                                                                       <  2:2159709/148‑1 (MQ=255)
                                                  gAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaaa                                                                                       <  2:496451/148‑1 (MQ=255)
                                                  gAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaaa                                                                                       <  1:3167649/148‑1 (MQ=255)
                                                  gAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaaa                                                                                       <  1:3167714/148‑1 (MQ=255)
                                                   aaGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAAc                                                                                      >  2:303661/1‑148 (MQ=255)
                                                      tAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGcc                                                                                   <  1:381204/148‑1 (MQ=255)
                                                        gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                 <  1:2428249/148‑1 (MQ=255)
                                                        gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                 >  1:2903731/1‑148 (MQ=255)
                                                        gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                 >  1:2903714/1‑148 (MQ=255)
                                                          agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                               <  1:2208996/148‑1 (MQ=255)
                                                          agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                               <  1:812748/148‑1 (MQ=255)
                                                           gagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAtt                                                                              >  2:1220233/1‑148 (MQ=255)
                                                           gagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAtt                                                                              >  1:1371513/1‑148 (MQ=255)
                                                             gCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTg                                                                            <  2:849862/148‑1 (MQ=255)
                                                              cTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGc                                                                           >  1:362768/1‑148 (MQ=255)
                                                              cTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGc                                                                           >  1:3270668/1‑148 (MQ=255)
                                                               tGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                          <  1:585236/148‑1 (MQ=255)
                                                               tGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                          >  2:816682/1‑148 (MQ=255)
                                                               tGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                          >  2:816590/1‑148 (MQ=255)
                                                                  aaCACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCaa                                                                       <  2:2377277/148‑1 (MQ=255)
                                                                  aaCACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCaa                                                                       <  2:2377114/148‑1 (MQ=255)
                                                                  aaCACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCaa                                                                       <  1:1667291/148‑1 (MQ=255)
                                                                  aaCACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCaa                                                                       <  1:3066877/148‑1 (MQ=255)
                                                                  aaCACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCaa                                                                       <  1:3066823/148‑1 (MQ=255)
                                                                     acaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACtt                                                                    <  1:2251587/148‑1 (MQ=255)
                                                                      caAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACttt                                                                   <  2:2837227/148‑1 (MQ=255)
                                                                       aaaTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTg                                                                  >  1:1056407/1‑148 (MQ=255)
                                                                       aaaTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTg                                                                  >  2:2556559/1‑148 (MQ=255)
                                                                         aTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaa                                                                <  1:1426723/148‑1 (MQ=255)
                                                                               tCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGc                                                          >  1:2235691/1‑148 (MQ=255)
                                                                                cgccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCa                                                         >  1:746570/1‑148 (MQ=255)
                                                                                cgccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCa                                                         <  1:1727759/148‑1 (MQ=255)
                                                                                cgccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCa                                                         <  1:819775/148‑1 (MQ=255)
                                                                                 gccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCAt                                                        <  2:2934003/148‑1 (MQ=255)
                                                                                 gccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCAt                                                        <  2:1445528/148‑1 (MQ=255)
                                                                                  ccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCagg                                                       >  1:2175990/1‑146 (MQ=255)
                                                                                  ccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATg                                                       >  1:2175942/1‑148 (MQ=255)
                                                                                   cgcgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGt                                                      <  1:3168038/148‑1 (MQ=255)
                                                                                     cgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTg                                                    <  2:669984/148‑1 (MQ=255)
                                                                                     cgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACGTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTg                                                    <  2:1870923/148‑1 (MQ=255)
                                                                                      gATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGa                                                   <  1:2784267/148‑1 (MQ=255)
                                                                                       atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                  <  1:1360197/148‑1 (MQ=255)
                                                                                         atTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCa                                                >  1:103947/1‑148 (MQ=255)
                                                                                            cTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAag                                             <  1:1804337/148‑1 (MQ=255)
                                                                                               gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                          <  2:3222924/148‑1 (MQ=255)
                                                                                                ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                         >  1:2883173/1‑148 (MQ=255)
                                                                                                ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                         <  1:2322201/148‑1 (MQ=255)
                                                                                                 tgatgaCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGgc                                        <  2:348269/148‑1 (MQ=255)
                                                                                                       cGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaa                                  >  2:117899/1‑148 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                               <  1:1220233/148‑1 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                               <  1:3013652/148‑1 (MQ=255)
                                                                                                           aCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaaa                              <  2:2697413/148‑1 (MQ=255)
                                                                                                            cGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACAATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATAAAAGAGGGCGCAGAAAAAAt                             >  2:221845/1‑148 (MQ=255)
                                                                                                                   tGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAg                      >  1:700890/1‑148 (MQ=255)
                                                                                                                    gAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGc                     <  2:1939273/148‑1 (MQ=255)
                                                                                                                     aGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCt                    >  2:1379863/1‑148 (MQ=255)
                                                                                                                     aGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCt                    >  1:504466/1‑148 (MQ=255)
                                                                                                                       ttCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTgg                  >  2:400999/1‑148 (MQ=255)
                                                                                                                        tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                 >  2:166691/1‑148 (MQ=255)
                                                                                                                        tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                 <  2:3124586/148‑1 (MQ=255)
                                                                                                                           aCGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTa              <  1:3169619/148‑1 (MQ=255)
                                                                                                                               ccAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCg          >  2:1534022/1‑148 (MQ=255)
                                                                                                                                cAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGa         >  2:315685/1‑148 (MQ=255)
                                                                                                                                   cgcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACt      >  1:1948477/1‑148 (MQ=255)
                                                                                                                                    gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt     <  2:504466/148‑1 (MQ=255)
                                                                                                                                       tCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTg  <  1:2983147/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
ACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTG  >  CP009273/1479590‑1479872

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: