Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,565 A→T 100% 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,5650AT95.4% 931.7 / 21.8 294K83I (AAA→ATA) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (4/10);  new base T (128/152);  total (132/162)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 305 total aligned reads displayed.

TGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTC  >  CP009273/1477418‑1477711
                                                                                                                                                   |                                                                                                                                                  
 gaagaaGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACa                                                                                                                                                 >  1:1991235/1‑150 (MQ=255)
  aagaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAg                                                                                                                                                >  2:1290495/1‑150 (MQ=255)
   agaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAgg                                                                                                                                               >  1:4950919/1‑150 (MQ=255)
   agaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAgg                                                                                                                                               <  1:1617291/150‑1 (MQ=255)
    gaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGa                                                                                                                                              >  1:169754/1‑150 (MQ=255)
      agGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACa                                                                                                                                            >  1:3562205/1‑150 (MQ=255)
      agGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACa                                                                                                                                            <  2:3504516/150‑1 (MQ=255)
       ggTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACAt                                                                                                                                           <  2:3394635/150‑1 (MQ=255)
       ggTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACAt                                                                                                                                           <  2:921934/150‑1 (MQ=255)
        gTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGTTTAGTCAGAAAATACAGGACAtt                                                                                                                                          <  2:2754991/150‑1 (MQ=255)
         ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtc                                                                                                                                         <  1:2539577/150‑1 (MQ=255)
         ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtc                                                                                                                                         >  2:2083681/1‑150 (MQ=255)
          tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtct                                                                                                                                        >  2:3143252/1‑150 (MQ=255)
          tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtct                                                                                                                                        <  1:954114/150‑1 (MQ=255)
          tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtct                                                                                                                                        <  1:5128148/150‑1 (MQ=255)
           aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtctc                                                                                                                                       <  2:2919063/150‑1 (MQ=255)
           aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtctc                                                                                                                                       >  2:2494972/1‑150 (MQ=255)
             aaaTCCTGATTCACAACAGGCCATTTTCCAGGAGATGTCGTAAGAGATTGACCAGGAGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCTGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGa                                                                                                                                     <  2:670242/150‑1 (MQ=255)
             aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGa                                                                                                                                     <  2:4551157/150‑1 (MQ=255)
              aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACACTTTACCGGTTAGTCAGAAAATACAGGACATTCTCGaa                                                                                                                                    >  1:2761391/1‑150 (MQ=255)
              aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGaa                                                                                                                                    <  1:539535/150‑1 (MQ=255)
               aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAg                                                                                                                                   <  2:97321/150‑1 (MQ=255)
               aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                                                                                                                                   <  2:1167004/150‑1 (MQ=255)
                tCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGc                                                                                                                                  <  1:3749891/150‑1 (MQ=255)
                 ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCg                                                                                                                                 >  1:414090/1‑150 (MQ=255)
                   tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGAt                                                                                                                               <  1:4248944/150‑1 (MQ=255)
                   tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGAt                                                                                                                               >  2:516001/1‑150 (MQ=255)
                     aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTc                                                                                                                             <  1:4346735/150‑1 (MQ=255)
                        cacaACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTg                                                                                                                          <  2:1860566/150‑1 (MQ=255)
                        cacaACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATGTACCGGTTAGTCAGAAAATACAGGACAGTCTCGAAGCGAGTCGTg                                                                                                                          >  1:3420499/1‑150 (MQ=255)
                         acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACACTTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGa                                                                                                                         >  1:643895/1‑150 (MQ=255)
                         acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGa                                                                                                                         >  2:1069106/1‑150 (MQ=255)
                           aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGATAAAACAGGACATTCTCGAAGCGATTCGTGATc                                                                                                                       <  1:4871555/150‑1 (MQ=255)
                            acaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCa                                                                                                                      <  1:2937623/150‑1 (MQ=255)
                              aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCAcc                                                                                                                    <  2:1025030/150‑1 (MQ=255)
                               ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   <  2:2832653/150‑1 (MQ=255)
                               ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   <  1:708579/150‑1 (MQ=255)
                               ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   <  1:3293673/150‑1 (MQ=255)
                               ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   <  2:758382/150‑1 (MQ=255)
                               ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   <  2:4432350/150‑1 (MQ=255)
                                 ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAgg                                                                                                                 >  1:1604133/1‑150 (MQ=255)
                                 ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAgg                                                                                                                 <  1:2710330/150‑1 (MQ=255)
                                 ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAgg                                                                                                                 >  2:2607511/1‑150 (MQ=255)
                                     tttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt                                                                                                             <  1:1044701/150‑1 (MQ=255)
                                         cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                         <  2:1455302/150‑1 (MQ=255)
                                         cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtc                                                                                                         >  2:5043610/1‑150 (MQ=255)
                                           ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgt                                                                                                       <  2:4334963/150‑1 (MQ=255)
                                             agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                     >  2:3934970/1‑150 (MQ=255)
                                             agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                     >  2:4710016/1‑150 (MQ=255)
                                                   cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgggg                                                                                               >  2:2565040/1‑150 (MQ=255)
                                                    gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                              <  2:2275177/150‑1 (MQ=255)
                                                     aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaa                                                                                             >  2:1574460/1‑150 (MQ=255)
                                                       agagATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAc                                                                                           <  1:2695578/150‑1 (MQ=255)
                                                        gagaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACg                                                                                          <  1:3860186/150‑1 (MQ=255)
                                                         agaTTGACCCGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACgg                                                                                         <  1:3897576/150‑1 (MQ=255)
                                                           aTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGt                                                                                       <  1:1560300/150‑1 (MQ=255)
                                                             tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                     >  2:2231004/1‑150 (MQ=255)
                                                             tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTc                                                                                     <  1:818851/150‑1 (MQ=255)
                                                               aCCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTg                                                                                   <  1:4724194/150‑1 (MQ=255)
                                                                ccAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                  <  1:3054508/150‑1 (MQ=255)
                                                                   ggcggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaa                                                                               >  2:3477726/1‑150 (MQ=255)
                                                                      ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTCTCCTGACAATGTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAAc                                                                            <  2:2518934/150‑1 (MQ=255)
                                                                      ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAAc                                                                            >  2:3565835/1‑150 (MQ=255)
                                                                         aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGAc                                                                         >  2:4909213/1‑150 (MQ=255)
                                                                          gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACt                                                                        >  1:1192380/1‑150 (MQ=255)
                                                                              aaGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGt                                                                    >  1:1788357/1‑150 (MQ=255)
                                                                                 tCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTAc                                                                 <  2:1883585/150‑1 (MQ=255)
                                                                                       tgcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGaaaa                                                           <  2:2418996/150‑1 (MQ=255)
                                                                                           tGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTg                                                       <  2:2012176/150‑1 (MQ=255)
                                                                                              aGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTAt                                                    >  2:2111188/1‑150 (MQ=255)
                                                                                                    aCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCt                                              >  2:3726083/1‑150 (MQ=255)
                                                                                                      gACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                                            >  2:148964/1‑150 (MQ=255)
                                                                                                       aCCGGAAATTACTTATCCTGACAATTTCCCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTggg                                           <  2:1979753/150‑1 (MQ=255)
                                                                                                        ccGTAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                          <  1:2730831/150‑1 (MQ=255)
                                                                                                        ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                          >  1:1936475/1‑150 (MQ=255)
                                                                                                         cGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGc                                               <  1:3221047/144‑1 (MQ=255)
                                                                                                         cGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgc                                         <  2:1092088/150‑1 (MQ=255)
                                                                                                           gAAATTAATTATCCTGAAAATTTACCGGTTAGTCAGAAAATACAGGAATTTCTCGAAGCGATTCGTGATCACCTGGTGGTGTTCGTCGCCGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                       <  1:2632604/150‑1 (MQ=255)
                                                                                                               ttACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGa                                   <  1:1189878/150‑1 (MQ=255)
                                                                                                                 aCTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAtt                                 <  2:1604133/150‑1 (MQ=255)
                                                                                                                 aCTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAtt                                 <  2:788037/150‑1 (MQ=255)
                                                                                                                       ccTGACACTTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGa                           >  2:4640052/1‑150 (MQ=255)
                                                                                                                       ccTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGa                           >  1:3661255/1‑150 (MQ=255)
                                                                                                                        cTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGAc                          >  2:411684/1‑150 (MQ=255)
                                                                                                                         tGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACt                         <  1:1290495/150‑1 (MQ=255)
                                                                                                                           aCAATTTACCGTTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGa                       <  1:313532/150‑1 (MQ=255)
                                                                                                                            cAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                      >  2:4243783/1‑150 (MQ=255)
                                                                                                                             aaTTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATc                     <  1:4696186/150‑1 (MQ=255)
                                                                                                                               tttACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCgg                   <  2:4766253/150‑1 (MQ=255)
                                                                                                                                   ccGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt               <  1:2089176/150‑1 (MQ=255)
                                                                                                                                   ccGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt               >  1:3707719/1‑150 (MQ=255)
                                                                                                                                     ggTTAGACATAAAATACTGGACATTCTCGACGCGATTCTTGATCACCTGGTGGTGATCGTCGCGGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATAc             <  1:606215/150‑1 (MQ=255)
                                                                                                                                         aGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAg         <  1:4240991/150‑1 (MQ=255)
                                                                                                                                            cAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcg      >  1:282079/1‑150 (MQ=255)
                                                                                                                                               aaaataCCGGACTTTCTCGACGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt   <  1:134728/150‑1 (MQ=255)
                                                                                                                                               aaaataCAGGACATTCTCGACGCTTTTCGTGATCACCTGGTGGTGATCGTCGCCGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt   <  1:3143252/150‑1 (MQ=255)
                                                                                                                                               aaaataCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt   >  1:4431266/1‑150 (MQ=255)
                                                                                                                                                aaataCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc  >  1:5065570/1‑150 (MQ=255)
                                                                                                                                                aaataCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc  >  2:3867750/1‑150 (MQ=255)
                                                                                                                                                aaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgtc  >  1:4544481/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTC  >  CP009273/1477418‑1477711

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: ‑