Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,698,875 G→T 100% 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,8750GT99.5% 784.0 / ‑4.3 239F92L (TTC→TTAdppBdipeptide/heme transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (119/118);  total (120/119)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 239 total aligned reads displayed.

AACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAG  >  CP009273/3698730‑3699019
                                                                                                                                                 |                                                                                                                                                
aaCCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTaa                                                                                                                                                <  2:1223107/148‑1 (MQ=255)
    tGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAActct                                                                                                                                            >  2:486732/1‑148 (MQ=255)
       cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                         <  1:2523003/148‑1 (MQ=255)
       cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                         <  1:2627711/148‑1 (MQ=255)
       cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                         >  1:850638/1‑148 (MQ=255)
        aGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTccc                                                                                                                                        >  2:1203555/1‑148 (MQ=255)
        aGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTccc                                                                                                                                        >  2:3009001/1‑148 (MQ=255)
        aGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTccc                                                                                                                                        >  1:171755/1‑148 (MQ=255)
              aGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCg                                                                                                                                  >  2:3036942/1‑148 (MQ=255)
              aGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCg                                                                                                                                  >  2:877555/1‑148 (MQ=255)
              aGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGACCTCTTCCCAAACCg                                                                                                                                  >  1:2076805/1‑148 (MQ=255)
               ggCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCgg                                                                                                                                 >  1:528261/1‑148 (MQ=255)
               ggCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCgg                                                                                                                                 >  1:213576/1‑148 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                                >  2:1013393/1‑148 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                                <  1:1612490/148‑1 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                                <  1:3257433/148‑1 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                            <  2:2528835/148‑1 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                            >  1:1927041/1‑148 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                            <  2:769512/148‑1 (MQ=255)
                     ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCg                                                                                                                           >  2:1362106/1‑148 (MQ=255)
                     ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCg                                                                                                                           >  1:1052789/1‑148 (MQ=255)
                      cGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                          <  2:1269283/148‑1 (MQ=255)
                       gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGc                                                                                                                         >  2:1536704/1‑148 (MQ=255)
                       gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGc                                                                                                                         <  1:1930323/148‑1 (MQ=255)
                            gtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTc                                                                                                                    <  1:250152/148‑1 (MQ=255)
                             tgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCa                                                                                                                   >  1:3436587/1‑148 (MQ=255)
                             tgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCa                                                                                                                   >  1:1205856/1‑148 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  <  1:304797/148‑1 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  >  1:3201970/1‑148 (MQ=255)
                                tCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTg                                                                                                                <  1:1235434/148‑1 (MQ=255)
                                 cGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGa                                                                                                               >  1:921279/1‑148 (MQ=255)
                                   aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaaa                                                                                                             >  2:1308510/1‑148 (MQ=255)
                                   aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaaa                                                                                                             <  1:2052400/148‑1 (MQ=255)
                                    aaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAAt                                                                                                            >  2:1333596/1‑148 (MQ=255)
                                     aaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATg                                                                                                           <  1:1338419/148‑1 (MQ=255)
                                      aTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGc                                                                                                          <  2:1794078/148‑1 (MQ=255)
                                       tGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTGTTCATTGAAATGcc                                                                                                         >  2:682757/1‑148 (MQ=255)
                                         gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTa                                                                                                       >  2:1181741/1‑148 (MQ=255)
                                         gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTa                                                                                                       >  2:3087635/1‑148 (MQ=255)
                                          aaCCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAg                                                                                                      <  1:1006253/148‑1 (MQ=255)
                                          aaCCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAg                                                                                                      <  1:2781095/148‑1 (MQ=255)
                                          aaCCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAg                                                                                                      <  1:1470594/148‑1 (MQ=255)
                                            ccGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGAt                                                                                                    >  1:857644/1‑148 (MQ=255)
                                               cgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGc                                                                                                 >  1:3228396/1‑148 (MQ=255)
                                               cgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGc                                                                                                 <  1:2646002/148‑1 (MQ=255)
                                                gTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGcc                                                                                                <  2:1382082/148‑1 (MQ=255)
                                                    aaCCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGc                                                                                            <  2:1837121/148‑1 (MQ=255)
                                                       cGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATa                                                                                         <  1:1181741/148‑1 (MQ=255)
                                                                cACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                                <  1:2261186/148‑1 (MQ=255)
                                                                cACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                                >  1:1594947/1‑148 (MQ=255)
                                                                  cGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAt                                                                              <  1:601643/148‑1 (MQ=255)
                                                                     cGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTa                                                                           <  2:2188477/148‑1 (MQ=255)
                                                                     cGACCGGAATACCGACTGACGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGTGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTa                                                                           >  2:2844283/1‑148 (MQ=255)
                                                                       aCCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                         >  2:2478336/1‑148 (MQ=255)
                                                                       aCCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                         <  2:2970423/148‑1 (MQ=255)
                                                                       aCCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                         <  1:726803/148‑1 (MQ=255)
                                                                         cGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATgg                                                                       <  2:259926/148‑1 (MQ=255)
                                                                          ggAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGga                                                                      <  1:3105749/148‑1 (MQ=255)
                                                                          ggAATACCGACTGCCGTAGCAAAAATCATCGATCAGACGCCAAGTTACAGAGTGCCCTGGAAGCGCGGCACTAACTCTTCCCAAACAGCTATGCGTCTTTTCATAGAAATTCCAATAACGCCATGCATAACGCCCCAAAAGAAATCga                                                                      >  1:2797729/1‑148 (MQ=255)
                                                                           gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                     <  2:3479804/148‑1 (MQ=255)
                                                                           gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                     <  1:3391932/148‑1 (MQ=255)
                                                                              tACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                                  <  1:2832430/148‑1 (MQ=255)
                                                                                ccGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATAct                                                                >  2:145730/1‑148 (MQ=255)
                                                                                ccGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATAct                                                                >  2:2163102/1‑148 (MQ=255)
                                                                                  gACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgc                                                              <  2:3436587/148‑1 (MQ=255)
                                                                                     tGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACtgctgc                                                           <  2:2729925/148‑1 (MQ=255)
                                                                                       ccGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCca                                                         <  2:3282459/148‑1 (MQ=255)
                                                                                          tAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACAt                                                      <  2:677121/148‑1 (MQ=255)
                                                                                            gCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCg                                                    <  1:1264418/148‑1 (MQ=255)
                                                                                              aaaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGt                                                  <  2:1499665/148‑1 (MQ=255)
                                                                                              aaaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGATTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGt                                                  >  2:326990/1‑148 (MQ=255)
                                                                                               aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                 >  1:1953845/1‑148 (MQ=255)
                                                                                               aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                 >  1:179906/1‑148 (MQ=255)
                                                                                               aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                 >  1:1498282/1‑148 (MQ=255)
                                                                                                aaaTCATCGCGCAGTAGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGAGGATTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCACATCTAATGGAGATACTGCTGCCACATCGGttt                                                <  1:2579283/148‑1 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAATCCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAGGTAATGGAGATACTGCGGCCACATCGGtgt                                                >  2:2910328/1‑146 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                                >  1:1239749/1‑148 (MQ=255)
                                                                                                  atcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCGTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTAt                                              <  2:840396/148‑1 (MQ=255)
                                                                                                   tcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATc                                             <  2:160872/148‑1 (MQ=255)
                                                                                                    catcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCt                                            >  2:1932299/1‑148 (MQ=255)
                                                                                                     atcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTa                                           <  2:2347865/148‑1 (MQ=255)
                                                                                                       cgcgCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAg                                         >  1:2288870/1‑148 (MQ=255)
                                                                                                            aGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAg                                    <  1:1532115/148‑1 (MQ=255)
                                                                                                             gACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGt                                   <  2:280664/148‑1 (MQ=255)
                                                                                                              aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                                  <  1:1581271/148‑1 (MQ=255)
                                                                                                               cGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCCTTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTc                                 >  2:2000944/1‑148 (MQ=255)
                                                                                                                  cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAACTGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                              <  2:171755/148‑1 (MQ=255)
                                                                                                                  cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                              >  2:356251/1‑148 (MQ=255)
                                                                                                                  cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                              <  1:1376214/148‑1 (MQ=255)
                                                                                                                           gCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTg                     <  2:2710134/148‑1 (MQ=255)
                                                                                                                            cGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTgc                    >  2:2547143/1‑148 (MQ=255)
                                                                                                                             gTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTgcg                   <  2:1147766/148‑1 (MQ=255)
                                                                                                                                 ccTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTg               <  1:102795/148‑1 (MQ=255)
                                                                                                                                     gAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGc           <  1:990692/148‑1 (MQ=255)
                                                                                                                                      aaGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGct          <  2:1239749/148‑1 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct        <  2:3344004/148‑1 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct        >  2:511063/1‑148 (MQ=255)
                                                                                                                                           cgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTggg     >  1:2837106/1‑148 (MQ=255)
                                                                                                                                           cgGCACTAACTCGTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTggg     <  1:3036942/148‑1 (MQ=255)
                                                                                                                                             gCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGga   >  1:2712619/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAG  >  CP009273/3698730‑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: