Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,858 G→T 39.7% V71V (GTC→GTA pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,8580GT39.7% 235.7 / 198.5 194V71V (GTC→GTApepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (58/59);  new base T (27/50);  total (85/109)
Fisher's exact test for biased strand distribution p-value = 5.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 194 total aligned reads displayed.

CTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAG  >  CP009273/4219712‑4220004
                                                                                                                                                  |                                                                                                                                                  
   cAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCGACaa                                                                                                                                                <  1:3360130/148‑1 (MQ=255)
      ccAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGc                                                                                                                                             >  1:706649/1‑148 (MQ=255)
      ccAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGc                                                                                                                                             >  1:1556726/1‑148 (MQ=255)
      ccAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGc                                                                                                                                             >  2:205133/1‑148 (MQ=255)
       cAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCGACAACGCt                                                                                                                                            <  1:458369/148‑1 (MQ=255)
       cAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCt                                                                                                                                            >  2:1590679/1‑148 (MQ=255)
       cAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCt                                                                                                                                            >  2:2150161/1‑148 (MQ=255)
         atatACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACACTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTAt                                                                                                                                          >  2:3012641/1‑148 (MQ=255)
           ataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGa                                                                                                                                        <  1:3082493/148‑1 (MQ=255)
            taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGaa                                                                                                                                       <  2:1981981/148‑1 (MQ=255)
            taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCAGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGaa                                                                                                                                       >  2:2780919/1‑148 (MQ=255)
              cAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATa                                                                                                                                     <  2:3394418/148‑1 (MQ=255)
                      cACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGAc                                                                                                                             <  2:981869/148‑1 (MQ=255)
                        cGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAg                                                                                                                           >  2:3051470/1‑148 (MQ=255)
                         gTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGa                                                                                                                          <  1:1589926/148‑1 (MQ=255)
                          tttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGaa                                                                                                                         <  2:3325253/148‑1 (MQ=255)
                          tttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGaa                                                                                                                         >  2:1315799/1‑148 (MQ=255)
                            tAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAac                                                                                                                       <  2:1932765/148‑1 (MQ=255)
                            tAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACAGGATCTACAACGCTATGAATACCGGTGACAGAAac                                                                                                                       <  2:131073/148‑1 (MQ=255)
                             aaCCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAaca                                                                                                                      <  1:455486/148‑1 (MQ=255)
                               ccACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAcc                                                                                                                    >  1:2211248/1‑148 (MQ=255)
                                cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATTCCGGTGACAGAAACAccc                                                                                                                   <  1:617787/148‑1 (MQ=255)
                                cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACAccc                                                                                                                   >  1:2541943/1‑148 (MQ=255)
                                  cGTCAGTAATTGGTGCCAGAAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAg                                                                                                                 <  1:2621184/148‑1 (MQ=255)
                                      aGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGa                                                                                                             <  2:3108678/148‑1 (MQ=255)
                                        tAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAATGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGc                                                                                                           >  1:2064966/1‑148 (MQ=255)
                                         aaTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCg                                                                                                          <  2:3086578/148‑1 (MQ=255)
                                          aTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCga                                                                                                         >  2:970650/1‑148 (MQ=255)
                                             ggTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGaa                                                                                                      <  1:2087970/148‑1 (MQ=255)
                                                gCCAGCAGCCCGCACTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCg                                                                                                   <  1:3166371/148‑1 (MQ=255)
                                                     cagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                                                                                              <  1:60999/148‑1 (MQ=255)
                                                        cccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                                                                                           >  2:1423258/1‑148 (MQ=255)
                                                        cccGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                                                                                           >  2:3544685/1‑148 (MQ=255)
                                                         ccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGttt                                                                                          <  2:336914/148‑1 (MQ=255)
                                                               tCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtg                                                                                    <  1:1923164/148‑1 (MQ=255)
                                                               tCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtg                                                                                    <  1:1840104/148‑1 (MQ=255)
                                                               tCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCtg                                                                                    <  1:35632/148‑1 (MQ=255)
                                                                  cgGCACTGTTTCCGCAACTGGAAAGTATTCCCGCCGCCGACACTCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTgtgt                                                                                 <  2:1730126/148‑1 (MQ=255)
                                                                    gCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTaa                                                                               >  1:1360920/1‑148 (MQ=255)
                                                                    gCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTaa                                                                               <  2:350188/148‑1 (MQ=255)
                                                                    gCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTaa                                                                               <  2:350195/148‑1 (MQ=255)
                                                                    gCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTaa                                                                               >  1:856965/1‑148 (MQ=255)
                                                                     cACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAat                                                                              >  1:803248/1‑148 (MQ=255)
                                                                      aCTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatc                                                                             >  1:3264367/1‑148 (MQ=255)
                                                                      aCTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatc                                                                             <  2:706649/148‑1 (MQ=255)
                                                                        tGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcat                                                                           >  1:497028/1‑148 (MQ=255)
                                                                         gTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcatc                                                                          <  1:2365506/148‑1 (MQ=255)
                                                                         gTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACTGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcatc                                                                          <  2:2381500/148‑1 (MQ=255)
                                                                         gTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcatc                                                                          >  2:150336/1‑148 (MQ=255)
                                                                            tCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCa                                                                       <  1:2155430/148‑1 (MQ=255)
                                                                              tgCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAgg                                                                     >  1:2098438/2‑148 (MQ=255)
                                                                                  aCTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTg                                                                 <  1:1116983/148‑1 (MQ=255)
                                                                                    tGGAAAGTATTCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCg                                                               <  1:724742/148‑1 (MQ=255)
                                                                                     ggAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGt                                                              <  2:658347/148‑1 (MQ=255)
                                                                                       aaaGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTa                                                            <  1:1219737/148‑1 (MQ=255)
                                                                                        aaGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTAc                                                           >  1:1637162/1‑148 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                     <  1:2447069/148‑1 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                     >  2:2736768/1‑148 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                     >  1:3079239/1‑148 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                     >  1:1722652/1‑148 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                     <  1:427845/148‑1 (MQ=255)
                                                                                              tCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGc                                                     <  2:1713299/148‑1 (MQ=255)
                                                                                                  gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaaa                                                 <  2:862474/148‑1 (MQ=255)
                                                                                                  gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaaa                                                 >  2:3443388/1‑148 (MQ=255)
                                                                                                  gccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGaaa                                                 <  1:2524248/148‑1 (MQ=255)
                                                                                                   ccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAg                                                <  2:3476405/148‑1 (MQ=255)
                                                                                                   ccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAg                                                <  2:1360865/148‑1 (MQ=255)
                                                                                                   ccgccgACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAg                                                >  1:1700595/1‑148 (MQ=255)
                                                                                                         aCAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaa                                          >  2:2863631/1‑148 (MQ=255)
                                                                                                         aCAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaa                                          >  1:206864/1‑148 (MQ=255)
                                                                                                          cAATCACGATCTCCGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATaaa                                         <  1:826158/148‑1 (MQ=255)
                                                                                                               aCGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCg                                    <  2:497028/148‑1 (MQ=255)
                                                                                                               aCGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCg                                    <  1:1555666/148‑1 (MQ=255)
                                                                                                                cGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGc                                   >  1:1073752/1‑148 (MQ=255)
                                                                                                                  aTCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTg                                 >  1:2730998/1‑148 (MQ=255)
                                                                                                                    ctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAg                               <  2:740320/148‑1 (MQ=255)
                                                                                                                    ctcAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAg                               <  2:3558470/148‑1 (MQ=255)
                                                                                                                          aTTTTTAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGa                         >  1:1644518/1‑148 (MQ=255)
                                                                                                                          aTTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGa                         >  1:2406478/1‑148 (MQ=255)
                                                                                                                            tttCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTGTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAcc                       >  1:1337257/1‑148 (MQ=255)
                                                                                                                            tttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAcc                       >  2:2664092/1‑148 (MQ=255)
                                                                                                                            tttCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGAcc                       >  2:1221481/1‑148 (MQ=255)
                                                                                                                              tCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCt                     >  1:377908/1‑148 (MQ=255)
                                                                                                                               cAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTg                    >  2:2310975/1‑148 (MQ=255)
                                                                                                                                aaTCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGc                   <  1:1926892/148‑1 (MQ=255)
                                                                                                                                  tCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCaa                 <  2:1432240/148‑1 (MQ=255)
                                                                                                                                     cgGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTg              <  1:1779754/148‑1 (MQ=255)
                                                                                                                                     cgGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTg              <  1:1590679/148‑1 (MQ=255)
                                                                                                                                       gCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGtt            <  2:2211248/148‑1 (MQ=255)
                                                                                                                                        cAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTc           <  1:578917/148‑1 (MQ=255)
                                                                                                                                           cGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGc        <  1:451623/148‑1 (MQ=255)
                                                                                                                                           cGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGc        <  1:2244852/148‑1 (MQ=255)
                                                                                                                                            gggATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCa       <  1:905169/148‑1 (MQ=255)
                                                                                                                                             ggATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCaa      <  2:2032013/148‑1 (MQ=255)
                                                                                                                                              gATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAt     >  2:316762/1‑148 (MQ=255)
                                                                                                                                              gATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAt     >  2:1888430/1‑148 (MQ=255)
                                                                                                                                               aTCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt    <  1:2294022/148‑1 (MQ=255)
                                                                                                                                               aTCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt    <  1:1315760/148‑1 (MQ=255)
                                                                                                                                               aTCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAAtt    <  1:92947/148‑1 (MQ=255)
                                                                                                                                                 cGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAg  <  2:2235576/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAGGTCTGCGTTACGCCAGCGAAAGGGATAAACACCGCTGAGCGGCGACCCTGCAACTGTTCAGCAATTAG  >  CP009273/4219712‑4220004

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: