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,8750GT98.6% 915.7 / 2.7 280F92L (TTC→TTAdppBdipeptide/heme transporter
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (135/140);  total (139/141)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
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 280 total aligned reads displayed.

TAACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGAT  >  CP009273/3698729‑3699021
                                                                                                                                                  |                                                                                                                                                  
tAACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTaa                                                                                                                                                  >  1:2773822/1‑149 (MQ=255)
tAACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTaa                                                                                                                                                  >  1:3417677/1‑149 (MQ=255)
  aCCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAAct                                                                                                                                                >  1:2175772/1‑149 (MQ=255)
   ccTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAActc                                                                                                                                               <  1:1603200/149‑1 (MQ=255)
    cTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAActct                                                                                                                                              >  2:1522031/1‑149 (MQ=255)
       tCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTcc                                                                                                                                           <  1:1551892/149‑1 (MQ=255)
        cAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTccc                                                                                                                                          <  1:1084231/149‑1 (MQ=255)
          gcgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCaa                                                                                                                                        >  1:2477579/1‑149 (MQ=255)
           cgcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCaaa                                                                                                                                       <  1:2505611/149‑1 (MQ=255)
            gcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAAc                                                                                                                                      >  1:2457442/1‑149 (MQ=255)
            gcCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGAGTTTAACCGCAGCCAGCACTCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAATCTCGGCACTAACTCTTCCCAAAc                                                                                                                                      <  2:2391986/149‑1 (MQ=255)
              cAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACACCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCg                                                                                                                                    <  2:3536092/149‑1 (MQ=255)
               aGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCgg                                                                                                                                   >  2:4054558/1‑149 (MQ=255)
                  ccAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGAt                                                                                                                                <  1:1708820/149‑1 (MQ=255)
                    aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                              >  2:826508/1‑149 (MQ=255)
                    aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                              >  2:2885898/1‑149 (MQ=255)
                    aaCCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGAACGCCGACAGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGc                                                                                                                              >  2:2731365/1‑149 (MQ=255)
                      ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                            >  1:3590976/1‑149 (MQ=255)
                      ccGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCgg                                                                                                                            >  2:1945425/1‑149 (MQ=255)
                        gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                          >  1:420136/1‑149 (MQ=255)
                        gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                          <  2:279627/149‑1 (MQ=255)
                        gCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCt                                                                                                                          >  1:2632346/1‑149 (MQ=255)
                           gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCGCGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                       <  2:1465413/149‑1 (MQ=255)
                           gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                       <  2:1500751/149‑1 (MQ=255)
                           gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                       <  2:1261684/149‑1 (MQ=255)
                           gtgtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCtttt                                                                                                                       <  1:2039689/149‑1 (MQ=255)
                             gtgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCa                                                                                                                     >  2:3272833/1‑149 (MQ=255)
                              tgATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGGAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGACGCGCGGCACTAACTCTTCCCAACCCGGGATGCGGCGTTTccg                                                                                                                    >  2:3439790/1‑147 (MQ=255)
                                 tCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGa                                                                                                                 <  1:2830678/149‑1 (MQ=255)
                                    aaaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAAt                                                                                                              >  1:1807474/1‑149 (MQ=255)
                                     aaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATg                                                                                                             >  1:3770002/1‑149 (MQ=255)
                                     aaaTGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATg                                                                                                             <  2:709161/149‑1 (MQ=255)
                                          gAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAg                                                                                                        >  1:1512629/1‑149 (MQ=255)
                                           aaCCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGa                                                                                                       <  2:3099388/149‑1 (MQ=255)
                                            aCCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGAt                                                                                                      >  2:1862498/1‑149 (MQ=255)
                                              cgcgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCg                                                                                                    >  1:1075026/1‑149 (MQ=255)
                                              cgcgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCg                                                                                                    >  1:3129047/1‑149 (MQ=255)
                                              cgcgTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCg                                                                                                    >  2:3743915/1‑149 (MQ=255)
                                                 gTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCa                                                                                                 <  2:1846496/149‑1 (MQ=255)
                                                 gTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCa                                                                                                 <  2:1437455/149‑1 (MQ=255)
                                                         gCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAAc                                                                                         <  1:3355461/149‑1 (MQ=255)
                                                            gccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGcc                                                                                      <  2:4114199/149‑1 (MQ=255)
                                                            gccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGcc                                                                                      >  1:2228715/1‑149 (MQ=255)
                                                            gccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGcc                                                                                      >  2:3918141/1‑149 (MQ=255)
                                                             ccagcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGccc                                                                                     <  1:1860372/149‑1 (MQ=255)
                                                                gcACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCaa                                                                                  >  1:453190/1‑149 (MQ=255)
                                                                  aCGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCTTTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAt                                                                                <  2:2007567/149‑1 (MQ=255)
                                                                  aCGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAAt                                                                                <  2:1389656/149‑1 (MQ=255)
                                                                   cGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATg                                                                               <  1:1115995/149‑1 (MQ=255)
                                                                   cGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATg                                                                               <  1:3306950/149‑1 (MQ=255)
                                                                      cGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTaa                                                                            >  1:2393710/1‑149 (MQ=255)
                                                                       gACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                           >  2:2087936/1‑149 (MQ=255)
                                                                       gACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAAt                                                                           >  1:3501106/1‑149 (MQ=255)
                                                                         ccGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATgg                                                                         >  1:2405885/1‑149 (MQ=255)
                                                                          cGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGga                                                                        <  1:1917786/149‑1 (MQ=255)
                                                                           ggAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGgag                                                                       <  1:3575296/149‑1 (MQ=255)
                                                                              aTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                                    <  1:4004706/149‑1 (MQ=255)
                                                                              aTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                                    <  2:2162676/149‑1 (MQ=255)
                                                                              aTACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCGTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATa                                                                    <  2:90640/149‑1 (MQ=255)
                                                                                aCCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTCACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATAct                                                                  <  2:3295666/149‑1 (MQ=255)
                                                                                 ccGACTGCCGTAGCAAAAATCATCGCGCAGTCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATTCATAACGCCCCAAATGTAATGGAGATActg                                                                 <  1:1160397/149‑1 (MQ=255)
                                                                                  cGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATTCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgc                                                                <  2:26867/149‑1 (MQ=255)
                                                                                  cGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgc                                                                <  1:662497/149‑1 (MQ=255)
                                                                                  cGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgc                                                                <  1:2445948/149‑1 (MQ=255)
                                                                                   gACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATActgct                                                               <  1:2664391/149‑1 (MQ=255)
                                                                                      tGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGcc                                                            <  2:2632394/149‑1 (MQ=255)
                                                                                       gCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCca                                                           >  1:383429/1‑149 (MQ=255)
                                                                                         cGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCcaca                                                         <  1:822091/149‑1 (MQ=255)
                                                                                             gCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCgg                                                     >  2:3408841/1‑149 (MQ=255)
                                                                                               aaaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                   >  1:516009/1‑149 (MQ=255)
                                                                                               aaaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGtt                                                   >  2:2636947/1‑149 (MQ=255)
                                                                                                aaaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGttt                                                  >  2:220209/1‑149 (MQ=255)
                                                                                                 aaaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTa                                                 <  1:4019693/149‑1 (MQ=255)
                                                                                                  aaTCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTAt                                                >  2:355620/1‑149 (MQ=255)
                                                                                                    tcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCt                                              <  1:253311/149‑1 (MQ=255)
                                                                                                    tcatcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCt                                              <  1:3593208/149‑1 (MQ=255)
                                                                                                     catcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTa                                             >  2:1011186/1‑149 (MQ=255)
                                                                                                      atcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTaa                                            <  1:3696923/149‑1 (MQ=255)
                                                                                                      atcGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTAAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTaa                                            <  2:3305953/149‑1 (MQ=255)
                                                                                                           gcAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGa                                       <  1:228960/149‑1 (MQ=255)
                                                                                                              gACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACTGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                                    <  2:3220248/149‑1 (MQ=255)
                                                                                                              gACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGtt                                    >  2:2224484/1‑149 (MQ=255)
                                                                                                               aCGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTc                                   <  1:1830377/149‑1 (MQ=255)
                                                                                                                   cAAGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCGTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagcc                               <  2:573021/149‑1 (MQ=255)
                                                                                                                    aaGTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagcca                              >  2:164653/1‑149 (MQ=255)
                                                                                                                      gTTCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTcagccagc                            >  1:650677/1‑149 (MQ=255)
                                                                                                                       ttCCAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCcagca                           <  2:599215/149‑1 (MQ=255)
                                                                                                                         ccAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCagcagc                         <  1:2885898/149‑1 (MQ=255)
                                                                                                                          cAGCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCt                        <  1:2730925/149‑1 (MQ=255)
                                                                                                                            gCGTGGCCTGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTgc                      <  2:509570/149‑1 (MQ=255)
                                                                                                                                    tGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACg              >  1:2343628/1‑149 (MQ=255)
                                                                                                                                    tGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACg              >  1:2117217/1‑149 (MQ=255)
                                                                                                                                    tGGAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACg              >  2:2351893/1‑149 (MQ=255)
                                                                                                                                     ggAAGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGc             >  2:3618941/1‑149 (MQ=255)
                                                                                                                                       aaGCGCGGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGctc           >  1:1380655/1‑149 (MQ=255)
                                                                                                                                           gcgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTggg       >  2:2506973/1‑149 (MQ=255)
                                                                                                                                           gcgGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTggg       >  2:3954341/1‑149 (MQ=255)
                                                                                                                                            cgGCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTgggg      <  1:3361866/149‑1 (MQ=255)
                                                                                                                                             ggCACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGga     <  1:1522031/149‑1 (MQ=255)
                                                                                                                                               cACTAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGgaga   >  2:465273/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAACCTGTCAGCGCCAGGCCAACCGCTGTGTGATCGAAAATGGAACCGCGTTTAACCGCAGCCAGCACGCCGACCGGAATACCGACTGCCGTAGCAAAAATCATCGCGCAGACGCCAAGTTCCAGCGTGGCCTGGAAGCGCGGCACGAACTCTTCCCAAACCGGGATGCGGCTTTTCATTGAAATGCCTAGATCGCCATGCATAACGCCCCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGAT  >  CP009273/3698729‑3699021

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: