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,8750GT95.9% 656.6 / 17.3 219F92L (TTC→TTAdppBdipeptide/heme transporter
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base T (98/111);  total (106/113)
Fisher's exact test for biased strand distribution p-value = 9.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 219 total aligned reads displayed.

AACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGG  >  CP009273/3698730‑3699015
                                                                                                                                                 |                                                                                                                                            
  ccTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAAct                                                                                                                                          >  1:1165650/1‑148 (MQ=255)
  ccTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAAct                                                                                                                                          >  1:359571/1‑148 (MQ=255)
  ccTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAAct                                                                                                                                          <  1:2696707/148‑1 (MQ=255)
   cTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAActc                                                                                                                                         <  2:609020/148‑1 (MQ=255)
   cTGTCAGCGCCAGGCCAACCGCGGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAActc                                                                                                                                         <  2:846347/148‑1 (MQ=255)
     gTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCtt                                                                                                                                       <  2:463796/148‑1 (MQ=255)
       cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                     <  1:1280034/148‑1 (MQ=255)
       cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                     <  1:167919/148‑1 (MQ=255)
         gcgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCa                                                                                                                                   >  1:1008536/1‑148 (MQ=255)
         gcgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCa                                                                                                                                   >  1:2623826/1‑148 (MQ=255)
          cgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCaa                                                                                                                                  >  2:3315539/1‑148 (MQ=255)
              aGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCg                                                                                                                              >  1:1094721/1‑148 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                            <  1:243145/148‑1 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                            >  2:3128382/1‑148 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                            >  2:581465/1‑148 (MQ=255)
                   aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATg                                                                                                                         <  2:578231/148‑1 (MQ=255)
                   aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATg                                                                                                                         >  2:237880/1‑148 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                        >  1:2881343/1‑148 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                        >  1:2999336/1‑148 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                        >  2:1315054/1‑148 (MQ=255)
                     ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCg                                                                                                                       <  1:1794949/148‑1 (MQ=255)
                       gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAACCCGGGATGCGGc                                                                                                                     <  2:1703735/148‑1 (MQ=255)
                        cTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                    >  1:3153779/1‑148 (MQ=255)
                        cTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                    <  1:3053998/148‑1 (MQ=255)
                          gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCttt                                                                                                                  <  1:3371365/148‑1 (MQ=255)
                           tgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                 <  2:10604/148‑1 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACAGGGATGCGGCTTTTCAt                                                                                                              >  1:1921326/1‑148 (MQ=255)
                               aTCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGACGCGCGGCACTAACTCTTCCCACACCGGGATGCGGCTTTTCAtt                                                                                                             >  2:2975131/1‑148 (MQ=255)
                               aTCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAtt                                                                                                             >  2:2542818/1‑148 (MQ=255)
                                tCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTg                                                                                                            >  2:2390594/1‑148 (MQ=255)
                                tCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTg                                                                                                            >  2:889165/1‑148 (MQ=255)
                                  gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                          <  2:3140698/148‑1 (MQ=255)
                                  gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                          <  1:2675364/148‑1 (MQ=255)
                                     aaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATg                                                                                                       <  1:327364/148‑1 (MQ=255)
                                      aTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGc                                                                                                      <  1:1315054/148‑1 (MQ=255)
                                        ggAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAAGACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAACCCGGGATGCGGCTTTTCATTGAAATGCCt                                                                                                    <  2:2767271/148‑1 (MQ=255)
                                         gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTa                                                                                                   <  1:3479349/148‑1 (MQ=255)
                                            ccGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGAt                                                                                                >  1:594741/1‑148 (MQ=255)
                                               cgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGc                                                                                             <  1:752036/148‑1 (MQ=255)
                                                gTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGTCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGcc                                                                                            <  1:664729/148‑1 (MQ=255)
                                                gTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGcc                                                                                            <  1:969854/148‑1 (MQ=255)
                                                     aCCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCa                                                                                       >  1:3247812/1‑148 (MQ=255)
                                                      ccGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCAt                                                                                      >  1:2280281/1‑148 (MQ=255)
                                                       cGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATa                                                                                     >  2:3455121/1‑148 (MQ=255)
                                                          agccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACg                                                                                  <  1:1437993/148‑1 (MQ=255)
                                                          agccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACg                                                                                  <  1:1439795/148‑1 (MQ=255)
                                                          agccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACg                                                                                  <  1:1598453/148‑1 (MQ=255)
                                                          agccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACg                                                                                  <  2:1371323/148‑1 (MQ=255)
                                                               gcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCa                                                                             >  1:2793049/1‑148 (MQ=255)
                                                               gcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCa                                                                             >  2:3430495/1‑148 (MQ=255)
                                                                cACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                            <  1:378376/148‑1 (MQ=255)
                                                                cACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                            >  1:716480/1‑148 (MQ=255)
                                                                   gCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATg                                                                         <  1:2978600/148‑1 (MQ=255)
                                                                    ccGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGt                                                                        >  2:2195710/1‑148 (MQ=255)
                                                                       aCCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                     <  2:2670707/148‑1 (MQ=255)
                                                                        ccGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATg                                                                    <  2:2061984/148‑1 (MQ=255)
                                                                        ccGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATg                                                                    <  2:967323/148‑1 (MQ=255)
                                                                           gAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                 >  1:2455784/1‑148 (MQ=255)
                                                                            aaTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgaga                                                                <  1:2131267/148‑1 (MQ=255)
                                                                            aaTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgaga                                                                <  1:2350305/148‑1 (MQ=255)
                                                                             aTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGAt                                                               <  2:607155/148‑1 (MQ=255)
                                                                              tACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGGGATa                                                              <  2:557720/148‑1 (MQ=255)
                                                                              tACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                              >  2:2890087/1‑148 (MQ=255)
                                                                                 cGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActg                                                           <  2:3530913/148‑1 (MQ=255)
                                                                                  gACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgc                                                          <  2:258091/148‑1 (MQ=255)
                                                                                         gTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCcaca                                                   <  1:1450258/148‑1 (MQ=255)
                                                                                             cAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCgg                                               >  2:3100385/1‑148 (MQ=255)
                                                                                             cAAAAATCATCGCGCAGACGCCAAGTTCAAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCgg                                               >  1:2501931/1‑148 (MQ=255)
                                                                                              aaaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGt                                              <  1:254008/148‑1 (MQ=255)
                                                                                               aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                             >  2:660895/1‑148 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                            >  1:1856826/1‑148 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                            <  2:911484/148‑1 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                            <  2:911625/148‑1 (MQ=255)
                                                                                                 aaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTa                                           <  1:1713125/148‑1 (MQ=255)
                                                                                                 aaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTa                                           >  1:531413/1‑148 (MQ=255)
                                                                                                  atcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTAt                                          >  2:547787/1‑148 (MQ=255)
                                                                                                  atcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTAt                                          <  2:2089639/148‑1 (MQ=255)
                                                                                                       cgcgCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAg                                     <  2:2300659/148‑1 (MQ=255)
                                                                                                         cgcAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGcc                                   >  2:2403438/1‑148 (MQ=255)
                                                                                                         cgcAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGcc                                   >  2:1207547/1‑148 (MQ=255)
                                                                                                            aGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTTCTGCCACATCGGTTTATCTAAGCCGAg                                >  1:1068844/1‑148 (MQ=255)
                                                                                                            aGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAg                                >  1:3488999/1‑148 (MQ=255)
                                                                                                              aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                              >  2:1770849/1‑148 (MQ=255)
                                                                                                              aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                              >  2:431026/1‑148 (MQ=255)
                                                                                                               cGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTc                             <  2:1685280/148‑1 (MQ=255)
                                                                                                                  cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTCATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                          >  1:2722694/1‑148 (MQ=255)
                                                                                                                  cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                          <  2:343537/148‑1 (MQ=255)
                                                                                                                   aaGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagcc                         >  2:500484/1‑148 (MQ=255)
                                                                                                                   aaGTTCCAGCGTGGCCGGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagcc                         <  2:3137825/148‑1 (MQ=255)
                                                                                                                        ccAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCcagcag                    <  2:1008536/148‑1 (MQ=255)
                                                                                                                          aGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCt                  >  1:2972834/1‑148 (MQ=255)
                                                                                                                          aGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCt                  >  1:1612191/1‑148 (MQ=255)
                                                                                                                               ggCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCAGAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGcgcg             <  2:1723443/148‑1 (MQ=255)
                                                                                                                                     gAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGc       >  2:175200/1‑148 (MQ=255)
                                                                                                                                      aaGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGct      >  1:13880/1‑148 (MQ=255)
                                                                                                                                       aGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctc     <  2:346932/148‑1 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct    >  2:3047474/1‑148 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct    <  1:1721698/148‑1 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct    <  1:1259018/148‑1 (MQ=255)
                                                                                                                                          gcgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTgg  <  1:1871757/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
AACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGG  >  CP009273/3698730‑3699015

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: