Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,570 G→A 46.4% D85N (GAC→AAC)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,477,5700GA46.4% 67.7 / inf 263D85N (GAC→AAC) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (75/66);  new base A (61/61);  total (136/127)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 264 total aligned reads displayed.

AGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGG  >  CP009273/1477421‑1477717
                                                                                                                                                     |                                                                                                                                                   
agaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAgg                                                                                                                                                     <  1:2966934/150‑1 (MQ=255)
   agGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACa                                                                                                                                                  >  1:4046742/1‑150 (MQ=255)
    ggTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAt                                                                                                                                                 <  1:1145454/150‑1 (MQ=255)
     gTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGGACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAtt                                                                                                                                                <  1:1384484/150‑1 (MQ=255)
        aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtctc                                                                                                                                             >  2:2072510/1‑150 (MQ=255)
        aaaaaTCCTGATGCACAACAGGCCATTGTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtctc                                                                                                                                             <  2:4163655/150‑1 (MQ=255)
          aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGa                                                                                                                                           >  2:130615/1‑150 (MQ=255)
               cTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGa                                                                                                                                      >  2:1204592/1‑150 (MQ=255)
               cTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGa                                                                                                                                      >  2:1204677/1‑150 (MQ=255)
                tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                                                                                                                                     <  1:828112/150‑1 (MQ=255)
                tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                                                                                                                                     <  1:827241/150‑1 (MQ=255)
                tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGAt                                                                                                                                     <  2:3525710/150‑1 (MQ=255)
                   tGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCg                                                                                                                                  <  1:438578/150‑1 (MQ=255)
                        aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTAACGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATc                                                                                                                             >  2:2491889/1‑150 (MQ=255)
                           aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc                                                                                                                          >  2:1013185/1‑150 (MQ=255)
                              ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAgg                                                                                                                       >  2:1503869/1‑150 (MQ=255)
                               cATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAggt                                                                                                                      >  1:3314358/1‑150 (MQ=255)
                                aTTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAggtg                                                                                                                     >  1:1193074/1‑150 (MQ=255)
                                 ttttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTAATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtgg                                                                                                                    >  2:2631895/1‑150 (MQ=255)
                                    tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                 >  1:2090693/1‑150 (MQ=255)
                                    tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                 >  1:4331674/1‑150 (MQ=255)
                                      cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                               >  1:2469819/1‑150 (MQ=255)
                                      cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                               <  1:2653627/150‑1 (MQ=255)
                                      cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                               <  1:703837/150‑1 (MQ=255)
                                       aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcg                                                                                                              >  2:1711090/1‑150 (MQ=255)
                                       aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcg                                                                                                              >  1:1222749/1‑150 (MQ=255)
                                       aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcg                                                                                                              >  2:1711061/1‑150 (MQ=255)
                                        ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAATCTGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgt                                                                                                             >  1:719857/1‑150 (MQ=255)
                                        ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgt                                                                                                             >  2:2338191/1‑150 (MQ=255)
                                          agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                           >  2:711029/1‑150 (MQ=255)
                                              ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgg                                                                                                       >  1:498213/1‑150 (MQ=255)
                                                cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgggg                                                                                                     <  1:3145159/150‑1 (MQ=255)
                                                  aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                                   <  2:2936287/150‑1 (MQ=255)
                                                  aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                                   <  1:810934/150‑1 (MQ=255)
                                                  aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                                   >  1:3908489/1‑150 (MQ=255)
                                                  aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                                   >  1:3907688/1‑150 (MQ=255)
                                                   aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCTGAAATTACTTATCCTGACACTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaaa                                                                                                  <  2:3845323/150‑1 (MQ=255)
                                                        aTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGt                                                                                             >  2:1837748/1‑150 (MQ=255)
                                                         ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                            >  1:3501408/1‑150 (MQ=255)
                                                         ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                            >  1:3500235/1‑150 (MQ=255)
                                                         ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                            <  2:3099606/150‑1 (MQ=255)
                                                          tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                           >  2:3576425/1‑150 (MQ=255)
                                                            aCCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTg                                                                                         <  2:783272/150‑1 (MQ=255)
                                                              cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                       <  1:1373573/150‑1 (MQ=255)
                                                               aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTa                                                                                      <  2:859742/150‑1 (MQ=255)
                                                                  cggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaaaa                                                                                   <  2:3628209/150‑1 (MQ=255)
                                                                   ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTCAAAc                                                                                  >  2:4294898/1‑150 (MQ=255)
                                                                    gcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACg                                                                                 <  1:1678494/150‑1 (MQ=255)
                                                                      aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGACACGGGTTCTGGTAAAACGAc                                                                               >  1:4350547/1‑150 (MQ=255)
                                                                       gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCACAAAAAAAAGGACATTCTCGAAGCGATTCTTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACt                                                                              <  1:1799408/150‑1 (MQ=255)
                                                                          aaaGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAg                                                                           >  1:2784454/1‑150 (MQ=255)
                                                                            aGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGtt                                                                         >  1:1061134/1‑150 (MQ=255)
                                                                             gTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTa                                                                        <  1:1522830/150‑1 (MQ=255)
                                                                              tCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAc                                                                       <  2:3041795/150‑1 (MQ=255)
                                                                              tCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAc                                                                       <  2:2579952/150‑1 (MQ=255)
                                                                                        tGAAGCGGCACGACCGTAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCTAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTg                                                             <  1:387149/150‑1 (MQ=255)
                                                                                         gAAGCGGCACGACCGTAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGt                                                            <  2:4070855/150‑1 (MQ=255)
                                                                                         gAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGACAATCTGt                                                            >  2:2837174/1‑150 (MQ=255)
                                                                                           aGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                                          <  2:1210098/150‑1 (MQ=255)
                                                                                              ggCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGa                                                       >  1:3972220/1‑150 (MQ=255)
                                                                                                  cGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTg                                                   >  2:2243494/1‑150 (MQ=255)
                                                                                                  cGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTg                                                   >  1:1986036/1‑150 (MQ=255)
                                                                                                     ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCATAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                                <  1:3148244/150‑1 (MQ=255)
                                                                                                     ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                                >  2:93197/1‑150 (MQ=255)
                                                                                                      cGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgc                                               <  2:1871702/150‑1 (MQ=255)
                                                                                                       ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                              >  2:4218815/1‑150 (MQ=255)
                                                                                                        gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                             >  1:1060028/1‑150 (MQ=255)
                                                                                                         aaaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGcgcg                                            >  1:3419651/1‑150 (MQ=255)
                                                                                                         aaaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGcgcg                                            <  1:4293256/150‑1 (MQ=255)
                                                                                                         aaaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGcgcg                                            >  2:4280627/1‑150 (MQ=255)
                                                                                                            ttACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGa                                         >  1:3800320/1‑150 (MQ=255)
                                                                                                               cTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTa                                      <  2:3652350/150‑1 (MQ=255)
                                                                                                               cTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTa                                      >  2:570696/1‑150 (MQ=255)
                                                                                                                 tATCCTGCCAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                    <  2:3290041/150‑1 (MQ=255)
                                                                                                                  aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGACACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                   >  2:3154998/1‑150 (MQ=255)
                                                                                                                  aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                   >  2:3936885/1‑150 (MQ=255)
                                                                                                                       gACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg                              >  1:2121131/1‑150 (MQ=255)
                                                                                                                       gACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg                              >  1:1709953/1‑150 (MQ=255)
                                                                                                                           aTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                          >  2:1943632/1‑150 (MQ=255)
                                                                                                                             ttACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGc                        >  2:4336016/1‑150 (MQ=255)
                                                                                                                                ccGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt                     >  1:3775445/1‑150 (MQ=255)
                                                                                                                                  ggTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                   >  2:3360099/1‑150 (MQ=255)
                                                                                                                                  ggTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                   <  1:2414437/150‑1 (MQ=255)
                                                                                                                                  ggTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                   >  2:1813669/1‑150 (MQ=255)
                                                                                                                                   gTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAcc                  <  2:3500240/150‑1 (MQ=255)
                                                                                                                                    ttAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAccc                 <  1:1941022/150‑1 (MQ=255)
                                                                                                                                       gTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGc              >  2:3801264/1‑150 (MQ=255)
                                                                                                                                       gTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGc              >  1:3657175/1‑150 (MQ=255)
                                                                                                                                            aaaaaaCAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt         >  1:1886666/1‑150 (MQ=255)
                                                                                                                                            aaaaaaCAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt         >  1:3235373/1‑150 (MQ=255)
                                                                                                                                             aaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc        >  1:3135464/1‑150 (MQ=255)
                                                                                                                                             aaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc        <  2:2662641/150‑1 (MQ=255)
                                                                                                                                              aaaaCAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtcg       >  1:2383777/1‑150 (MQ=255)
                                                                                                                                              aaaaCAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtcg       >  1:2264923/1‑150 (MQ=255)
                                                                                                                                              aaaaCAGAACATACTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtcg       >  2:2865279/1‑150 (MQ=255)
                                                                                                                                                aaCAGACCATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCgtcgtc     <  1:4288842/150‑1 (MQ=255)
                                                                                                                                                aaCAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCgtcgtc     <  1:4288047/150‑1 (MQ=255)
                                                                                                                                                  cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTg   >  2:1780049/1‑150 (MQ=255)
                                                                                                                                                  cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTATGGTAAAACGAATCATTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACAGATCGGCCATACCCAGCCGCgtcgtcgg   >  1:1907802/1‑148 (MQ=255)
                                                                                                                                                   aGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTgg  >  2:4383096/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGG  >  CP009273/1477421‑1477717

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: