Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,565 A→T 68.9% K83I (AAA→ATA)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,477,5650AT68.9% 448.0 / 279.9 324K83I (AAA→ATA) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (48/53);  new base T (106/117);  total (154/170)
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

Only 100 of 334 total aligned reads displayed.

GTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGT  >  CP009273/1477417‑1477713
                                                                                                                                                    |                                                                                                                                                    
gTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaaaa                                                                                                                                                     >  2:4107149/1‑150 (MQ=255)
 tGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTCTCCTGACAATTTACCGGTTAGTCAGAAAAAAc                                                                                                                                                    <  1:1705216/150‑1 (MQ=255)
 tGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATAc                                                                                                                                                    <  2:2859069/150‑1 (MQ=255)
   aagaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAg                                                                                                                                                  >  2:3277589/1‑150 (MQ=255)
    agaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAgg                                                                                                                                                 >  2:2936582/1‑150 (MQ=255)
      aagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGAc                                                                                                                                               <  1:2664560/150‑1 (MQ=255)
       agGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACa                                                                                                                                              <  2:4048213/150‑1 (MQ=255)
        ggTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACTATTTACCGGTTAGTCAGAAAATACAGGACAt                                                                                                                                             >  1:1183225/1‑150 (MQ=255)
         gTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGTGATGGCGAAAGAGATTTACCCGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGATATTTCTTATCCTGACTATTTACCGGTTAGTCAGAAAAAACAGGACAtt                                                                                                                                            <  1:5072665/150‑1 (MQ=255)
          ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtc                                                                                                                                           <  1:1004847/150‑1 (MQ=255)
           tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtct                                                                                                                                          <  2:1130467/150‑1 (MQ=255)
              aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGa                                                                                                                                       <  2:2533297/150‑1 (MQ=255)
                 tCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGc                                                                                                                                    >  2:385811/1‑150 (MQ=255)
                  ccTGATGCACAACAGGCCAGTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAAGGACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCg                                                                                                                                   <  2:130180/150‑1 (MQ=255)
                    tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGAt                                                                                                                                 >  1:783713/1‑150 (MQ=255)
                    tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGAt                                                                                                                                 >  1:3230694/1‑150 (MQ=255)
                     gATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAtt                                                                                                                                <  1:1211886/150‑1 (MQ=255)
                     gATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAtt                                                                                                                                >  1:4321231/1‑150 (MQ=255)
                         cacaACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTg                                                                                                                            >  1:3822078/1‑150 (MQ=255)
                           caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAt                                                                                                                          <  2:4525670/150‑1 (MQ=255)
                            aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATc                                                                                                                         <  2:1762687/150‑1 (MQ=255)
                            aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATc                                                                                                                         <  1:4407815/150‑1 (MQ=255)
                               aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc                                                                                                                      <  2:379348/150‑1 (MQ=255)
                               aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc                                                                                                                      <  2:771688/150‑1 (MQ=255)
                                ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                     >  2:4562870/1‑150 (MQ=255)
                                 gCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAg                                                                                                                    <  2:104420/150‑1 (MQ=255)
                                       ttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCTGgtggtg                                                                                                              <  1:4126714/150‑1 (MQ=255)
                                       ttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGgtggtg                                                                                                              >  2:1868776/1‑150 (MQ=255)
                                       ttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGgtggtg                                                                                                              >  1:1465339/1‑150 (MQ=255)
                                        tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                             <  1:1629538/150‑1 (MQ=255)
                                          cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                           <  2:1532654/150‑1 (MQ=255)
                                           aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcg                                                                                                          >  1:3397456/1‑150 (MQ=255)
                                             gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtc                                                                                                        <  1:1594314/150‑1 (MQ=255)
                                              agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                       >  1:4056687/1‑150 (MQ=255)
                                              agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                       <  1:3644771/150‑1 (MQ=255)
                                               gaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGc                                                                                                      >  1:3115629/1‑150 (MQ=255)
                                                aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGcc                                                                                                     <  2:4817985/150‑1 (MQ=255)
                                                  ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgg                                                                                                   >  1:460477/1‑150 (MQ=255)
                                                    cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgggg                                                                                                 <  2:2345295/150‑1 (MQ=255)
                                                     gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                                >  1:5290953/1‑150 (MQ=255)
                                                      aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                               <  2:1222051/150‑1 (MQ=255)
                                                      aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                               >  2:3549871/1‑150 (MQ=255)
                                                       aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaaa                                                                                              >  2:3744376/1‑150 (MQ=255)
                                                        agagATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAc                                                                                             <  2:1372661/150‑1 (MQ=255)
                                                         gagaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAca                                                                                            <  2:5208750/150‑2 (MQ=255)
                                                         gagaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACg                                                                                            >  2:2116734/1‑150 (MQ=255)
                                                          agaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACgg                                                                                           >  1:739134/1‑150 (MQ=255)
                                                           gaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACggg                                                                                          >  1:3924178/1‑150 (MQ=255)
                                                             ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                        >  1:2083525/1‑150 (MQ=255)
                                                             ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                        <  1:2371202/150‑1 (MQ=255)
                                                              tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                       <  2:1612871/150‑1 (MQ=255)
                                                               gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCt                                                                                      <  2:4955226/150‑1 (MQ=255)
                                                                  cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGt                                                                                   <  1:3029154/150‑1 (MQ=255)
                                                                   aGGCGGCTGGTAAAGTCCTTCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTa                                                                                  <  1:4492948/150‑1 (MQ=255)
                                                                       ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAAc                                                                              >  1:2268793/1‑150 (MQ=255)
                                                                           gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACt                                                                          <  2:2996863/150‑1 (MQ=255)
                                                                             gAAAGTCCTGCTGCGTGAAGCGGCACTACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAGGGGTTCTGGTAAAACGACTCa                                                                        <  2:263239/150‑1 (MQ=255)
                                                                             gAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCa                                                                        >  1:2603767/1‑150 (MQ=255)
                                                                                 gTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTa                                                                    >  2:4344468/1‑150 (MQ=255)
                                                                                     tgctgcGTGAAGCGGCACGACCGGAAATTACTTATCCCGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGACCCGGGTTCTGGTAATACGACTCAGTTACCGa                                                                >  2:4005199/1‑150 (MQ=255)
                                                                                        tgcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                                             >  2:4255346/1‑150 (MQ=255)
                                                                                           gTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCt                                                          >  1:587719/1‑150 (MQ=255)
                                                                                            tGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTg                                                         >  2:5264570/1‑150 (MQ=255)
                                                                                                gCGGCACGACCGTAAATTACTTATCCTGACAATTTACCGTTTAGTAATAAAATAAAGGACATTCTCTAAGCGATTCGTTATCACCTGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGTTAAAACGACTCAGTTACCGAAAATCTGTATg                                                     <  1:2162606/150‑1 (MQ=255)
                                                                                                     aCGACCGGAAATTACTTATCCTGAAAATTTACCGTTTAGACACAAAAAACCGGACATTCTCTAAGCTATTCGTTATCACCAGGTGGTGATCGTCGCCGGGTAAAGGGGTTCTGTTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCt                                                <  1:1914102/150‑1 (MQ=255)
                                                                                                        aCCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTggg                                             >  2:2776822/1‑150 (MQ=255)
                                                                                                         ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                            <  2:402018/150‑1 (MQ=255)
                                                                                                          ggtaaaTTACTTTTCCTGACAATTTACCGTTTTTTAATAAATTAACGTACTTCCCCTAAGCTTTGGGTTTTCACCTGGTGGTGATCGTCGCGGGGGAATCGGTTTCTGGTAACACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgc                                           <  2:2280043/147‑1 (MQ=255)
                                                                                                             aaaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGcgcg                                        >  1:3542906/1‑150 (MQ=255)
                                                                                                               aTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCggg                                      >  1:1532383/1‑150 (MQ=255)
                                                                                                                 tACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCACCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAt                                    >  2:1691514/1‑150 (MQ=255)
                                                                                                                    ttATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaa                                 >  2:1954467/1‑150 (MQ=255)
                                                                                                                     taaccTGACAATTTCCCGGTTAGTCAGAAAAAACAGGACATTCTCGACGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                <  1:4849895/147‑1 (MQ=255)
                                                                                                                     tATCCTGACAATTTACCGGTTAGTCAGAAAATACAGTACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaaa                                >  2:3436928/1‑150 (MQ=255)
                                                                                                                     tATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGCGCGCGGGATTaaa                                >  1:1818011/1‑150 (MQ=255)
                                                                                                                        ccTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGa                             >  2:4991603/1‑150 (MQ=255)
                                                                                                                         cTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCGGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc                            <  2:3618348/150‑1 (MQ=255)
                                                                                                                         cTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc                            >  2:2229767/1‑150 (MQ=255)
                                                                                                                          tGACAATTTCCCGTTTAGTCATAAAATACCGGAATTTCTCGACGCGATTCGTTATCACCTGGTGGTGATCGTCGCCGGGTAAACGGTTTCTGTTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACt                           <  1:648860/150‑1 (MQ=255)
                                                                                                                            aCAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGa                         >  1:5042291/1‑150 (MQ=255)
                                                                                                                             cAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                        <  1:5266956/150‑1 (MQ=255)
                                                                                                                              aaTTTACCGTTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATc                       <  1:1295558/150‑1 (MQ=255)
                                                                                                                              aaTTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATc                       >  2:4712510/1‑150 (MQ=255)
                                                                                                                               aTTTACCGGTTAGTCTGAAACAACAGGCCATTCTCGACGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGACACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                      <  2:1745126/150‑1 (MQ=255)
                                                                                                                               aTTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                      >  2:2915742/1‑150 (MQ=255)
                                                                                                                               aTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                      <  2:582362/150‑1 (MQ=255)
                                                                                                                                    ccGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt                 <  2:4205314/150‑1 (MQ=255)
                                                                                                                                    ccGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt                 >  2:3627737/1‑150 (MQ=255)
                                                                                                                                      tgTTAGCAAGAAAATACAGGCCATTCTCGACGCGATTCTTTATCACCAGGTGGTGATCTTCGCGGGTGAAACTGGTTATGTTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc               <  1:270127/149‑1 (MQ=255)
                                                                                                                                      ggTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc               >  1:5115566/1‑150 (MQ=255)
                                                                                                                                      ggTTAGACAGAAAATACAGGACATTCTCGAAGCTATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc               <  2:3374184/150‑1 (MQ=255)
                                                                                                                                         tAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCa            >  2:673940/1‑150 (MQ=255)
                                                                                                                                         tAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCa            >  1:4300889/1‑150 (MQ=255)
                                                                                                                                         tAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAACACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCa            <  2:5025561/150‑1 (MQ=255)
                                                                                                                                          aGTCAGAAAAAACTGTACATTCTCGAAGCTATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGTTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAg           <  1:2929741/150‑1 (MQ=255)
                                                                                                                                             cAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg        >  1:1351708/1‑150 (MQ=255)
                                                                                                                                             aataaaataCAGGAAATTCTCGAAGCGATTCGTTATCACCAGGTGGTGATCGTCGCCGGGTAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg        <  1:4169470/147‑1 (MQ=255)
                                                                                                                                              aGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgc       >  2:533196/1‑150 (MQ=255)
                                                                                                                                                aaaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt     >  1:3615468/1‑150 (MQ=255)
                                                                                                                                                 aaataCAGGACATTCTCGAAGCGATACGTGATCACCTGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc    <  1:4898128/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTCGT  >  CP009273/1477417‑1477713

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