Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 159,666 G→T 94.7% G359C (GGT→TGT)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273159,6660GT94.7% 485.9 / 18.7 170G359C (GGT→TGT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (6/3);  new base T (84/77);  total (90/80)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 171 total aligned reads displayed.

AACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCG  >  CP009273/159520‑159809
                                                                                                                                                  |                                                                                                                                               
caCGCGTGAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATATAATCGCCAACGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACCAAGCGATCTTGCCgg                                                                                                                                                <  2:737356/147‑1 (MQ=255)
 aCGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGt                                                                                                                                               <  1:2156476/148‑1 (MQ=255)
    cgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTg                                                                                                                                            <  2:2000747/148‑1 (MQ=255)
    cgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTg                                                                                                                                            <  2:388547/148‑1 (MQ=255)
       tAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCCAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTg                                                                                                                                         >  1:474461/1‑148 (MQ=255)
        aGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGa                                                                                                                                        <  1:1685261/148‑1 (MQ=255)
              gCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAt                                                                                                                                  >  2:452874/1‑148 (MQ=255)
              gCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAt                                                                                                                                  >  1:1969112/1‑148 (MQ=255)
               cATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAtt                                                                                                                                 <  2:653084/148‑1 (MQ=255)
               cATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAtt                                                                                                                                 >  2:1132739/1‑148 (MQ=255)
                aTCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTa                                                                                                                                <  1:1444234/148‑1 (MQ=255)
                 tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTAc                                                                                                                               <  1:1265238/148‑1 (MQ=255)
                  ccATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACt                                                                                                                              <  2:2363343/148‑1 (MQ=255)
                     tGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCa                                                                                                                           <  2:2801801/148‑1 (MQ=255)
                     tGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCa                                                                                                                           >  1:947516/1‑148 (MQ=255)
                         gCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATgg                                                                                                                       <  1:236249/148‑1 (MQ=255)
                         gCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATgg                                                                                                                       <  1:3029849/148‑1 (MQ=255)
                          cAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATggg                                                                                                                      <  1:1083980/148‑1 (MQ=255)
                             cGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGAt                                                                                                                   <  2:1764506/148‑1 (MQ=255)
                             cGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGAt                                                                                                                   >  2:367175/1‑148 (MQ=255)
                              gTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATg                                                                                                                  <  2:2053621/148‑1 (MQ=255)
                                          ccGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGgc                                                                                                      >  1:2582921/1‑148 (MQ=255)
                                              gCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAg                                                                                                  <  2:2708886/148‑1 (MQ=255)
                                              gCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAg                                                                                                  <  2:2708889/148‑1 (MQ=255)
                                              gCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAg                                                                                                  >  1:2251267/1‑148 (MQ=255)
                                                 tCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATg                                                                                               >  1:2393364/1‑148 (MQ=255)
                                                 tCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATg                                                                                               <  2:1837194/148‑1 (MQ=255)
                                                  cTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGa                                                                                              >  1:2462709/1‑148 (MQ=255)
                                                  cTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGAGCCGCGGCGCAAAGTGAACCGTAGATATTACAAAGCGATCTTTCATGTTTGCTGATGGAATAACTGCAATGTGGATGCAGCGATCCGGCGCAGATGa                                                                                              >  2:2573165/1‑148 (MQ=255)
                                                     gAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAgct                                                                                           <  1:2236561/148‑1 (MQ=255)
                                                       gCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAgctgg                                                                                         <  1:2495365/148‑1 (MQ=255)
                                                          gggTATCAGCCTGCATTTAATCGCCCAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAgctggctg                                                                                      >  2:2665271/1‑148 (MQ=255)
                                                             tATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGAt                                                                                   <  2:1194060/148‑1 (MQ=255)
                                                                cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGGCAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTGTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGTTCaa                                                                                >  1:538179/1‑148 (MQ=255)
                                                                cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCaa                                                                                <  1:1698882/148‑1 (MQ=255)
                                                                cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCaa                                                                                <  2:2933233/148‑1 (MQ=255)
                                                                  gCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAAGGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAAcc                                                                              >  1:714236/1‑148 (MQ=255)
                                                                  gCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTAATGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCGGGATCAAcc                                                                              >  1:714160/1‑148 (MQ=255)
                                                                      gCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCa                                                                          >  1:660009/1‑148 (MQ=255)
                                                                      gCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCa                                                                          >  1:660013/1‑148 (MQ=255)
                                                                         tttAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCagta                                                                       >  2:706084/1‑148 (MQ=255)
                                                                          ttAATCGCCAAACAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCagtag                                                                      <  1:1132588/148‑1 (MQ=255)
                                                                            aaTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTg                                                                    >  2:2572935/1‑148 (MQ=255)
                                                                             aTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGTAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGa                                                                   <  2:310676/148‑1 (MQ=255)
                                                                              tCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGaa                                                                  <  1:1361652/148‑1 (MQ=255)
                                                                              tCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGaa                                                                  <  2:1056716/148‑1 (MQ=255)
                                                                              tCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGaa                                                                  <  1:3008597/148‑1 (MQ=255)
                                                                               cGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAAt                                                                 <  1:2167003/148‑1 (MQ=255)
                                                                                gCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATc                                                                <  1:819027/148‑1 (MQ=255)
                                                                                  cAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTa                                                              <  1:264143/148‑1 (MQ=255)
                                                                                   aaaGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTAc                                                             >  1:342697/1‑148 (MQ=255)
                                                                                   aaaGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTAc                                                             >  1:2388864/1‑148 (MQ=255)
                                                                                   aaaGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTAc                                                             >  1:342094/1‑148 (MQ=255)
                                                                                      gAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTgg                                                          >  2:2329668/1‑148 (MQ=255)
                                                                                      gAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTgg                                                          >  1:3085896/1‑148 (MQ=255)
                                                                                        aCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGcc                                                        >  2:1227937/1‑148 (MQ=255)
                                                                                         cAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCcg                                                       >  2:1398679/1‑148 (MQ=255)
                                                                                           aGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCcgcg                                                     >  2:630764/1‑148 (MQ=255)
                                                                                            gCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGa                                                    >  1:1052596/1‑148 (MQ=255)
                                                                                            gCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGa                                                    >  2:1543527/1‑148 (MQ=255)
                                                                                              aGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGaaa                                                  <  1:1566473/148‑1 (MQ=255)
                                                                                              aGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGaaa                                                  <  1:2448297/148‑1 (MQ=255)
                                                                                                 aCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGt                                               <  1:877819/148‑1 (MQ=255)
                                                                                                     gccgcGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGt                                           <  1:1614474/148‑1 (MQ=255)
                                                                                                     gccgcGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGt                                           <  1:3253207/148‑1 (MQ=255)
                                                                                                          ggcgCAAAGTGAACCGGAGATCTTACACAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaa                                      <  2:1689044/148‑1 (MQ=255)
                                                                                                          ggcgCAAAGTGAACCGGAGATCTTACACAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaa                                      <  2:1689140/148‑1 (MQ=255)
                                                                                                          ggcgCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaa                                      <  2:2125959/148‑1 (MQ=255)
                                                                                                           gcgcAAAGTTAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaaa                                     <  2:741302/148‑1 (MQ=255)
                                                                                                            cgcAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAAt                                    >  1:3229147/1‑148 (MQ=255)
                                                                                                            cgcAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAAt                                    <  1:127860/148‑1 (MQ=255)
                                                                                                             gcAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATg                                   <  1:2735036/148‑1 (MQ=255)
                                                                                                             gcAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATg                                   <  1:1016706/148‑1 (MQ=255)
                                                                                                              cAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGc                                  <  1:382921/148‑1 (MQ=255)
                                                                                                              cAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTAAATCTACTGGCCGCGAAACGTCTGTTACAAATGc                                  >  1:360308/1‑148 (MQ=255)
                                                                                                              cAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTAAATCTACTGGCCGCGAAACGTCTGTTACAAATGc                                  >  1:359129/1‑148 (MQ=255)
                                                                                                                aaGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTg                                <  2:420734/148‑1 (MQ=255)
                                                                                                                 aGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTgg                               >  1:3036525/1‑148 (MQ=255)
                                                                                                                   tGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTgggg                             >  1:245222/1‑148 (MQ=255)
                                                                                                                         ggAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTg                       >  2:802476/1‑148 (MQ=255)
                                                                                                                          gagaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTgg                      >  1:2918180/1‑148 (MQ=255)
                                                                                                                           agaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATTGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCTCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGa                     <  2:2469938/148‑1 (MQ=255)
                                                                                                                           agaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGa                     <  1:1141358/148‑1 (MQ=255)
                                                                                                                              tCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAggg                  >  2:1463415/1‑148 (MQ=255)
                                                                                                                                 tACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGa               >  2:3048338/1‑148 (MQ=255)
                                                                                                                                 tACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGa               <  1:3306118/148‑1 (MQ=255)
                                                                                                                                  acacAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa              <  1:818345/148‑1 (MQ=255)
                                                                                                                                  aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa              <  2:1969112/148‑1 (MQ=255)
                                                                                                                                  aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa              >  1:2021626/1‑148 (MQ=255)
                                                                                                                                  aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa              >  1:2021659/1‑148 (MQ=255)
                                                                                                                                  aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa              <  2:301982/148‑1 (MQ=255)
                                                                                                                                  aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa              <  2:1389887/148‑1 (MQ=255)
                                                                                                                                   cAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAAc             >  1:1597936/1‑148 (MQ=255)
                                                                                                                                     aaGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACgg           >  2:1267074/1‑148 (MQ=255)
                                                                                                                                     aaGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACgg           <  1:1024723/148‑1 (MQ=255)
                                                                                                                                      aGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGc          >  1:2868059/1‑148 (MQ=255)
                                                                                                                                      aGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGc          >  2:69668/1‑148 (MQ=255)
                                                                                                                                       gCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCt         <  1:2099504/148‑1 (MQ=255)
                                                                                                                                         gATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTa       <  1:1497479/148‑1 (MQ=255)
                                                                                                                                          aTCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAg      >  1:449177/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCG  >  CP009273/159520‑159809

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: