Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,520,652 C→A 100% I116I (ATC→ATA yfeH → putative inorganic ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,520,6520CA98.2% 689.6 / 3.1 216I116I (ATC→ATAyfeHputative inorganic ion transporter
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (97/115);  total (97/119)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-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 216 total aligned reads displayed.

TTTGTGGGTAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATCGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGT  >  CP009273/2520508‑2520795
                                                                                                                                                |                                                                                                                                               
tttGTGGGTAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGcc                                                                                                                                              <  2:1074420/148‑1 (MQ=255)
tttGTGGGTAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGcc                                                                                                                                              <  2:1074410/148‑1 (MQ=255)
   gtgGGTAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTc                                                                                                                                           <  1:1471880/148‑1 (MQ=255)
   gtgGGTAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTc                                                                                                                                           <  1:2871940/148‑1 (MQ=255)
        tAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTc                                                                                                                                      >  1:1574481/1‑148 (MQ=255)
        tAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTc                                                                                                                                      >  1:3327096/1‑148 (MQ=255)
        tAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTc                                                                                                                                      >  1:3326992/1‑148 (MQ=255)
         aaTGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCa                                                                                                                                     <  1:2613985/148‑1 (MQ=255)
         aaTGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCa                                                                                                                                     >  1:2436738/1‑148 (MQ=255)
         aaTGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCa                                                                                                                                     <  1:894774/148‑1 (MQ=255)
          aTGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATCGCCTTCACGTCaa                                                                                                                                    <  2:2391726/148‑1 (MQ=255)
            gtgCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATg                                                                                                                                  <  1:1667376/148‑1 (MQ=255)
              gcagcaCCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATggc                                                                                                                                <  2:3347911/148‑1 (MQ=255)
                 gcaCCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATggcggg                                                                                                                             >  1:3426871/1‑148 (MQ=255)
                  caCCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATggcgggc                                                                                                                            >  2:164309/1‑148 (MQ=255)
                  caCCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATggcgggc                                                                                                                            >  1:235669/1‑148 (MQ=255)
                            cTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTcg                                                                                                                  >  1:3088957/1‑148 (MQ=255)
                            cTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTcg                                                                                                                  >  1:2978637/1‑148 (MQ=255)
                             tGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTcgc                                                                                                                 >  1:2784514/1‑148 (MQ=255)
                               tttCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCgg                                                                                                               >  1:3276956/1‑148 (MQ=255)
                               tttCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCgg                                                                                                               >  1:1613514/1‑148 (MQ=255)
                               tttCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCgg                                                                                                               <  1:2508173/148‑1 (MQ=255)
                                ttCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGc                                                                                                              >  1:2816992/1‑148 (MQ=255)
                                 tCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCa                                                                                                             >  2:2544307/1‑148 (MQ=255)
                                   cGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGc                                                                                                           >  1:613490/1‑148 (MQ=255)
                                       tCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGtt                                                                                                       >  1:2979202/1‑148 (MQ=255)
                                       tCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGtt                                                                                                       <  1:2016681/148‑1 (MQ=255)
                                       tCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGtt                                                                                                       >  2:1202016/1‑148 (MQ=255)
                                         tGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTg                                                                                                     <  2:907180/148‑1 (MQ=255)
                                           ggTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGtt                                                                                                   >  1:2436365/1‑148 (MQ=255)
                                           ggTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGtt                                                                                                   >  2:2393975/1‑148 (MQ=255)
                                           ggTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGtt                                                                                                   <  2:3048870/148‑1 (MQ=255)
                                               tACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGc                                                                                               >  2:1782309/1‑148 (MQ=255)
                                                 cTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGt                                                                                             <  2:3522116/148‑1 (MQ=255)
                                                   gTTTGCCTGGTGGACACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCg                                                                                           >  2:2375688/1‑148 (MQ=255)
                                                    tttGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCgg                                                                                          >  2:207163/1‑148 (MQ=255)
                                                    tttGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCgg                                                                                          >  2:1663376/1‑148 (MQ=255)
                                                     ttGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGc                                                                                         <  1:1947510/148‑1 (MQ=255)
                                                        ccTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCCACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATc                                                                                      <  2:1684726/148‑1 (MQ=255)
                                                            gtggAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGc                                                                                  <  1:533268/148‑1 (MQ=255)
                                                            gtggAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGc                                                                                  <  1:2370723/148‑1 (MQ=255)
                                                             tggAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGcc                                                                                 <  1:2141633/148‑1 (MQ=255)
                                                               gAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCctg                                                                               <  2:848635/148‑1 (MQ=255)
                                                                aaaCCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGc                                                                                  >  1:1578987/1‑144 (MQ=255)
                                                                aaaCCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGc                                                                                  <  2:1578137/144‑1 (MQ=255)
                                                                 aaCCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCctgct                                                                             <  1:367154/148‑1 (MQ=255)
                                                                 aaCCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCctgct                                                                             >  2:1420885/1‑148 (MQ=255)
                                                                  aCCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCctgctg                                                                            <  2:2312889/148‑1 (MQ=255)
                                                                   ccGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTgg                                                                           <  1:2148410/148‑1 (MQ=255)
                                                                   ccGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCGTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTgg                                                                           <  2:1352242/148‑1 (MQ=255)
                                                                      gTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGa                                                                        <  1:3140207/148‑1 (MQ=255)
                                                                      gTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGa                                                                        <  1:2544307/148‑1 (MQ=255)
                                                                         aaTGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGAttt                                                                     <  2:319069/148‑1 (MQ=255)
                                                                         aaTGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGAttt                                                                     <  2:1003086/148‑1 (MQ=255)
                                                                              cGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCtt                                                                <  1:478535/148‑1 (MQ=255)
                                                                              cGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCtt                                                                <  1:3415504/148‑1 (MQ=255)
                                                                               gACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCttt                                                               <  1:500324/148‑1 (MQ=255)
                                                                                aCCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTc                                                              >  2:2923108/1‑148 (MQ=255)
                                                                                aCCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTc                                                              >  2:1535982/1‑148 (MQ=255)
                                                                                   cGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCAcc                                                           <  1:473437/148‑1 (MQ=255)
                                                                                   cGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCAcc                                                           <  1:2261116/148‑1 (MQ=255)
                                                                                    gATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCa                                                          >  2:2598767/1‑148 (MQ=255)
                                                                                     aTGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCAt                                                         <  2:2352136/148‑1 (MQ=255)
                                                                                       gCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTg                                                       >  1:752840/1‑148 (MQ=255)
                                                                                        ctctACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGc                                                      >  2:1742091/1‑148 (MQ=255)
                                                                                          ctACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCtg                                                    >  2:3311395/1‑148 (MQ=255)
                                                                                            aCTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCtggt                                                  >  2:615084/1‑148 (MQ=255)
                                                                                            aCTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCtggt                                                  >  2:791242/1‑148 (MQ=255)
                                                                                             cTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCtggtt                                                 >  1:1100895/1‑148 (MQ=255)
                                                                                               ccGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCtggttgg                                               >  1:1061271/1‑148 (MQ=255)
                                                                                                    tttCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTgg                                          >  1:177393/1‑148 (MQ=255)
                                                                                                    tttCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTgg                                          >  1:346414/1‑148 (MQ=255)
                                                                                                      tctcTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGtg                                        <  1:1549711/148‑1 (MQ=255)
                                                                                                       ctctACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGtga                                       <  2:2978637/148‑1 (MQ=255)
                                                                                                        tctACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGtgat                                      <  1:121620/148‑1 (MQ=255)
                                                                                                          tACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGtgatga                                    <  2:1672761/148‑1 (MQ=255)
                                                                                                           aCTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGaa                                   >  2:1732175/1‑148 (MQ=255)
                                                                                                             ttGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATg                                 <  2:3111771/148‑1 (MQ=255)
                                                                                                              tgtgCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGt                                <  1:130046/148‑1 (MQ=255)
                                                                                                              tgtgCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGt                                <  2:3426871/148‑1 (MQ=255)
                                                                                                               gtgCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGtt                               <  1:2927213/148‑1 (MQ=255)
                                                                                                               gtgCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGtt                               <  1:1154883/148‑1 (MQ=255)
                                                                                                                tgCATTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTc                              <  2:2691832/148‑1 (MQ=255)
                                                                                                                tgCAGTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTc                              <  2:2270963/148‑1 (MQ=255)
                                                                                                                   aTTCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACg                           <  2:1035526/148‑1 (MQ=255)
                                                                                                                    ttCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACgg                          <  2:1218132/148‑1 (MQ=255)
                                                                                                                    ttCTCCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACgg                          <  1:3046308/148‑1 (MQ=255)
                                                                                                                     tctcCCGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGt                         >  2:3220467/1‑148 (MQ=255)
                                                                                                                         ccGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAg                     >  1:1867474/1‑148 (MQ=255)
                                                                                                                          cGGCTACCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAgg                    <  2:2441383/148‑1 (MQ=255)
                                                                                                                               aCCGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTGGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCa               >  2:197173/1‑148 (MQ=255)
                                                                                                                                ccGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAg              >  1:2817231/1‑148 (MQ=255)
                                                                                                                                 cGTGCAGTCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGc             >  1:2303193/1‑148 (MQ=255)
                                                                                                                                        tCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGc      >  2:2029741/1‑148 (MQ=255)
                                                                                                                                        tCTGCAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGc      >  2:1603293/1‑148 (MQ=255)
                                                                                                                                            cAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGt  <  1:1240810/148‑1 (MQ=255)
                                                                                                                                            cAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGt  <  2:2784514/148‑1 (MQ=255)
                                                                                                                                            cAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGt  >  2:1630563/1‑148 (MQ=255)
                                                                                                                                            cAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGt  <  1:2393975/148‑1 (MQ=255)
                                                                                                                                            cAATAGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGt  <  2:106636/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTTGTGGGTAATGTGCAGCACCTTCGTGCTGTTTCCGATTCTGGGTGTACTGTTTGCCTGGTGGAAACCGGTAAATGTCGACCCGATGCTCTACTCCGGTTTTCTCTACTTGTGCATTCTCCCGGCTACCGTGCAGTCTGCAATCGCCTTCACGTCAATGGCGGGCGGTAACGTCGCGGCAGCGGTTTGTTCTGCGTCGGCATCCAGCCTGCTGGGGATTTTCCTTTCACCATTGCTGGTTGGTCTGGTGATGAATGTTCACGGTGCAGGGGGCAGCCTTGAGCAGGT  >  CP009273/2520508‑2520795

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: