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,7330GT97.3% 561.0 / 6.7 182D806Y (GAT→TAT) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (77/100);  total (78/104)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 182 total aligned reads displayed.

ACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCAT  >  CP009273/1479590‑1479876
                                                                                                                                               |                                                                                                                                               
  tGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTaa                                                                                                                                           <  1:2072195/148‑1 (MQ=255)
   ggCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAt                                                                                                                                          <  1:450563/148‑1 (MQ=255)
   ggCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAt                                                                                                                                          <  1:448052/148‑1 (MQ=255)
     cAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAtct                                                                                                                                        <  2:3298751/148‑1 (MQ=255)
        aCGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCg                                                                                                                                     <  1:2793262/148‑1 (MQ=255)
         cgcgTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGc                                                                                                                                    >  1:3195622/1‑148 (MQ=255)
         cgcgTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGc                                                                                                                                    >  1:3195771/1‑148 (MQ=255)
                cGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCAc                                                                                                                             <  2:588298/148‑1 (MQ=255)
                 gCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACt                                                                                                                            <  2:2294527/148‑1 (MQ=255)
                   aTTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTc                                                                                                                          >  2:2688353/1‑148 (MQ=255)
                   aTTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTc                                                                                                                          >  1:579182/1‑148 (MQ=255)
                    ttcttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCg                                                                                                                         >  1:765029/1‑148 (MQ=255)
                     tcttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGa                                                                                                                        >  1:1262890/1‑148 (MQ=255)
                      cttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGAc                                                                                                                       <  2:219398/148‑1 (MQ=255)
                          cGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCt                                                                                                                   <  2:2978582/148‑1 (MQ=255)
                           gTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtg                                                                                                                  <  2:1854301/148‑1 (MQ=255)
                            tGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCtgg                                                                                                                 <  1:249351/148‑1 (MQ=255)
                               aaaCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtggtgg                                                                                                              >  2:1584762/1‑148 (MQ=255)
                               aaaCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtggtgg                                                                                                              >  1:1869332/1‑148 (MQ=255)
                                 aCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaa                                                                                                            >  2:1168819/1‑148 (MQ=255)
                                  ccTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaa                                                                                                           >  2:1856023/1‑148 (MQ=255)
                                  ccTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaa                                                                                                           >  2:1970388/1‑148 (MQ=255)
                                    tGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaaa                                                                                                         >  1:2175550/1‑148 (MQ=255)
                                         cTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCa                                                                                                    >  2:1153007/1‑148 (MQ=255)
                                          tACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAg                                                                                                   <  2:3493214/148‑1 (MQ=255)
                                               gCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGTGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCcgcg                                                                                              >  2:3005202/1‑148 (MQ=255)
                                                      tAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGcc                                                                                       <  1:1149498/148‑1 (MQ=255)
                                                       agaagaGCTGTAACACAACTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCt                                                                                      <  2:2219498/148‑1 (MQ=255)
                                                        gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                     >  2:2940270/1‑148 (MQ=255)
                                                         aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                    >  1:2496688/1‑148 (MQ=255)
                                                         aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                    <  1:2352563/148‑1 (MQ=255)
                                                         aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                    <  2:560768/148‑1 (MQ=255)
                                                          agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                                   <  1:61525/148‑1 (MQ=255)
                                                          agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                                   <  1:1157985/148‑1 (MQ=255)
                                                           gagCAGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAtt                                                                                  <  1:1963976/148‑1 (MQ=255)
                                                               tGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                              <  1:1200272/148‑1 (MQ=255)
                                                               tGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCt                                                                              <  2:765029/148‑1 (MQ=255)
                                                                 gAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCa                                                                            >  1:2615998/1‑148 (MQ=255)
                                                                 gAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCa                                                                            >  2:3227889/1‑148 (MQ=255)
                                                                    cacaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGATTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACt                                                                         <  2:1123728/148‑1 (MQ=255)
                                                                    cacaAATCACGACGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCACCTCGCCACTTCTACAGCTGGTGGAAAAAAGTCAGCCGCGAAACTCCTGATTTGCTCAACt                                                                         >  1:1522551/1‑148 (MQ=255)
                                                                      caAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACttt                                                                       >  2:3277882/1‑148 (MQ=255)
                                                                       aaaTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTg                                                                      <  2:2253561/148‑1 (MQ=255)
                                                                       aaaTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTg                                                                      <  2:574505/148‑1 (MQ=255)
                                                                         aTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaa                                                                    >  2:1714791/1‑148 (MQ=255)
                                                                            aCGTCGCTGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaaaa                                                                 <  2:853723/148‑1 (MQ=255)
                                                                            aCGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaaaa                                                                 <  2:2732976/148‑1 (MQ=255)
                                                                               tCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGc                                                              <  2:3048925/148‑1 (MQ=255)
                                                                                  ccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATg                                                           <  2:2485127/148‑1 (MQ=255)
                                                                                     cgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTg                                                        <  2:1287212/148‑1 (MQ=255)
                                                                                     cgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTg                                                        <  2:795729/148‑1 (MQ=255)
                                                                                      gATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGa                                                       <  1:2997469/148‑1 (MQ=255)
                                                                                       atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                      <  2:724372/148‑1 (MQ=255)
                                                                                         atTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCa                                                    >  2:486817/1‑148 (MQ=255)
                                                                                          ttCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGATAAAAGCAGGTTGATCaa                                                   >  1:2035685/1‑148 (MQ=255)
                                                                                            cTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAag                                                 <  2:38226/148‑1 (MQ=255)
                                                                                             tGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAaga                                                <  2:1297700/148‑1 (MQ=255)
                                                                                              ggTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAagag                                               <  1:1327561/148‑1 (MQ=255)
                                                                                               gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                              >  1:2295952/1‑148 (MQ=255)
                                                                                               gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                              >  1:123322/1‑148 (MQ=255)
                                                                                               gTTGATGACGAAACGTTGTTTGAGTTCGACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAACAAGCATGTTGATCAAAGAgg                                              <  2:3348636/148‑1 (MQ=255)
                                                                                                ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                             <  1:70112/148‑1 (MQ=255)
                                                                                                     gaCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAg                                        <  1:416540/148‑1 (MQ=255)
                                                                                                      aCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                       >  2:3146924/1‑148 (MQ=255)
                                                                                                      aCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                       >  1:643856/1‑148 (MQ=255)
                                                                                                       cGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaa                                      >  1:3057174/1‑148 (MQ=255)
                                                                                                         aaaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaa                                    >  2:1555154/1‑148 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                   >  2:2379492/1‑148 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                   >  2:545446/1‑148 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                   <  2:3151914/148‑1 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                   >  1:1061899/1‑148 (MQ=255)
                                                                                                          aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                   >  2:2227967/1‑148 (MQ=255)
                                                                                                           aCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaaa                                  >  2:1105984/1‑148 (MQ=255)
                                                                                                             gttgttTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATc                                >  2:19662/1‑148 (MQ=255)
                                                                                                              ttgttTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCa                               <  1:2940270/148‑1 (MQ=255)
                                                                                                                 tttGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCa                            <  1:2525502/148‑1 (MQ=255)
                                                                                                                    gAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGc                         <  1:2140127/148‑1 (MQ=255)
                                                                                                                    gAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGc                         <  1:2324743/148‑1 (MQ=255)
                                                                                                                     aGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCACGCt                        >  2:3317349/1‑148 (MQ=255)
                                                                                                                      gTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTg                       <  1:3479303/148‑1 (MQ=255)
                                                                                                                       ttCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTgg                      <  1:2601681/148‑1 (MQ=255)
                                                                                                                         cTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGAt                    <  1:1105959/148‑1 (MQ=255)
                                                                                                                           aCGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTa                  <  1:139805/148‑1 (MQ=255)
                                                                                                                           aCGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTa                  <  2:2496688/148‑1 (MQ=255)
                                                                                                                            cGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGACACGCCTGATTTGCTCAACTTGGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTAc                 <  2:1888366/148‑1 (MQ=255)
                                                                                                                             gACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTAcc                >  2:720742/1‑148 (MQ=255)
                                                                                                                                cAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGa             >  1:1741704/1‑148 (MQ=255)
                                                                                                                                 aGCGCATCAGCCACTATGTAATCTCCGCTCTCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGaa            <  1:1885895/148‑1 (MQ=255)
                                                                                                                                 aGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGaa            <  2:1403754/148‑1 (MQ=255)
                                                                                                                                  gcgcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAAc           >  2:802225/1‑148 (MQ=255)
                                                                                                                                    gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt         <  1:2975875/148‑1 (MQ=255)
                                                                                                                                    gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt         <  1:835618/148‑1 (MQ=255)
                                                                                                                                      aTCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCt       <  1:252226/148‑1 (MQ=255)
                                                                                                                                      aTCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCt       <  2:2264535/148‑1 (MQ=255)
                                                                                                                                      aTCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCt       >  2:1042479/1‑148 (MQ=255)
                                                                                                                                       tCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTg      >  1:2365162/1‑148 (MQ=255)
                                                                                                                                       tCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTg      >  1:414472/1‑148 (MQ=255)
                                                                                                                                        cAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTgg     >  1:523698/1‑148 (MQ=255)
                                                                                                                                        cAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATAACCCGAACTTCTgg     >  2:693574/1‑148 (MQ=255)
                                                                                                                                           ccACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCAt  >  2:263495/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
ACTGGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCAT  >  CP009273/1479590‑1479876

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: