Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,644,550 C→A 17.7% pseudogene (216/765 nt) ydfE → Qin prophage; pseudogene;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,644,5500CA17.7% 646.6 / 78.2 222pseudogene (216/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (96/86);  new base A (21/19);  total (117/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

GTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTC  >  CP009273/1644404‑1644697
                                                                                                                                                  |                                                                                                                                                   
gTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCg                                                                                                                                                    <  2:1943774/148‑1 (MQ=255)
 tAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGa                                                                                                                                                   >  1:2177843/1‑148 (MQ=255)
  aGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGAt                                                                                                                                                  <  1:2737547/148‑1 (MQ=255)
    ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTc                                                                                                                                                <  2:3142697/148‑1 (MQ=255)
    ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTc                                                                                                                                                >  1:2953870/1‑148 (MQ=255)
    ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTc                                                                                                                                                >  1:2953922/1‑148 (MQ=255)
    ccACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTc                                                                                                                                                <  2:1912191/148‑1 (MQ=255)
       caGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAGTTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATg                                                                                                                                             >  1:1491931/1‑148 (MQ=255)
       caGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATg                                                                                                                                             <  2:3618149/148‑1 (MQ=255)
       caGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATg                                                                                                                                             <  1:956118/148‑1 (MQ=255)
       caGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCAtt                                                                                                                                             <  1:1437570/148‑2 (MQ=255)
       caGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATg                                                                                                                                             <  1:3340159/148‑1 (MQ=255)
            aTTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGTACTGCGCTTATAAATATATTGTCTGCGAGGATAAACCTTGCATACCCCGCCCTGCCCTCGATTCATGGGATg                                                                                                                                        >  2:2604307/1‑148 (MQ=255)
            aTTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATg                                                                                                                                        >  2:1223659/1‑148 (MQ=255)
             ttGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGc                                                                                                                                       <  1:1436353/148‑1 (MQ=255)
             ttGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGc                                                                                                                                       <  1:2855899/148‑1 (MQ=255)
              tGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCt                                                                                                                                      <  1:3582909/148‑1 (MQ=255)
              tGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCt                                                                                                                                      >  2:1406356/1‑148 (MQ=255)
              tGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCt                                                                                                                                      >  2:600260/1‑148 (MQ=255)
              tGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCt                                                                                                                                      >  2:1070232/1‑148 (MQ=255)
              tGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCt                                                                                                                                      >  1:3035695/1‑148 (MQ=255)
                  gCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  <  1:85007/148‑1 (MQ=255)
                  gCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAat                                                                                                                                  >  2:3587211/1‑148 (MQ=255)
                   cTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCTCCCTGCCCTCGATTCATGGGATGCTGAata                                                                                                                                 >  1:2109396/1‑148 (MQ=255)
                     ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAatata                                                                                                                               <  1:237007/148‑1 (MQ=255)
                     ggCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAatata                                                                                                                               <  1:3140692/148‑1 (MQ=255)
                         acgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCa                                                                                                                           >  1:3370877/1‑148 (MQ=255)
                         acgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCa                                                                                                                           >  1:3327520/1‑148 (MQ=255)
                         acgaACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCa                                                                                                                           >  2:3218299/1‑148 (MQ=255)
                              cGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGaaa                                                                                                                      <  1:2533976/148‑1 (MQ=255)
                               gCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGaaaa                                                                                                                     <  1:1249160/148‑1 (MQ=255)
                                    gCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCt                                                                                                                <  1:437288/148‑1 (MQ=255)
                                    gCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCt                                                                                                                <  1:1024076/148‑1 (MQ=255)
                                     cTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTg                                                                                                               <  1:2660014/148‑1 (MQ=255)
                                      tCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTgg                                                                                                              <  2:2898684/148‑1 (MQ=255)
                                           cAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCAGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATga                                                                                                         >  2:221713/1‑148 (MQ=255)
                                               aTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATg                                                                                                          >  1:1999328/1‑143 (MQ=255)
                                                   cACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCATCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTg                                                                                                 <  1:1415318/148‑1 (MQ=255)
                                                    aCTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGc                                                                                                >  1:1859393/1‑148 (MQ=255)
                                                    aCTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGc                                                                                                <  2:3736735/148‑1 (MQ=255)
                                                      tGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGctt                                                                                              <  1:2140077/148‑1 (MQ=255)
                                                       ggCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttc                                                                                             <  2:1037764/148‑1 (MQ=255)
                                                       ggCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttc                                                                                             <  2:1604821/148‑1 (MQ=255)
                                                        gCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  1:1526841/148‑1 (MQ=255)
                                                        gCAAGTCCGGGAAGATTATACGGGTGCTCTGGTCTGCGGTGATAAATTTATTGTATGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATGCATGGTTTGATTAATTTATGCAGGAAAACCGCTGGGATGAGGAGTCTGcttct                                                                                            <  1:2604307/148‑1 (MQ=255)
                                                             tCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGt                                                                                       >  2:2159672/1‑148 (MQ=255)
                                                               cTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTcc                                                                                     <  2:838814/148‑1 (MQ=255)
                                                               cTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTcc                                                                                     <  1:321969/148‑1 (MQ=255)
                                                                  gAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCgg                                                                                  >  2:979697/1‑148 (MQ=255)
                                                                   aaCATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGt                                                                                 >  1:1209592/1‑148 (MQ=255)
                                                                       ttATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTgag                                                                             >  2:2599929/1‑148 (MQ=255)
                                                                        tATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTgaga                                                                            >  2:754557/1‑148 (MQ=255)
                                                                        tATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTgaga                                                                            >  2:754758/1‑148 (MQ=255)
                                                                           ccGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTg                                                                         >  1:3753646/1‑148 (MQ=255)
                                                                           ccGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTg                                                                         >  2:3431216/1‑148 (MQ=255)
                                                                                tgctCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatcc                                                                    >  1:3163342/1‑148 (MQ=255)
                                                                                tgctCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatcc                                                                    <  2:648074/148‑1 (MQ=255)
                                                                                   tCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccgat                                                                 <  1:2359278/148‑1 (MQ=255)
                                                                                    cAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccgatc                                                                >  1:1413694/1‑148 (MQ=255)
                                                                                      ggACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccgatccg                                                              >  1:416510/1‑148 (MQ=255)
                                                                                      ggACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccgatccg                                                              >  1:416420/1‑148 (MQ=255)
                                                                                        aCTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATccgatccgat                                                            >  2:1485409/1‑148 (MQ=255)
                                                                                          tGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGa                                                          >  1:1798478/1‑148 (MQ=255)
                                                                                               ttATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAGCCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTc                                                     <  1:2158018/148‑1 (MQ=255)
                                                                                                tataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCa                                                    <  2:1578773/148‑1 (MQ=255)
                                                                                                tataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCa                                                    >  2:3505247/1‑148 (MQ=255)
                                                                                                tataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCa                                                    >  1:1310216/1‑148 (MQ=255)
                                                                                                 ataAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCAc                                                   <  1:1679977/148‑1 (MQ=255)
                                                                                                  taAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACt                                                  >  2:2613328/1‑148 (MQ=255)
                                                                                                   aaaTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACtt                                                 <  2:2727965/148‑1 (MQ=255)
                                                                                                   aaaTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACtt                                                 >  1:2944139/1‑148 (MQ=255)
                                                                                                    aaTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACttt                                                <  2:2051996/148‑1 (MQ=255)
                                                                                                     aTTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACtttt                                               >  1:825370/1‑148 (MQ=255)
                                                                                                      tttATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGTTGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTg                                              <  2:305594/148‑1 (MQ=255)
                                                                                                      tttATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTg                                              >  1:1101365/1‑148 (MQ=255)
                                                                                                       ttATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGa                                             >  1:1917009/1‑148 (MQ=255)
                                                                                                         aTTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTGTTGACa                                           >  2:783523/1‑148 (MQ=255)
                                                                                                           tGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAg                                         >  1:1651128/1‑148 (MQ=255)
                                                                                                            gTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGc                                        <  1:78628/148‑1 (MQ=255)
                                                                                                            gTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGc                                        <  2:2917569/148‑1 (MQ=255)
                                                                                                                 cGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGcc                                   >  2:2393830/1‑148 (MQ=255)
                                                                                                                    ggATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGATTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACCTCACTTTTGACAATCTGGCCCCt                                >  1:2410688/1‑148 (MQ=255)
                                                                                                                    ggATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCt                                >  2:2830680/1‑148 (MQ=255)
                                                                                                                    ggATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCt                                >  1:489804/1‑148 (MQ=255)
                                                                                                                     gATAAACCTGGCATACCCCGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTg                               >  2:3157526/1‑148 (MQ=255)
                                                                                                                         aaCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGAGGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGt                           <  2:1795069/148‑1 (MQ=255)
                                                                                                                          aCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCTCGGTTGAGACTGAatccga                                                                  >  1:1469465/1‑108 (MQ=255)
                                                                                                                          aCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccga                                                                  <  2:3371649/108‑1 (MQ=255)
                                                                                                                          aCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccga                                                                  >  1:3371649/1‑108 (MQ=255)
                                                                                                                          aCCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAatccga                                                                  <  2:1469479/108‑1 (MQ=255)
                                                                                                                             tGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAg                       >  1:703270/1‑148 (MQ=255)
                                                                                                                                 aTACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACg                   <  2:1322844/148‑1 (MQ=255)
                                                                                                                                   aCCCCGCCCTGCCCTCGATTCAGGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCt                 <  2:2528501/148‑1 (MQ=255)
                                                                                                                                    ccccGCCCTGCCCTAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTg                >  2:2994367/1‑148 (MQ=255)
                                                                                                                                      ccGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTc              >  2:7928/1‑148 (MQ=255)
                                                                                                                                       cGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCa             <  2:151608/148‑1 (MQ=255)
                                                                                                                                            tgccctAGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGtt        >  2:3305498/1‑148 (MQ=255)
                                                                                                                                              ccctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTaa      >  2:3547856/1‑148 (MQ=255)
                                                                                                                                              ccctCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTCAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTaa      <  1:2251695/148‑1 (MQ=255)
                                                                                                                                                  cGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTc  >  2:3100211/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGATTCATGGGATGCTGAATATATGCAGGAAAACCGCTGGGATGAGGAGTCTGCTTCTTTTGTCCCGGTTGAGACTGAATCCGATCCGATGAACGTCACTTTTGACAAGCTGGCCCCTGAAGTACAGAACGCTGTCATGGTTAAGTTC  >  CP009273/1644404‑1644697

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: