Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 159,666 G→T 100% 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,6660GT99.4% 628.9 / ‑4.2 196G359C (GGT→TGT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (89/105);  total (90/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 196 total aligned reads displayed.

ACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGC  >  CP009273/159521‑159810
                                                                                                                                                 |                                                                                                                                                
aCGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGt                                                                                                                                                <  1:2181988/148‑1 (MQ=255)
 cgcgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAACGCGATCTTTCCTGtt                                                                                                                                               >  2:2686377/1‑148 (MQ=255)
 cgcgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGtt                                                                                                                                               <  1:923631/148‑1 (MQ=255)
 cgcgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGtt                                                                                                                                               >  2:1487012/1‑148 (MQ=255)
 cgcgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGtt                                                                                                                                               <  2:1848796/148‑1 (MQ=255)
 cgcgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGtt                                                                                                                                               >  1:1818190/1‑148 (MQ=255)
 cgcgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGtt                                                                                                                                               <  1:2835424/148‑1 (MQ=255)
   cgTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCGTACAAAGCGATCTTTCCTGTTTg                                                                                                                                             <  2:2408992/148‑1 (MQ=255)
    gTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGc                                                                                                                                            <  2:757079/148‑1 (MQ=255)
     ttAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCt                                                                                                                                           >  2:3522243/1‑148 (MQ=255)
       aGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGa                                                                                                                                         <  2:3595025/148‑1 (MQ=255)
        gCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGAt                                                                                                                                        <  1:3137876/148‑1 (MQ=255)
        gCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGAt                                                                                                                                        <  2:1809964/148‑1 (MQ=255)
        gCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATATTACAAAGCGATCTTTCCTGTTTGCTGAt                                                                                                                                        >  1:1887571/1‑148 (MQ=255)
            ggCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGaa                                                                                                                                    <  1:279712/148‑1 (MQ=255)
            ggCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGaa                                                                                                                                    <  1:925940/148‑1 (MQ=255)
             gCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAt                                                                                                                                   <  2:1292977/148‑1 (MQ=255)
              cATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAtt                                                                                                                                  <  1:3445086/148‑1 (MQ=255)
               aTCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTa                                                                                                                                 <  2:2785823/148‑1 (MQ=255)
                tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTAc                                                                                                                                <  1:481141/148‑1 (MQ=255)
                tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTAc                                                                                                                                >  2:2795077/1‑148 (MQ=255)
                tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTAc                                                                                                                                <  1:288938/148‑1 (MQ=255)
                tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTAc                                                                                                                                <  1:1769844/148‑1 (MQ=255)
                 ccATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACt                                                                                                                               >  2:237177/1‑148 (MQ=255)
                  cATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACAGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTg                                                                                                                              <  2:3123101/148‑1 (MQ=255)
                   aTGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGc                                                                                                                             >  1:1159962/1‑148 (MQ=255)
                   aTGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGc                                                                                                                             >  1:3080272/1‑148 (MQ=255)
                    tGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTACTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCa                                                                                                                            <  1:3514700/148‑1 (MQ=255)
                    tGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCa                                                                                                                            <  2:1262267/148‑1 (MQ=255)
                      aCGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAAt                                                                                                                          >  2:3312960/1‑148 (MQ=255)
                      aCGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAAt                                                                                                                          >  2:1407081/1‑148 (MQ=255)
                       cGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATg                                                                                                                         <  1:882540/148‑1 (MQ=255)
                          aGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATgggg                                                                                                                      <  1:2735046/148‑1 (MQ=255)
                               gCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCa                                                                                                                 <  1:1352255/148‑1 (MQ=255)
                                  gggCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCg                                                                                                              <  2:2922353/148‑1 (MQ=255)
                                            ggCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCa                                                                                                    >  2:2381482/1‑148 (MQ=255)
                                             gCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAg                                                                                                   >  2:316080/1‑148 (MQ=255)
                                             gCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCAGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAg                                                                                                   <  2:1825413/148‑1 (MQ=255)
                                               gTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGAt                                                                                                 <  1:1971208/148‑1 (MQ=255)
                                                tCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCCATGGGGATGCAGCGATCCGGCGCAGATg                                                                                                >  1:1665784/1‑148 (MQ=255)
                                                tCTGGAGCCGGGTATCAGCCTGCATGTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTAAAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATg                                                                                                >  1:1665833/1‑148 (MQ=255)
                                                 cTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGa                                                                                               >  1:3369104/1‑148 (MQ=255)
                                                  tGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAg                                                                                              <  2:390940/148‑1 (MQ=255)
                                                  tGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAg                                                                                              <  1:184257/148‑1 (MQ=255)
                                                   ggAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGACTTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAgc                                                                                             >  2:731587/1‑148 (MQ=255)
                                                      gCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAgctgg                                                                                          <  1:1870767/148‑1 (MQ=255)
                                                          ggTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGctggctgg                                                                                      >  2:2996909/1‑148 (MQ=255)
                                                           gTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGa                                                                                     <  2:2888978/148‑1 (MQ=255)
                                                           gTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGa                                                                                     <  1:167238/148‑1 (MQ=255)
                                                              tCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCa                                                                                  >  2:2280071/1‑148 (MQ=255)
                                                              tCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCa                                                                                  >  2:957174/1‑148 (MQ=255)
                                                               cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCACTGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCaa                                                                                 >  1:2084367/1‑148 (MQ=255)
                                                               cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATTAGCTGTCTGGATaaa                                                                                 >  1:2939106/1‑148 (MQ=255)
                                                               cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCaa                                                                                 <  1:950725/148‑1 (MQ=255)
                                                               cAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCaa                                                                                 >  2:133897/1‑148 (MQ=255)
                                                                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAAc                                                                                <  1:395907/148‑1 (MQ=255)
                                                                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCag                                                                          <  1:139705/148‑1 (MQ=255)
                                                                         ttAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCagtag                                                                       >  2:1687216/1‑148 (MQ=255)
                                                                            aTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGa                                                                    >  1:1093838/1‑148 (MQ=255)
                                                                            aTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGa                                                                    >  1:1096916/1‑148 (MQ=255)
                                                                            aTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGa                                                                    <  2:363348/148‑1 (MQ=255)
                                                                               gCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATc                                                                 <  2:1786304/148‑1 (MQ=255)
                                                                                 cAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTa                                                               <  1:769676/148‑1 (MQ=255)
                                                                                 cAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTa                                                               <  1:2371988/148‑1 (MQ=255)
                                                                                 cAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTa                                                               <  2:2344785/148‑1 (MQ=255)
                                                                                 cAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTa                                                               <  1:633621/148‑1 (MQ=255)
                                                                                    aGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTGTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTg                                                            >  2:2460229/1‑148 (MQ=255)
                                                                                     gAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTgg                                                           >  2:1886864/1‑148 (MQ=255)
                                                                                      aaCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGc                                                          <  2:539156/148‑1 (MQ=255)
                                                                                      aaCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGc                                                          >  2:2568756/1‑148 (MQ=255)
                                                                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGaa                                                    >  1:247959/1‑148 (MQ=255)
                                                                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGaa                                                    <  2:2900347/148‑1 (MQ=255)
                                                                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGaa                                                    <  2:1709536/148‑1 (MQ=255)
                                                                                               aaCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACg                                                 >  1:2262584/1‑148 (MQ=255)
                                                                                               aaCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACg                                                 >  1:258710/1‑148 (MQ=255)
                                                                                                  gcgcCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCt                                              <  2:1629370/148‑1 (MQ=255)
                                                                                                  gcgcCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCt                                              <  2:1650529/148‑1 (MQ=255)
                                                                                                         ggcgCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaa                                       <  2:1818190/148‑1 (MQ=255)
                                                                                                          gcgcAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaaa                                      <  2:3450652/148‑1 (MQ=255)
                                                                                                                aGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTgg                                >  1:1734653/1‑148 (MQ=255)
                                                                                                                 gTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTggg                               <  1:676801/148‑1 (MQ=255)
                                                                                                                     aCCGGAGATCTTACCAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCa                           <  1:2865017/148‑1 (MQ=255)
                                                                                                                     aCCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCa                           >  2:1434766/1‑148 (MQ=255)
                                                                                                                       cGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACt                         >  1:3544932/1‑148 (MQ=255)
                                                                                                                         gagaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTgg                       >  1:2311112/1‑148 (MQ=255)
                                                                                                                         gagaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTgg                       <  1:420790/148‑1 (MQ=255)
                                                                                                                          agaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGa                      >  2:1457038/1‑148 (MQ=255)
                                                                                                                          agaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGa                      >  2:270472/1‑148 (MQ=255)
                                                                                                                            aTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAgg                    >  1:1565863/1‑148 (MQ=255)
                                                                                                                             tCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAggg                   <  2:3193592/148‑1 (MQ=255)
                                                                                                                               ttACAACGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTg                 <  2:1887604/148‑1 (MQ=255)
                                                                                                                                tACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGa                <  2:195387/148‑1 (MQ=255)
                                                                                                                                 aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa               <  1:3700304/148‑1 (MQ=255)
                                                                                                                                 aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGaa               >  1:3667840/1‑148 (MQ=255)
                                                                                                                                         aTCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAg       >  1:2947526/1‑148 (MQ=255)
                                                                                                                                         aTCTTTCCTGATGGCTGATGGAAGTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACACATGCTGGGGGCACTGGAGGGTGAACGGCTTAg       <  2:83850/148‑1 (MQ=255)
                                                                                                                                          tCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGt      >  2:2485754/1‑148 (MQ=255)
                                                                                                                                          tCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGt      >  2:3379111/1‑148 (MQ=255)
                                                                                                                                              tCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTgcgc  >  2:3605653/1‑148 (MQ=255)
                                                                                                                                              tCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTgcgc  >  1:1512688/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGC  >  CP009273/159521‑159810

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: