Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,495,859 C→A 6.8% E702D (GAG→GAT fryA ← putative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,495,8590CA 6.8% 660.5 / 16.9 176E702D (GAG→GATfryAputative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (72/92);  new base A (3/9);  total (75/101)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 181 total aligned reads displayed.

CTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACA  >  CP009273/2495712‑2496005
                                                                                                                                                   |                                                                                                                                                  
cTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGc                                                                                                                                                    >  2:1155284/1‑148 (MQ=255)
 tGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGct                                                                                                                                                   <  2:64951/148‑1 (MQ=255)
     aTTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCttt                                                                                                                                               <  2:924546/148‑1 (MQ=255)
     aTTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCttt                                                                                                                                               >  2:2389641/1‑148 (MQ=255)
       ttGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTAt                                                                                                                                             <  2:3297434/148‑1 (MQ=255)
         ggTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTg                                                                                                                                           <  2:15951/148‑1 (MQ=255)
         ggTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTg                                                                                                                                           >  1:401873/1‑148 (MQ=255)
           tgtgCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGACGCACAGGAACTGGATCGCCTGCTCTTTATTGCt                                                                                                                                         >  1:2507487/1‑148 (MQ=255)
             tgCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTaa                                                                                                                                       <  1:895944/148‑1 (MQ=255)
              gCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTaaa                                                                                                                                      <  1:2968944/148‑1 (MQ=255)
              gCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTaaa                                                                                                                                      <  2:1887756/148‑1 (MQ=255)
                 ggATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATc                                                                                                                                   >  2:1295154/1‑148 (MQ=255)
                 ggATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATc                                                                                                                                   <  1:1831351/148‑1 (MQ=255)
                 ggATCGCTACGCAAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATc                                                                                                                                   <  1:3080321/148‑1 (MQ=255)
                    tCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTg                                                                                                                                <  2:3058443/148‑1 (MQ=255)
                    tCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTg                                                                                                                                <  2:1420748/148‑1 (MQ=255)
                     cGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGa                                                                                                                               <  1:587879/148‑1 (MQ=255)
                          cGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaac                                                                                                                          >  2:3101732/1‑148 (MQ=255)
                          cGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaac                                                                                                                          >  1:3285796/1‑148 (MQ=255)
                          cGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATcaac                                                                                                                          >  1:1831627/1‑148 (MQ=255)
                                 aCCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGata                                                                                                                   >  2:2323170/1‑148 (MQ=255)
                                 aCCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGata                                                                                                                   <  2:2686354/148‑1 (MQ=255)
                                  ccAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGatat                                                                                                                  <  1:1374124/148‑1 (MQ=255)
                                   cAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATAtt                                                                                                                 <  1:2449451/148‑1 (MQ=255)
                                     aCGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATtct                                                                                                               >  1:3219483/1‑148 (MQ=255)
                                      cGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATtctc                                                                                                              >  1:2748488/1‑148 (MQ=255)
                                       gCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTcc                                                                                                             <  1:2621431/148‑1 (MQ=255)
                                         cGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGNCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAg                                                                                                           <  1:1679955/148‑1 (MQ=255)
                                             ggtAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGcc                                                                                                       <  2:1316832/148‑1 (MQ=255)
                                                 aCAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCcagca                                                                                                   <  1:1295154/148‑1 (MQ=255)
                                                 aCAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCcagca                                                                                                   <  1:1546475/148‑1 (MQ=255)
                                                 aCAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCcagca                                                                                                   <  1:2942743/148‑1 (MQ=255)
                                                  cAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTATCCAGCGCcagcag                                                                                                  >  2:2192496/1‑148 (MQ=255)
                                                   aaTCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGt                                                                                                 >  2:1735119/1‑148 (MQ=255)
                                                   aaTCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGt                                                                                                 <  2:3104948/148‑1 (MQ=255)
                                                    aTCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTg                                                                                                <  2:2434376/148‑1 (MQ=255)
                                                     tctcTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTgg                                                                                               <  1:433587/148‑1 (MQ=255)
                                                        ctTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCg                                                                                            >  2:3004066/1‑148 (MQ=255)
                                                          tCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGaa                                                                                          >  2:2927872/1‑148 (MQ=255)
                                                          tCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGaa                                                                                          >  1:1942670/1‑148 (MQ=255)
                                                           cccGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAc                                                                                         >  1:2718945/1‑148 (MQ=255)
                                                              gCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtc                                                                                      <  2:2409202/148‑1 (MQ=255)
                                                              gCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGtc                                                                                      <  2:1668339/148‑1 (MQ=255)
                                                                tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtctt                                                                                    <  2:1469081/148‑1 (MQ=255)
                                                                tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtctt                                                                                    <  2:1469014/148‑1 (MQ=255)
                                                                  ccAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttct                                                                                  <  2:2268213/148‑1 (MQ=255)
                                                                   cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTcttctt                                                                                 <  2:675089/148‑1 (MQ=255)
                                                                   cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTcttctt                                                                                 <  1:792187/148‑1 (MQ=255)
                                                                       cATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCgg                                                                             >  2:1229671/1‑148 (MQ=255)
                                                                          cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGt                                                                          >  2:1184026/1‑148 (MQ=255)
                                                                          cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGt                                                                          <  1:1911597/148‑1 (MQ=255)
                                                                           ttcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGGTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTa                                                                         <  1:596063/148‑1 (MQ=255)
                                                                            tcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                        <  1:954176/148‑1 (MQ=255)
                                                                            tcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                        <  2:2738850/148‑1 (MQ=255)
                                                                             cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                       <  2:567603/148‑1 (MQ=255)
                                                                             cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                       >  2:3286020/1‑148 (MQ=255)
                                                                             cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                       <  1:3119435/148‑1 (MQ=255)
                                                                             cttcCACCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                       <  1:572961/148‑1 (MQ=255)
                                                                                    cTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                >  2:2273526/1‑148 (MQ=255)
                                                                                    cTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTAGTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                <  2:485395/148‑1 (MQ=255)
                                                                                      cAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAg                                                              >  2:3217354/1‑148 (MQ=255)
                                                                                       aaaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGGACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGt                                                             >  1:1986440/1‑148 (MQ=255)
                                                                                       aaaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGt                                                             >  1:318954/1‑148 (MQ=255)
                                                                                        aaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTa                                                            >  1:1012385/1‑148 (MQ=255)
                                                                                         aCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTaa                                                           <  1:3062986/148‑1 (MQ=255)
                                                                                           ggATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACg                                                         <  2:458972/148‑1 (MQ=255)
                                                                                            gATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGc                                                        <  2:1358213/148‑1 (MQ=255)
                                                                                             aTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCt                                                       >  2:632399/1‑148 (MQ=255)
                                                                                             aTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCt                                                       >  2:2129469/1‑148 (MQ=255)
                                                                                              tGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCtt                                                      <  1:1610880/148‑1 (MQ=255)
                                                                                              tGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCtt                                                      <  2:1675145/148‑1 (MQ=255)
                                                                                               gTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTc                                                     >  2:2672557/1‑148 (MQ=255)
                                                                                                ttCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCa                                                    <  1:526076/148‑1 (MQ=255)
                                                                                                 tCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCaa                                                   <  1:3029392/148‑1 (MQ=255)
                                                                                                 tCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCaa                                                   >  2:3303737/1‑148 (MQ=255)
                                                                                                 tCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCGTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCaa                                                   >  1:1332611/1‑148 (MQ=255)
                                                                                                     tGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTc                                               <  2:314633/148‑1 (MQ=255)
                                                                                                        gcCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTg                                            <  2:1909234/148‑1 (MQ=255)
                                                                                                          ccGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTgcg                                          >  1:930635/1‑148 (MQ=255)
                                                                                                            gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCa                                        <  1:2848306/148‑1 (MQ=255)
                                                                                                            gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCa                                        >  2:2063051/1‑148 (MQ=255)
                                                                                                                         ttGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcac                           <  1:3101732/148‑1 (MQ=255)
                                                                                                                          tGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcaca                          <  1:330485/148‑1 (MQ=255)
                                                                                                                             cGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGc                       <  1:207007/148‑1 (MQ=255)
                                                                                                                             cGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGc                       <  1:1155284/148‑1 (MQ=255)
                                                                                                                             cGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGc                       <  1:2753632/148‑1 (MQ=255)
                                                                                                                              gCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGcc                      >  2:381670/1‑148 (MQ=255)
                                                                                                                                acaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTg                    <  1:1075990/148‑1 (MQ=255)
                                                                                                                                   ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                 >  2:3002704/1‑148 (MQ=255)
                                                                                                                                   ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                 <  2:401873/148‑1 (MQ=255)
                                                                                                                                   ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                 >  1:3017310/1‑148 (MQ=255)
                                                                                                                                   ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCAGTCACATGCCTGAcg                 >  1:3017298/1‑148 (MQ=255)
                                                                                                                                   ggAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                 <  1:222435/148‑1 (MQ=255)
                                                                                                                                    gAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcgc                >  2:2624080/1‑148 (MQ=255)
                                                                                                                                      aCTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACgcgc              <  1:1170047/148‑1 (MQ=255)
                                                                                                                                      aCTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACgcgc              <  2:2485223/148‑1 (MQ=255)
                                                                                                                                              gCCAGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGATCGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACTCGCCAGTTccc      <  2:1361256/148‑1 (MQ=255)
                                                                                                                                                cTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGa    >  2:1076550/1‑148 (MQ=255)
                                                                                                                                                cTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGa    >  2:1077821/1‑148 (MQ=255)
                                                                                                                                                  gCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCCGTTCCCGaca  >  2:2346224/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGAGATTTGGTGTGCGGGATCGCTACGCCAAAACCAACGCCGGTGGTAACAATCTCTTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACA  >  CP009273/2495712‑2496005

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: