Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,478,971 G→A 100% E552K (GAG→AAG)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,478,9710GA98.6% 711.3 / 0.1 220E552K (GAG→AAG) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (116/100);  total (118/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 223 total aligned reads displayed.

CGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCC  >  CP009273/1478827‑1479115
                                                                                                                                                |                                                                                                                                                
cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                               <  1:463644/148‑1 (MQ=255)
cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                               >  1:1205893/1‑148 (MQ=255)
cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                               >  1:2578126/1‑148 (MQ=255)
  tATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCt                                                                                                                                             >  2:1230322/1‑148 (MQ=255)
  tATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCt                                                                                                                                             >  2:1230632/1‑148 (MQ=255)
   aTGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTg                                                                                                                                            <  2:2779955/148‑1 (MQ=255)
   aTGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTg                                                                                                                                            <  2:1589009/148‑1 (MQ=255)
          tGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtct                                                                                                                                     <  2:1444967/148‑1 (MQ=255)
              aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                 <  1:1051411/148‑1 (MQ=255)
              aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                 <  1:1876047/148‑1 (MQ=255)
              aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                 <  2:2117738/148‑1 (MQ=255)
               gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                                <  2:2054581/148‑1 (MQ=255)
                    aaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtgtg                                                                                                                           <  1:2640504/148‑1 (MQ=255)
                     aaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGa                                                                                                                          <  2:1869965/148‑1 (MQ=255)
                     aaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGa                                                                                                                          <  1:1933733/148‑1 (MQ=255)
                       aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAAt                                                                                                                        <  1:1228833/148‑1 (MQ=255)
                       aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAAt                                                                                                                        >  1:267605/1‑148 (MQ=255)
                        cATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATc                                                                                                                       <  2:1983257/148‑1 (MQ=255)
                          tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtg                                                                                                                     <  1:877907/148‑1 (MQ=255)
                          tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtg                                                                                                                     <  1:2634765/148‑1 (MQ=255)
                           ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgt                                                                                                                    >  1:3462183/1‑148 (MQ=255)
                           ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgt                                                                                                                    <  1:3263366/148‑1 (MQ=255)
                           ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgt                                                                                                                    >  1:2956580/1‑148 (MQ=255)
                            gCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgtg                                                                                                                   <  2:906820/148‑1 (MQ=255)
                               gcgtgcgtGAGGCGATGATTCTCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGCTGGACAAACAGCAGGCATCGTACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTGTCTCGCGTTTGTGAATCTGTGGaa                                                                                                                <  2:3200808/148‑1 (MQ=255)
                               gcgtgcgtGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGaa                                                                                                                <  1:470602/148‑1 (MQ=255)
                                   gcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTAt                                                                                                            <  1:874841/148‑1 (MQ=255)
                                    cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATc                                                                                                           <  2:2565649/148‑1 (MQ=255)
                                    cgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATc                                                                                                           <  1:1082622/148‑1 (MQ=255)
                                     gtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCt                                                                                                          <  2:2338179/148‑1 (MQ=255)
                                      tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                         >  2:1325544/1‑148 (MQ=255)
                                      tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                         >  2:3155312/1‑148 (MQ=255)
                                      tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCtt                                                                                                         <  1:2233945/148‑1 (MQ=255)
                                       gAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTg                                                                                                        <  2:1150672/148‑1 (MQ=255)
                                        aGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTgg                                                                                                       >  2:574154/1‑148 (MQ=255)
                                        aGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTgg                                                                                                       >  2:323670/1‑148 (MQ=255)
                                        aGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACAACGACGCTTCCACGACAAAAAGACTGACTTTCTCGCGATTGTGAATCTGTGGAATTATCTTgg                                                                                                       >  2:2433226/1‑148 (MQ=255)
                                          gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCg                                                                                                     >  2:1036323/1‑148 (MQ=255)
                                             atgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGCATTATCTTGGCGagc                                                                                                  >  1:2133515/1‑148 (MQ=255)
                                              tgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                                                                                 >  2:2902897/1‑148 (MQ=255)
                                                      aCGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAgg                                                                                         <  1:591535/148‑1 (MQ=255)
                                                          ccGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCt                                                                                     >  2:1695765/1‑148 (MQ=255)
                                                           cgcgCTCTCCATTCTGGATCCGCGCGAACGTCCGTTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGGGTTTGTGATTCTGTGGAATTATCTTGGCGAGCAGCAAAATGCGCtt                                                                                    <  2:3466802/148‑1 (MQ=255)
                                                           cgcgCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCtt                                                                                    <  2:1100706/148‑1 (MQ=255)
                                                             cgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttc                                                                                  <  1:1901120/148‑1 (MQ=255)
                                                             cgcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttc                                                                                  >  1:2834800/1‑148 (MQ=255)
                                                                tctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttcttc                                                                               >  1:968692/1‑148 (MQ=255)
                                                                 ctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTcc                                                                              >  2:826439/1‑148 (MQ=255)
                                                                 ctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTcc                                                                              >  1:1015216/1‑148 (MQ=255)
                                                                  tcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCa                                                                             >  1:2315628/1‑148 (MQ=255)
                                                                   ccATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCaa                                                                            <  2:928586/148‑1 (MQ=255)
                                                                    cATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAAc                                                                           >  2:3059634/1‑148 (MQ=255)
                                                                     aTTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGACTTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACg                                                                          >  2:631236/1‑148 (MQ=255)
                                                                     aTTCAGGATCCGCACGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACg                                                                          <  2:773398/148‑1 (MQ=255)
                                                                     aTTCAGGATCCGCACGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACg                                                                          <  2:1466400/148‑1 (MQ=255)
                                                                       tCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGcc                                                                        >  2:3077924/1‑148 (MQ=255)
                                                                        cAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCt                                                                       >  1:488045/1‑148 (MQ=255)
                                                                         aGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCtt                                                                      <  2:2578994/148‑1 (MQ=255)
                                                                          ggATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTc                                                                     >  1:1543287/1‑148 (MQ=255)
                                                                          ggATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTc                                                                     >  1:1541781/1‑148 (MQ=255)
                                                                           gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGACTTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTcc                                                                    >  2:951369/1‑148 (MQ=255)
                                                                              ccGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGAATTTCACGCGTTTGTGAATCTGTGGAATTAACTTTGCGAGCAGCAAAATGCGCTTTCTTCCAACGCCTTCCGTc                                                                 >  1:846857/1‑148 (MQ=255)
                                                                                 cgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGcc                                                              >  1:2317538/1‑148 (MQ=255)
                                                                                 cgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGcc                                                              >  2:1797820/1‑148 (MQ=255)
                                                                                         tCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGt                                                      >  1:3011729/1‑148 (MQ=255)
                                                                                          ccGATGGACAAACAGCAGGCATCGGACGAAAAACAGCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATATTGGCGAGCAGCAAAAGGCGCTTTCGTCCAACGCCTTCCGTCGCCTGTGTCGTa                                                     <  2:2569521/148‑1 (MQ=255)
                                                                                             aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTTAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCg                                                  <  2:1692678/148‑1 (MQ=255)
                                                                                              tGGACAAACAGCAGGCATCGGACGAAACACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGa                                                 <  2:3399027/148‑1 (MQ=255)
                                                                                                gACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAtt                                               >  2:3415294/1‑148 (MQ=255)
                                                                                                   aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtc                                            >  2:1753923/1‑148 (MQ=255)
                                                                                                     aCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtctc                                          >  2:1237918/1‑148 (MQ=255)
                                                                                                      cagcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCa                                         <  1:3169165/148‑1 (MQ=255)
                                                                                                      cagcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCa                                         >  2:2804419/1‑148 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      <  1:507264/148‑1 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      >  2:471373/1‑148 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      <  1:797457/148‑1 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      >  1:2262962/1‑148 (MQ=255)
                                                                                                           ggCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTACAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTAt                                    >  1:2655530/1‑148 (MQ=255)
                                                                                                            gCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATc                                   <  2:479449/148‑1 (MQ=255)
                                                                                                            gCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATc                                   <  1:3077924/148‑1 (MQ=255)
                                                                                                              aTCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtg                                 >  1:1778252/1‑148 (MQ=255)
                                                                                                                  gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcg                             >  1:3355234/1‑148 (MQ=255)
                                                                                                                  gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcg                             >  1:171230/1‑148 (MQ=255)
                                                                                                                  gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcg                             <  2:125008/148‑1 (MQ=255)
                                                                                                                  gACGAAAAACATCGTCGCTTACACGACAAAAAGTCTGACTTTCTCGGGTTTGGGAAACTGTGGAATTATAGTGGCGAGCAGCAAAACGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTGGTACCGATTATCTAAACTAACtgcgcg                             >  1:3287967/1‑148 (MQ=255)
                                                                                                                   aCGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgt                            <  2:3225642/148‑1 (MQ=255)
                                                                                                                      aaaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGGTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcg                         >  2:3296673/1‑148 (MQ=255)
                                                                                                                         aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                      <  2:1283115/148‑1 (MQ=255)
                                                                                                                         aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                      <  2:899269/148‑1 (MQ=255)
                                                                                                                         aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                      <  1:1978503/148‑1 (MQ=255)
                                                                                                                         aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                      <  1:486421/148‑1 (MQ=255)
                                                                                                                          aCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGaa                     <  2:267605/148‑1 (MQ=255)
                                                                                                                              cgtcgCTTCCACGACAAAAAGTCTGACTTTATCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTATTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATTGc                 >  2:2844844/1‑148 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg              >  1:215186/1‑148 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg              >  1:1400357/1‑148 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg              >  2:2411568/1‑148 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg              <  2:1205893/148‑1 (MQ=255)
                                                                                                                                            cAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAcc   >  1:2585570/1‑148 (MQ=255)
                                                                                                                                            cAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAcc   >  1:769760/1‑148 (MQ=255)
                                                                                                                                             aaaaAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAccc  >  2:1137727/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCC  >  CP009273/1478827‑1479115

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