Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,565 A→T 36.4% 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,5650AT36.4% 203.6 / 192.4 186K83I (AAA→ATA) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (56/62);  new base T (36/32);  total (92/94)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

TGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTC  >  CP009273/1477418‑1477711
                                                                                                                                                   |                                                                                                                                                  
tGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaaa                                                                                                                                                    <  1:104066/148‑1 (MQ=255)
tGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaaa                                                                                                                                                    <  2:2040342/148‑1 (MQ=255)
 gaagaaGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaaaa                                                                                                                                                   >  1:2107433/1‑148 (MQ=255)
  aagaagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATAc                                                                                                                                                  <  2:576285/148‑1 (MQ=255)
     aagGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAgg                                                                                                                                               >  2:1313485/1‑148 (MQ=255)
      agGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGa                                                                                                                                              <  2:2424742/148‑1 (MQ=255)
         ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACAt                                                                                                                                           >  1:423818/1‑148 (MQ=255)
         ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAt                                                                                                                                           <  1:194569/148‑1 (MQ=255)
         ttAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAt                                                                                                                                           >  1:2499888/1‑148 (MQ=255)
          tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACAtt                                                                                                                                          >  1:1352383/1‑148 (MQ=255)
          tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAtt                                                                                                                                          <  1:2100500/148‑1 (MQ=255)
           aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATtc                                                                                                                                         <  2:1127007/148‑1 (MQ=255)
           aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                                                                                                                         <  1:3305255/148‑1 (MQ=255)
            aaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtct                                                                                                                                        <  2:3250668/148‑1 (MQ=255)
             aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGACCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtctc                                                                                                                                       >  1:2630299/1‑148 (MQ=255)
              aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCg                                                                                                                                      >  2:1215064/1‑148 (MQ=255)
              aaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCg                                                                                                                                      >  2:2667038/1‑148 (MQ=255)
               aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGa                                                                                                                                     >  1:1877953/1‑148 (MQ=255)
                 ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGACGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGacg                                                                                                                                   >  2:1794824/1‑146 (MQ=255)
                 ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                                                                                                                                   <  2:1775170/148‑1 (MQ=255)
                 ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                                                                                                                                   >  2:1794974/1‑148 (MQ=255)
                   tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCg                                                                                                                                 >  1:1957380/1‑148 (MQ=255)
                    gATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGa                                                                                                                                >  2:3153124/1‑148 (MQ=255)
                     aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                                                                                                                               >  2:2875581/1‑148 (MQ=255)
                      tGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGAtt                                                                                                                              >  2:910732/1‑148 (MQ=255)
                       gCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTc                                                                                                                             >  2:2835441/1‑148 (MQ=255)
                        cacaACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCg                                                                                                                            <  1:230255/148‑1 (MQ=255)
                          caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTg                                                                                                                          >  1:2201402/1‑148 (MQ=255)
                           aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGa                                                                                                                         <  2:1360516/148‑1 (MQ=255)
                            acaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAt                                                                                                                        <  1:2369576/148‑1 (MQ=255)
                             caGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATc                                                                                                                       <  2:509317/148‑1 (MQ=255)
                               ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAc                                                                                                                     <  2:1880700/148‑1 (MQ=255)
                                 ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   <  1:1321519/148‑1 (MQ=255)
                                 ccATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCa                                                                                                                   >  2:279475/1‑148 (MQ=255)
                                       tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt                                                                                                             <  1:821599/148‑1 (MQ=255)
                                       tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt                                                                                                             <  2:2309038/148‑1 (MQ=255)
                                         cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                           <  1:1470038/148‑1 (MQ=255)
                                         cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGa                                                                                                           <  1:1467247/148‑1 (MQ=255)
                                               aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                     <  1:2043706/148‑1 (MQ=255)
                                               aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAtcgtcg                                                                                                     <  1:649608/148‑1 (MQ=255)
                                                 ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTAATCCTGACAATTTACCGGTAAGACAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGcc                                                                                                   >  1:2758922/1‑148 (MQ=255)
                                                   cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCgg                                                                                                 >  1:1220503/1‑148 (MQ=255)
                                                    gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCggg                                                                                                <  1:433420/148‑1 (MQ=255)
                                                      aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGa                                                                                              >  1:2955324/1‑148 (MQ=255)
                                                        gagaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGaaa                                                                                            <  2:1216697/148‑1 (MQ=255)
                                                         agaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAc                                                                                           >  1:1367193/1‑148 (MQ=255)
                                                          gaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACg                                                                                          <  1:916857/148‑1 (MQ=255)
                                                           aTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACgg                                                                                         <  2:1948139/148‑1 (MQ=255)
                                                             tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGt                                                                                       >  1:1532965/1‑148 (MQ=255)
                                                              taCCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                      <  2:2848806/147‑1 (MQ=255)
                                                              gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                      <  1:2744372/148‑1 (MQ=255)
                                                              gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGtt                                                                                      <  1:1311636/148‑1 (MQ=255)
                                                                 ccggcggcCGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTg                                                                                   <  1:3219762/146‑1 (MQ=255)
                                                                  aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                  >  2:1439974/1‑148 (MQ=255)
                                                                  aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTgg                                                                                  >  2:3082746/1‑148 (MQ=255)
                                                                      ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaaa                                                                              <  1:1840747/148‑1 (MQ=255)
                                                                      ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTaaa                                                                              <  2:1171523/148‑1 (MQ=255)
                                                                         aGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACg                                                                           <  1:732300/148‑1 (MQ=255)
                                                                          gggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGa                                                                          >  1:2845878/1‑148 (MQ=255)
                                                                           ggAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGAc                                                                         >  1:2610543/1‑148 (MQ=255)
                                                                                   ctgctgCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGACGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAAGGGGTTCTGGTAAAACGACTCAGTTAc                                                                 <  2:2563252/148‑1 (MQ=255)
                                                                                     gctgcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCg                                                               >  1:1674952/1‑148 (MQ=255)
                                                                                     gctgcGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCg                                                               <  2:1653704/148‑1 (MQ=255)
                                                                                        gcGTTAAGCGGCACGACGGTAAATTACTTATCCTGACATTTTACCGTTTAGTCATAAAATACAGGAAATTCTCGACGCGATTCGTGATCACCAGGTGTTTATGGTCGCCGGGGAAACGGGTTGTGGTAAAACGACTCAGTTACCGaaa                                                            <  2:649714/148‑1 (MQ=255)
                                                                                          gTGAAGCGGCCCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAAt                                                          <  2:2731871/148‑1 (MQ=255)
                                                                                           tGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATc                                                         >  2:125239/1‑148 (MQ=255)
                                                                                           tGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATc                                                         >  2:1022917/1‑148 (MQ=255)
                                                                                           tGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATGCTCGAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATc                                                         <  2:1626912/148‑1 (MQ=255)
                                                                                               gCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTa                                                     >  1:1270800/1‑148 (MQ=255)
                                                                                                   cACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGACGCGATTCGTGATCACCCGGTGTTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGa                                                 <  2:730780/148‑1 (MQ=255)
                                                                                                   cACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGa                                                 >  2:1673421/1‑148 (MQ=255)
                                                                                                     cGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGc                                               >  1:1580606/1‑148 (MQ=255)
                                                                                                        ccGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgg                                            >  2:3292515/1‑148 (MQ=255)
                                                                                                          ggAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTgggg                                          >  1:919539/1‑148 (MQ=255)
                                                                                                            aaaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcg                                        >  1:1752802/1‑148 (MQ=255)
                                                                                                             aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                       >  1:2453293/1‑148 (MQ=255)
                                                                                                             aaTTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGgcgc                                       >  2:2896325/1‑148 (MQ=255)
                                                                                                                tACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCggg                                    >  2:1684174/1‑148 (MQ=255)
                                                                                                                 aCTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGa                                   >  2:2310171/1‑148 (MQ=255)
                                                                                                                  cTTATCCTGACAATTTCCCGGTTAGTCAGAAAAAACAGGACATTCTCTAAGCGATTCGTTATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAt                                  <  2:2921293/148‑1 (MQ=255)
                                                                                                                  cTTATCCTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGCTTCTGGTACAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGAt                                  >  2:3197170/1‑148 (MQ=255)
                                                                                                                     aTCCTGACAATTTACCGGTTAGTCATAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTaa                               <  1:547300/148‑1 (MQ=255)
                                                                                                                       ccTGACAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAg                             >  1:1807734/1‑148 (MQ=255)
                                                                                                                           aCAATTTACCGGTTAGTCAGAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACt                         >  2:2887226/1‑148 (MQ=255)
                                                                                                                           aCAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACt                         >  2:1310906/1‑148 (MQ=255)
                                                                                                                            cAATTTCCCGGTTAGTCACAAAATACCGGAAATTCTCTACGCTTTTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAACCGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg                        <  1:996358/148‑1 (MQ=255)
                                                                                                                            cAATTTACCGTTTAGCAACAAAATACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGGCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg                        <  2:2991013/148‑1 (MQ=255)
                                                                                                                            cAATTTACCGGTTAGTCAGAAAATACCGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTg                        <  2:674235/148‑1 (MQ=255)
                                                                                                                             aaTTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGa                       >  2:45508/1‑148 (MQ=255)
                                                                                                                             aaTTTACCGGTTAGACACAAAATACAGGACATTCTCGAAGCGATTCGTTATCACCTGGTGGTTATGGTCGCGGGGTAAACGGGTTCTGGAAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGa                       <  1:1173749/148‑1 (MQ=255)
                                                                                                                              aTTTACCGTTTAGTCATAAAAAACAGGACATTCTCTAAGCTATTCGTTATCACCCGGTGGTGATCGTCGCCGGGTAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                      <  1:279475/148‑1 (MQ=255)
                                                                                                                              aTTTACCGGTTATTCATAAAAAACAGGACTTTCTCGAAGCGATTCGTGATCACCCGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGAt                      <  1:278787/148‑1 (MQ=255)
                                                                                                                                ttACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCg                    >  1:986610/1‑148 (MQ=255)
                                                                                                                                   ccGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGcc                 >  2:2367440/1‑148 (MQ=255)
                                                                                                                                   ccGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGcc                 >  2:2231154/1‑148 (MQ=255)
                                                                                                                                     ggTTAGTCACAAAAAACAGGACATTCTCGACGCGATTCGTGATCACCCGGTGGGGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCAt               <  1:1363223/148‑1 (MQ=255)
                                                                                                                                            cacaaaataCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTTATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGc        <  1:1313485/145‑1 (MQ=255)
                                                                                                                                               aaaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgc     >  1:740406/1‑148 (MQ=255)
                                                                                                                                               aaaaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCcgc     >  1:739735/1‑148 (MQ=255)
                                                                                                                                                 aaaaCAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGcgt   >  1:1450829/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAAGAAGGTTAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGATCGTCGCCGGGGAAACGGGTTCTGGTAAAACGACTCAGTTACCGAAAATCTGTATGGAGCTGGGGCGCGGGATTAAAGGACTGATCGGCCATACCCAGCCGCGTC  >  CP009273/1477418‑1477711

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: