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.1% 700.6 / 22.3 235F92L (TTC→TTAdppBdipeptide/heme transporter
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (114/109);  total (118/117)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 235 total aligned reads displayed.

AACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAG  >  CP009273/3698730‑3699019
                                                                                                                                                 |                                                                                                                                                
aaCCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTaa                                                                                                                                                <  2:1728332/148‑1 (MQ=255)
aaCCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGaa                                                                                                                                                >  2:60401/1‑148 (MQ=255)
  ccTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAAct                                                                                                                                              >  1:1747650/1‑148 (MQ=255)
     gTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCtt                                                                                                                                           <  2:2543034/148‑1 (MQ=255)
       cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                         <  2:116462/148‑1 (MQ=255)
         gcgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCa                                                                                                                                       >  2:1295427/1‑148 (MQ=255)
           gcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCaaa                                                                                                                                     <  2:3455107/148‑1 (MQ=255)
            ccAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAAc                                                                                                                                    >  2:558977/1‑148 (MQ=255)
             cAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAAcc                                                                                                                                   <  2:783163/148‑1 (MQ=255)
             cAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAAcc                                                                                                                                   >  1:386097/1‑148 (MQ=255)
               ggCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCgg                                                                                                                                 >  2:1694313/1‑148 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                                >  2:1723772/1‑148 (MQ=255)
                gCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCggg                                                                                                                                <  2:1471338/148‑1 (MQ=255)
                   aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATg                                                                                                                             >  1:2447223/1‑148 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                            <  1:2731944/148‑1 (MQ=255)
                    aCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                            >  1:3374351/1‑148 (MQ=255)
                     ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCg                                                                                                                           >  2:1201844/1‑148 (MQ=255)
                      cGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                          <  1:743588/148‑1 (MQ=255)
                      cGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                          <  1:1394913/148‑1 (MQ=255)
                       gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGc                                                                                                                         >  1:669493/1‑148 (MQ=255)
                        cTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                        <  2:931466/148‑1 (MQ=255)
                        cTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                        <  1:2295500/148‑1 (MQ=255)
                         tgtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtt                                                                                                                       <  2:137227/148‑1 (MQ=255)
                          gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCttt                                                                                                                      <  1:3257969/148‑1 (MQ=255)
                          gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCttt                                                                                                                      >  2:778665/1‑148 (MQ=255)
                           tgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                     <  2:3227933/148‑1 (MQ=255)
                            gtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTc                                                                                                                    <  1:2035468/148‑1 (MQ=255)
                             tgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCa                                                                                                                   >  2:2682334/1‑148 (MQ=255)
                             tgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCa                                                                                                                   >  2:3049635/1‑148 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  <  1:1893439/148‑1 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  <  1:1528476/148‑1 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGACACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  >  2:1601011/1‑148 (MQ=255)
                              gATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGCAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAt                                                                                                                  >  1:853267/1‑148 (MQ=255)
                               aTCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAtt                                                                                                                 >  2:3525521/1‑148 (MQ=255)
                                  gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                              <  1:3061527/148‑1 (MQ=255)
                                  gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                              <  2:2302416/148‑1 (MQ=255)
                                  gAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAAGACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaa                                                                                                              <  2:2035587/148‑1 (MQ=255)
                                   aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGaaa                                                                                                             >  2:320042/1‑148 (MQ=255)
                                    aaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAAt                                                                                                            >  1:3173354/1‑148 (MQ=255)
                                      aTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGc                                                                                                          <  1:1315529/148‑1 (MQ=255)
                                       tGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGcc                                                                                                         <  1:2390200/148‑1 (MQ=255)
                                        ggAACCACGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCt                                                                                                        >  1:907748/1‑148 (MQ=255)
                                         gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTa                                                                                                       >  1:1608874/1‑148 (MQ=255)
                                         gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTa                                                                                                       >  2:1123840/1‑148 (MQ=255)
                                            ccGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGAt                                                                                                    <  1:1619021/148‑1 (MQ=255)
                                            ccGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGAt                                                                                                    <  1:2702614/148‑1 (MQ=255)
                                               cgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTTGATGCGCGGCACTAACTCTTCCCTTACCTGGATGCTGCTTTTCATTGAAATGCCTAGATCGc                                                                                                 <  2:1312504/148‑1 (MQ=255)
                                               cgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGc                                                                                                 <  2:1190156/148‑1 (MQ=255)
                                               cgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGc                                                                                                 <  2:2968582/148‑1 (MQ=255)
                                                      ccGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCAt                                                                                          >  2:153662/1‑148 (MQ=255)
                                                      ccGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCAt                                                                                          >  2:2956695/1‑148 (MQ=255)
                                                        gCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATaa                                                                                        <  1:3543784/148‑1 (MQ=255)
                                                          agccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACg                                                                                      >  1:2271999/1‑148 (MQ=255)
                                                           gccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGc                                                                                     >  1:1063561/1‑148 (MQ=255)
                                                           gccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCAGTGAAATGCCTAGATCGCCATGCATAACGc                                                                                     >  1:2168478/1‑148 (MQ=255)
                                                            ccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGcc                                                                                    <  1:1723772/148‑1 (MQ=255)
                                                                cACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                                <  1:2077581/148‑1 (MQ=255)
                                                                  cGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAt                                                                              <  2:456700/148‑1 (MQ=255)
                                                                  cGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAt                                                                              >  1:802841/1‑148 (MQ=255)
                                                                     cGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTa                                                                           <  2:1888935/148‑1 (MQ=255)
                                                                       aCCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                         >  2:3144179/1‑148 (MQ=255)
                                                                       aCCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                         <  1:2185317/148‑1 (MQ=255)
                                                                        ccGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATg                                                                        >  2:2106965/1‑148 (MQ=255)
                                                                        ccGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATg                                                                        >  2:2106917/1‑148 (MQ=255)
                                                                          ggAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGga                                                                      <  2:1747650/148‑1 (MQ=255)
                                                                            aaTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgaga                                                                    <  1:2258079/148‑1 (MQ=255)
                                                                             aTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGAt                                                                   <  1:1049614/148‑1 (MQ=255)
                                                                              tACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                                  <  2:284337/148‑1 (MQ=255)
                                                                               aCCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAAGGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATAc                                                                 <  2:2011255/148‑1 (MQ=255)
                                                                                 cGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGACACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActg                                                               <  1:1600981/148‑1 (MQ=255)
                                                                                   aCTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgct                                                             >  1:1883081/1‑148 (MQ=255)
                                                                                      gCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGcc                                                          >  2:1585727/1‑148 (MQ=255)
                                                                                            gCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCg                                                    >  1:41113/1‑148 (MQ=255)
                                                                                               aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                 >  1:417713/1‑148 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                                >  1:3258955/1‑148 (MQ=255)
                                                                                                aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                                <  1:384480/148‑1 (MQ=255)
                                                                                                 aaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTa                                               >  2:2977746/1‑148 (MQ=255)
                                                                                                   tcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATc                                             <  1:1123840/148‑1 (MQ=255)
                                                                                                     atcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCtc                                           >  1:2595213/1‑147 (MQ=255)
                                                                                                           cAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGa                                     >  2:126913/1‑148 (MQ=255)
                                                                                                            aGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAg                                    <  2:1623236/148‑1 (MQ=255)
                                                                                                            aGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAg                                    >  1:2368572/1‑148 (MQ=255)
                                                                                                             gACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGt                                   <  2:907748/148‑1 (MQ=255)
                                                                                                              aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                                  >  2:383366/1‑148 (MQ=255)
                                                                                                              aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                                  <  1:3405718/148‑1 (MQ=255)
                                                                                                               cGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTc                                 >  2:2462262/1‑148 (MQ=255)
                                                                                                                  cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagc                              <  1:1705262/148‑1 (MQ=255)
                                                                                                                   aaGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagcc                             >  1:428265/1‑148 (MQ=255)
                                                                                                                   aaGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagcc                             >  2:1982354/1‑148 (MQ=255)
                                                                                                                     gTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagccag                           <  2:231435/148‑1 (MQ=255)
                                                                                                                     gTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagccag                           >  1:2809923/1‑148 (MQ=255)
                                                                                                                     gTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagccag                           >  1:2809088/1‑148 (MQ=255)
                                                                                                                           gCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTg                     >  1:2536244/1‑148 (MQ=255)
                                                                                                                            cGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTgc                    >  1:1681986/1‑148 (MQ=255)
                                                                                                                                      aaGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGct          <  2:261594/148‑1 (MQ=255)
                                                                                                                                      aaGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGct          >  1:29962/1‑148 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct        <  1:1271619/148‑1 (MQ=255)
                                                                                                                                        gcgcGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctct        <  1:3049635/148‑1 (MQ=255)
                                                                                                                                          gcgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTgg      <  2:2193473/148‑1 (MQ=255)
                                                                                                                                             gCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGga   <  1:558977/148‑1 (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: