Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,698,875 G→T 90.0% F92L (TTC→TTA)  dppB ← dipeptide/heme transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,698,8750GT90.0% 566.3 / 59.7 220F92L (TTC→TTA) dppBdipeptide/heme transporter
Reads supporting (aligned to +/- strand):  ref base G (10/12);  new base T (92/106);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 220 total aligned reads displayed.

CTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAG  >  CP009273/3698733‑3699019
                                                                                                                                              |                                                                                                                                                
 tGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTCGCCAAAATCATCGCGCAGACGCCCAGTTCCAGCGTGGCCTGGAAGCGCGGCCCTAActct                                                                                                                                            >  2:1911759/1‑148 (MQ=255)
   tCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTc                                                                                                                                          >  1:2149522/1‑148 (MQ=255)
      gcgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCa                                                                                                                                       >  2:2337141/1‑148 (MQ=255)
       cgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCaa                                                                                                                                      >  1:1539774/1‑148 (MQ=255)
          cagccaaaCCGATGAGTGATCGAAAATGCAACCGAGTTTAAACGCAGACAGCAAGCCGAACGTAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGATCCATCGAGGCCTGGAAGCGCGGCACTAACTCTTCCCAAAca                                                                                                                                   >  2:2635347/7‑147 (MQ=255)
          cAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAAcc                                                                                                                                   >  2:232707/1‑148 (MQ=255)
           aGGCCCACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCg                                                                                                                                  <  2:3094869/148‑1 (MQ=255)
           aGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCg                                                                                                                                  >  2:217237/1‑148 (MQ=255)
              ccAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGa                                                                                                                               <  1:2334456/148‑1 (MQ=255)
              ccAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGa                                                                                                                               <  1:3162383/148‑1 (MQ=255)
               cAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGAt                                                                                                                              <  1:1461167/148‑1 (MQ=255)
                aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATg                                                                                                                             <  2:240838/148‑1 (MQ=255)
                aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATg                                                                                                                             >  1:813020/1‑148 (MQ=255)
                aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATg                                                                                                                             <  2:2845147/148‑1 (MQ=255)
                  ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCg                                                                                                                           >  2:2553357/1‑148 (MQ=255)
                   cGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                          <  1:1901415/148‑1 (MQ=255)
                   cGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                          <  2:2431566/148‑1 (MQ=255)
                     cTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                        >  1:3078799/1‑148 (MQ=255)
                     cTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                        >  1:2785662/1‑148 (MQ=255)
                        tgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAACAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                     <  2:2801798/148‑1 (MQ=255)
                        tgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                     <  2:700104/148‑1 (MQ=255)
                         gtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTc                                                                                                                    <  2:861302/148‑1 (MQ=255)
                           gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  <  1:2699270/148‑1 (MQ=255)
                            aTCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAtt                                                                                                                 >  1:970800/1‑148 (MQ=255)
                            aTCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAtt                                                                                                                 >  1:1279267/1‑148 (MQ=255)
                              cGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGa                                                                                                               >  2:1935906/1‑148 (MQ=255)
                               gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                              <  1:867285/148‑1 (MQ=255)
                               gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                              <  1:3155859/148‑1 (MQ=255)
                                aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaaa                                                                                                             >  1:2509257/1‑148 (MQ=255)
                                aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaaa                                                                                                             <  1:862575/148‑1 (MQ=255)
                                aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaaa                                                                                                             >  2:2619457/1‑148 (MQ=255)
                                   aTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGc                                                                                                          <  2:1286454/148‑1 (MQ=255)
                                   aTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGc                                                                                                          <  1:2552659/148‑1 (MQ=255)
                                   aTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGc                                                                                                          <  1:2380936/148‑1 (MQ=255)
                                    tGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGcc                                                                                                         >  1:2850830/1‑148 (MQ=255)
                                    tGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGcc                                                                                                         >  2:857356/1‑148 (MQ=255)
                                    tGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGcc                                                                                                         <  1:2036838/148‑1 (MQ=255)
                                      gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTa                                                                                                       >  1:1166825/1‑148 (MQ=255)
                                       aaCCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAg                                                                                                      >  1:380589/1‑148 (MQ=255)
                                         ccGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGAt                                                                                                    <  1:1634237/148‑1 (MQ=255)
                                               ttAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCAt                                                                                              >  2:1553799/1‑148 (MQ=255)
                                                tAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCGTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATg                                                                                             <  2:2987063/148‑1 (MQ=255)
                                                 aaCCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGc                                                                                            >  1:1803143/1‑148 (MQ=255)
                                                 aaCCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGc                                                                                            <  2:1835655/148‑1 (MQ=255)
                                                     gCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATaa                                                                                        >  1:147027/1‑148 (MQ=255)
                                                     gCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATaa                                                                                        <  2:1539774/148‑1 (MQ=255)
                                                     gCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATaa                                                                                        <  2:3181016/148‑1 (MQ=255)
                                                       agccagcACGCCGACCGGAATACCGACTGCCGTACCAACAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACg                                                                                      <  2:1768201/148‑1 (MQ=255)
                                                        gccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGc                                                                                     >  2:1620111/1‑148 (MQ=255)
                                                         ccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGcc                                                                                    <  1:727603/148‑1 (MQ=255)
                                                          cagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGccc                                                                                   <  1:2399613/148‑1 (MQ=255)
                                                            gcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCa                                                                                 >  1:2576943/1‑148 (MQ=255)
                                                            gcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCa                                                                                 >  1:2577146/1‑148 (MQ=255)
                                                            gcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCa                                                                                 >  1:1115366/1‑148 (MQ=255)
                                                             cACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                                <  1:1137005/148‑1 (MQ=255)
                                                               cGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAt                                                                              <  2:2928261/148‑1 (MQ=255)
                                                                  cGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCTAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTa                                                                           <  2:2637728/148‑1 (MQ=255)
                                                                   gACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTaa                                                                          <  1:483712/148‑1 (MQ=255)
                                                                        gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                     >  1:2350015/1‑148 (MQ=255)
                                                                        gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                     <  1:607130/148‑1 (MQ=255)
                                                                        gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                     <  1:2268968/148‑1 (MQ=255)
                                                                        gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                     >  1:3190072/1‑148 (MQ=255)
                                                                          aTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGAt                                                                   >  1:1386292/1‑148 (MQ=255)
                                                                           tACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                                  <  2:104790/148‑1 (MQ=255)
                                                                            aCCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATAc                                                                 <  1:1104347/148‑1 (MQ=255)
                                                                              cGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActg                                                               >  1:1941296/1‑148 (MQ=255)
                                                                                  tGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACtgctgc                                                           <  2:2248882/148‑1 (MQ=255)
                                                                                    ccGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCca                                                         >  2:3105157/1‑148 (MQ=255)
                                                                                     cGTAGCAAAAATCCTCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCcac                                                        <  1:2535591/148‑1 (MQ=255)
                                                                                     cGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCcac                                                        >  2:599493/1‑148 (MQ=255)
                                                                                       tAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACAt                                                      <  1:318546/148‑1 (MQ=255)
                                                                                       tAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACAt                                                      <  2:1950483/148‑1 (MQ=255)
                                                                                       tAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACAt                                                      >  2:402553/1‑148 (MQ=255)
                                                                                       tAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACAt                                                      <  2:32022/148‑1 (MQ=255)
                                                                                           aaaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGt                                                  <  1:2300296/148‑1 (MQ=255)
                                                                                            aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                 <  1:1576576/148‑1 (MQ=255)
                                                                                             aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                                <  2:2248328/148‑1 (MQ=255)
                                                                                             aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                                >  1:499113/1‑148 (MQ=255)
                                                                                               atcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTAt                                              >  1:1034590/1‑148 (MQ=255)
                                                                                               atcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTAt                                              <  2:650846/148‑1 (MQ=255)
                                                                                                   tcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTaa                                          <  1:557128/148‑1 (MQ=255)
                                                                                                   tcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTaa                                          >  2:1207349/1‑148 (MQ=255)
                                                                                                   tcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTaa                                          >  2:3289701/1‑148 (MQ=255)
                                                                                                    cgcgCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAg                                         >  2:3295455/1‑148 (MQ=255)
                                                                                                    cgcgCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAg                                         >  1:2072897/1‑148 (MQ=255)
                                                                                                     gcgcAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGc                                        <  2:1447829/148‑1 (MQ=255)
                                                                                                           aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                                  <  2:2049493/148‑1 (MQ=255)
                                                                                                               cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                              <  2:2528559/148‑1 (MQ=255)
                                                                                                                  gTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagccag                           >  2:3252484/1‑148 (MQ=255)
                                                                                                                      cAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCagcagc                       <  2:2362271/148‑1 (MQ=255)
                                                                                                                              ccTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTg               <  1:2961714/148‑1 (MQ=255)
                                                                                                                              ccTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTg               <  2:1135916/148‑1 (MQ=255)
                                                                                                                                tGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGAc             >  1:683416/1‑148 (MQ=255)
                                                                                                                                 ggAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACg            <  1:1935906/148‑1 (MQ=255)
                                                                                                                                      cgcgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTg       <  2:2509257/148‑1 (MQ=255)
                                                                                                                                      cgcgGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTg       <  2:3191948/148‑1 (MQ=255)
                                                                                                                                        cgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTggg     <  2:2785662/148‑1 (MQ=255)
                                                                                                                                         ggCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTgggg    <  2:486153/148‑1 (MQ=255)
                                                                                                                                         ggCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTgggg    <  2:2149522/148‑1 (MQ=255)
                                                                                                                                           cACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGgag  <  1:1308622/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAG  >  CP009273/3698733‑3699019

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