Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,570 G→A 5.3% 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,5700GA 5.3% 1335.5 / 42.4 364D85N (GAC→AAC) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (172/171);  new base A (9/10);  total (182/182)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 374 total aligned reads displayed.

AGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCG  >  CP009273/1477421‑1477719
                                                                                                                                                     |                                                                                                                                                     
agaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAgg                                                                                                                                                       >  2:4010578/1‑150 (MQ=255)
  aagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGAc                                                                                                                                                     <  1:1851733/150‑1 (MQ=255)
      ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                                                                                                                                 <  1:2049959/150‑1 (MQ=255)
       tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtct                                                                                                                                                >  2:817637/1‑150 (MQ=255)
         aaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCg                                                                                                                                              <  2:3471450/150‑1 (MQ=255)
         aaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCg                                                                                                                                              <  1:139407/150‑1 (MQ=255)
          aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGa                                                                                                                                             <  1:4976219/150‑1 (MQ=255)
          aaaTCCTGATGCACAACAGGCCATTTACCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCATGCTGCGTGAAGCGTCACGACCGGAAATTACTTATACTGACAATGAACCGGTTAGTCAGAAAAAACATGACATtcacta                                                                                                                                             >  1:1118200/1‑146 (MQ=255)
           aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGaa                                                                                                                                            >  2:3616569/1‑150 (MQ=255)
            aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                                                                                                                                           <  2:4001345/150‑1 (MQ=255)
            aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                                                                                                                                           >  1:2175342/1‑150 (MQ=255)
              ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCg                                                                                                                                         <  1:225236/150‑1 (MQ=255)
                       caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAt                                                                                                                                <  1:1042243/150‑1 (MQ=255)
                         acaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCa                                                                                                                              >  1:2680940/1‑150 (MQ=255)
                           aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc                                                                                                                            >  1:4723000/1‑150 (MQ=255)
                              ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAgg                                                                                                                         >  1:4482763/1‑150 (MQ=255)
                                    tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                   <  1:2323244/150‑1 (MQ=255)
                                    tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                   >  1:2463076/1‑150 (MQ=255)
                                    tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                   >  1:3544380/1‑150 (MQ=255)
                                    tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                                   >  1:927598/1‑150 (MQ=255)
                                      cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                                 <  1:4336766/150‑1 (MQ=255)
                                       aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcg                                                                                                                >  2:4110459/1‑150 (MQ=255)
                                                cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgggg                                                                                                       <  2:2701320/150‑1 (MQ=255)
                                                cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgggg                                                                                                       >  2:3936590/1‑150 (MQ=255)
                                                 gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                                      <  2:299933/150‑1 (MQ=255)
                                                 gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                                      >  1:3443581/1‑150 (MQ=255)
                                                 gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                                      >  1:768449/1‑150 (MQ=255)
                                                   aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaaa                                                                                                    >  1:4626611/1‑150 (MQ=255)
                                                   aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaaa                                                                                                    <  1:2631762/150‑1 (MQ=255)
                                                    agagATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAc                                                                                                   <  2:354156/150‑1 (MQ=255)
                                                      agaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTCTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACgg                                                                                                 <  2:876555/150‑1 (MQ=255)
                                                       gaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACggg                                                                                                >  2:1376944/1‑150 (MQ=255)
                                                       gaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACggg                                                                                                <  2:3527499/150‑1 (MQ=255)
                                                         ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                              <  1:1866928/150‑1 (MQ=255)
                                                          tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                             <  1:3282461/150‑1 (MQ=255)
                                                             ccAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                          >  1:4107717/1‑150 (MQ=255)
                                                              cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                         >  2:1551528/1‑150 (MQ=255)
                                                              cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                         >  2:1950520/1‑150 (MQ=255)
                                                              cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                         <  1:1008799/150‑1 (MQ=255)
                                                              cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                         <  1:3662589/150‑1 (MQ=255)
                                                               aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTa                                                                                        >  1:2762129/1‑150 (MQ=255)
                                                                ggcggcAGGGAAAGCCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaa                                                                                       <  2:3315362/150‑1 (MQ=255)
                                                                  cggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaaaa                                                                                     >  1:2934059/1‑150 (MQ=255)
                                                                    gcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACg                                                                                   >  1:3344259/1‑150 (MQ=255)
                                                                    gcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACg                                                                                   <  1:1462184/150‑1 (MQ=255)
                                                                      aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGAc                                                                                 <  2:503321/150‑1 (MQ=255)
                                                                      aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGAc                                                                                 >  2:3066847/1‑150 (MQ=255)
                                                                        ggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTc                                                                               >  2:4712744/1‑150 (MQ=255)
                                                                        ggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTc                                                                               >  2:2158847/1‑150 (MQ=255)
                                                                            aGTCCTGCTGCGTGAAGCGGCACGACCGTAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGtt                                                                           <  2:4433788/150‑1 (MQ=255)
                                                                            aGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGtt                                                                           <  2:2571893/150‑1 (MQ=255)
                                                                                    tgcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                                                   >  1:1006025/1‑150 (MQ=255)
                                                                                        tGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTg                                                               <  1:3021969/150‑1 (MQ=255)
                                                                                          aaGCGGCACGACCGTAAATTACTTATCCTGACAATTTCCCGTTTTGTCATAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                                             <  1:4749164/150‑1 (MQ=255)
                                                                                          aaGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                                             <  2:253877/150‑1 (MQ=255)
                                                                                          aaGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                                             >  1:702860/1‑150 (MQ=255)
                                                                                           aGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                                            >  1:1479070/1‑150 (MQ=255)
                                                                                             cGGCACGACCGGAATTTACTTCTCCTGACTATTTCCCGGTTAGTCAGAAAAAACCGGACTTGCTCGAAGCGATTCGGGATCACCTGGTGGTGCTCGTCGCCGGGTAAACGGGTGCGGGTAACACGACTCAGTTACCGAAAATCTGTATgg                                                          <  2:1761372/150‑1 (MQ=255)
                                                                                             cGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATgg                                                          >  1:1225472/1‑150 (MQ=255)
                                                                                               gCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTAAGAAAAAACAGGACATTCTCGAAGCTATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAg                                                        <  1:586011/150‑1 (MQ=255)
                                                                                                  cGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTg                                                     >  1:2883159/1‑150 (MQ=255)
                                                                                                   gACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                                                    <  2:1651432/150‑1 (MQ=255)
                                                                                                   gACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                                                    >  2:1715464/1‑150 (MQ=255)
                                                                                                   gACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                                                    >  2:1884391/1‑150 (MQ=255)
                                                                                                     ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                                  >  2:449521/1‑150 (MQ=255)
                                                                                                       tgAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                                <  2:2872964/149‑1 (MQ=255)
                                                                                                       ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCGGGGgcg                                                >  2:2221709/1‑150 (MQ=255)
                                                                                                        gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                               >  1:3807660/1‑150 (MQ=255)
                                                                                                        gAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                               >  2:1057278/1‑150 (MQ=255)
                                                                                                          aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCgg                                             >  1:1597369/1‑150 (MQ=255)
                                                                                                            ttACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGa                                           <  1:3223526/150‑1 (MQ=255)
                                                                                                              aCTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAtt                                         <  1:4317089/150‑1 (MQ=255)
                                                                                                                ttATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaa                                       <  2:4626360/150‑1 (MQ=255)
                                                                                                                ttATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaa                                       <  2:37719/150‑1 (MQ=255)
                                                                                                                 tATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                      <  2:1225472/150‑1 (MQ=255)
                                                                                                                 tATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                      >  2:5109597/1‑150 (MQ=255)
                                                                                                                  aTCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                                     >  2:1042631/1‑150 (MQ=255)
                                                                                                                    ccTGACAATTTACCGGTTTGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGa                                   <  1:3558815/150‑1 (MQ=255)
                                                                                                                      tGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACt                                 >  2:4567008/1‑150 (MQ=255)
                                                                                                                         cAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGACGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                              <  2:1512112/150‑1 (MQ=255)
                                                                                                                           aTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                            <  1:4135450/150‑1 (MQ=255)
                                                                                                                             ttACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGc                          >  2:4960371/1‑150 (MQ=255)
                                                                                                                              tACCGGTTAGTCAGAAAAAACATGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGcc                         >  2:2493532/1‑150 (MQ=255)
                                                                                                                               aCCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCa                        <  1:594494/150‑1 (MQ=255)
                                                                                                                                ccGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt                       >  1:4231422/1‑150 (MQ=255)
                                                                                                                                ccGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt                       >  2:2841464/1‑150 (MQ=255)
                                                                                                                                 cGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATa                      >  2:586609/1‑150 (MQ=255)
                                                                                                                                  ggTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                     <  2:3443581/150‑1 (MQ=255)
                                                                                                                                  ggTTAGTCAGAAAAAACAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc                     <  1:1754781/150‑1 (MQ=255)
                                                                                                                                   gTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAcc                    >  1:983934/1‑150 (MQ=255)
                                                                                                                                    ttAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAccc                   >  2:2366081/1‑150 (MQ=255)
                                                                                                                                     tAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCa                  <  1:4940493/150‑1 (MQ=255)
                                                                                                                                      aGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAg                 <  1:4154402/150‑1 (MQ=255)
                                                                                                                                         cAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg              <  1:817637/150‑1 (MQ=255)
                                                                                                                                          aGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgc             >  1:3867522/1‑150 (MQ=255)
                                                                                                                                                aaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCgtcgtc       >  1:3978147/1‑150 (MQ=255)
                                                                                                                                                 aCAGAACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCt      >  2:1904609/1‑150 (MQ=255)
                                                                                                                                                  cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTg     >  2:1987406/1‑150 (MQ=255)
                                                                                                                                                  cAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTg     >  1:1035068/1‑150 (MQ=255)
                                                                                                                                                    ggACATTCTCGAAGCGATTCGTGATCCCCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTggc   <  2:4203072/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGTCTGGCG  >  CP009273/1477421‑1477719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: