Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,551,620 C→A 100% D226Y (GAC→TAC)  ddpF ← D,D‑dipeptide permease system, ATP‑binding component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,551,6200CA97.8% 717.7 / 6.2 227D226Y (GAC→TAC) ddpFD,D‑dipeptide permease system, ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (101/121);  total (105/122)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-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 227 total aligned reads displayed.

CAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCC  >  CP009273/1551475‑1551762
                                                                                                                                                 |                                                                                                                                              
cAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGc                                                                                                                                              >  2:36731/1‑148 (MQ=255)
cAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGc                                                                                                                                              >  1:3441151/1‑148 (MQ=255)
        cGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCt                                                                                                                                      >  2:2228256/1‑148 (MQ=255)
           cAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACa                                                                                                                                   >  2:2247742/1‑148 (MQ=255)
             aCGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAAt                                                                                                                                 <  1:1142325/148‑1 (MQ=255)
             aCGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAAt                                                                                                                                 >  1:2673734/1‑148 (MQ=255)
              cGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATc                                                                                                                                >  2:1746428/1‑148 (MQ=255)
              cGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATc                                                                                                                                >  2:3083680/1‑148 (MQ=255)
                tGCGGTTTCCCGGCAGATCCGTGTTACGTAATTCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTg                                                                                                                              >  1:638069/1‑148 (MQ=255)
                 gCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGc                                                                                                                             >  1:2444807/1‑148 (MQ=255)
                    gTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCg                                                                                                                          >  1:2045266/1‑148 (MQ=255)
                    gTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATAAATGGCGGGGAGGGAATCCAGCAATAATCGGGGGTATGGATGTGCAGGTGCCGTCCGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCg                                                                                                                          >  1:72991/1‑148 (MQ=255)
                      ttCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCgag                                                                                                                        >  2:3193533/1‑148 (MQ=255)
                      ttCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCgag                                                                                                                        >  2:2062170/1‑148 (MQ=255)
                       tCCCGTCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCgaga                                                                                                                       <  1:1810694/148‑1 (MQ=255)
                       tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCgtga                                                                                                                       <  1:3047704/148‑1 (MQ=255)
                        cccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGAt                                                                                                                      >  1:1962830/1‑148 (MQ=255)
                        cccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGAt                                                                                                                      <  2:3409077/148‑1 (MQ=255)
                        cccGGCAGATCCGTTTTACGTAATGCCCAGTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCGGCCCGAGAt                                                                                                                      <  2:2507200/148‑1 (MQ=255)
                              aGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCa                                                                                                                <  1:2370321/148‑1 (MQ=255)
                              aGATCCGTTTTAAGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGTCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTTCCCTCAGCACCTGCAGCGCGTAGCCCAGTTATACAATATGCCCGAGATacaaca                                                                                                                >  1:2793623/1‑148 (MQ=255)
                                aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGCTACATCAcc                                                                                                              >  2:3151656/1‑148 (MQ=255)
                                    gTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCa                                                                                                          >  2:613597/1‑148 (MQ=255)
                                      tttACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAcc                                                                                                        <  1:2111481/148‑1 (MQ=255)
                                        tACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCg                                                                                                      <  1:1818736/148‑1 (MQ=255)
                                         aCGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGa                                                                                                     <  2:430629/148‑1 (MQ=255)
                                          cGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGAt                                                                                                    >  1:2467725/1‑148 (MQ=255)
                                            tAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                                                                                  <  1:3256798/148‑1 (MQ=255)
                                            tAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCg                                                                                                  <  2:1271717/148‑1 (MQ=255)
                                              aTGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCt                                                                                                >  2:2902984/1‑148 (MQ=255)
                                                gCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                                                                                              <  1:1178208/148‑1 (MQ=255)
                                                  ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCata                                                                                            <  1:1056552/148‑1 (MQ=255)
                                                  ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCata                                                                                            <  1:962126/148‑1 (MQ=255)
                                                     ttcttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGa                                                                                         >  1:3177020/1‑148 (MQ=255)
                                                     ttcttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGa                                                                                         <  2:2438233/148‑1 (MQ=255)
                                                        ttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGt                                                                                      <  2:2045266/148‑1 (MQ=255)
                                                              cAGCGGTTTATCAATGGCGGGGAGGTAAACCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTAAAATCTGCCCGAGATACATCACCGCCACACGATCGCTCATATGACTTATCAca                                                                                >  1:1424665/1‑147 (MQ=255)
                                                              cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCAcc                                                                                <  2:966427/148‑1 (MQ=255)
                                                               aGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCg                                                                               >  1:2709008/1‑148 (MQ=255)
                                                               aGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCg                                                                               >  1:2595471/1‑148 (MQ=255)
                                                               aGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCg                                                                               >  1:2276284/1‑148 (MQ=255)
                                                                 cGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgag                                                                             >  2:2894314/1‑148 (MQ=255)
                                                                 cGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCgag                                                                             <  1:1375757/148‑1 (MQ=255)
                                                                     ttATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt                                                                         <  2:2556064/148‑1 (MQ=255)
                                                                     ttATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt                                                                         <  2:1148159/148‑1 (MQ=255)
                                                                       aTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtg                                                                       <  2:2074141/148‑1 (MQ=255)
                                                                       aTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtg                                                                       <  2:2114217/148‑1 (MQ=255)
                                                                       aTCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtg                                                                       <  2:342935/148‑1 (MQ=255)
                                                                         cAAGGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtg                                                                     <  2:1860962/148‑1 (MQ=255)
                                                                           aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgtg                                                                   >  1:681519/1‑148 (MQ=255)
                                                                           aTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgtgtg                                                                   >  2:1230669/1‑148 (MQ=255)
                                                                             ggCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGaa                                                                 <  1:1347572/148‑1 (MQ=255)
                                                                             ggCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGaa                                                                 <  1:36728/148‑1 (MQ=255)
                                                                              gCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGCCATTGTGTGaaa                                                                >  1:2700300/1‑148 (MQ=255)
                                                                               cGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAAt                                                               <  1:911830/148‑1 (MQ=255)
                                                                               cGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAAt                                                               <  1:2399951/148‑1 (MQ=255)
                                                                                ggggAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATc                                                              <  1:1621928/148‑1 (MQ=255)
                                                                                 gggagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCa                                                             >  2:1476522/1‑148 (MQ=255)
                                                                                   gagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGc                                                           >  2:2531319/1‑148 (MQ=255)
                                                                                    agggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGca                                                          <  1:1096132/148‑1 (MQ=255)
                                                                                    agggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGca                                                          <  2:1739444/148‑1 (MQ=255)
                                                                                        aaTCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATa                                                      <  2:3442165/148‑1 (MQ=255)
                                                                                           ccAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGt                                                   >  2:690915/1‑148 (MQ=255)
                                                                                            cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTc                                                  <  1:1856886/148‑1 (MQ=255)
                                                                                             aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCa                                                 <  1:2848097/148‑1 (MQ=255)
                                                                                              gCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAg                                                >  1:1547416/1‑148 (MQ=255)
                                                                                               cAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGc                                               >  2:1673996/1‑148 (MQ=255)
                                                                                               cAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGc                                               >  2:3411684/1‑148 (MQ=255)
                                                                                                aataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGcc                                              <  1:3324114/148‑1 (MQ=255)
                                                                                                  taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCg                                            >  1:454694/1‑148 (MQ=255)
                                                                                                   aatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGt                                           >  1:3442142/1‑148 (MQ=255)
                                                                                                    atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTg                                          >  1:2940050/1‑148 (MQ=255)
                                                                                                     tCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa                                         >  2:1913981/1‑148 (MQ=255)
                                                                                                     tCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa                                         <  1:195913/148‑1 (MQ=255)
                                                                                                      cGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAt                                        <  1:2228022/148‑1 (MQ=255)
                                                                                                        ggTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAttt                                      <  2:1384789/148‑1 (MQ=255)
                                                                                                         gtgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGAtttt                                     >  1:3031917/1‑148 (MQ=255)
                                                                                                          tgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTc                                    <  2:909692/148‑1 (MQ=255)
                                                                                                          tgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTc                                    <  2:1129209/148‑1 (MQ=255)
                                                                                                           gtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTcc                                   <  1:2062170/148‑1 (MQ=255)
                                                                                                            tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCt                                  >  1:320561/1‑148 (MQ=255)
                                                                                                             aTGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTg                                 <  2:2803277/148‑1 (MQ=255)
                                                                                                                gATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGTCCGTGATTTTCCTGTAg                              <  2:1729831/148‑1 (MQ=255)
                                                                                                                gATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAg                              <  1:1630445/148‑1 (MQ=255)
                                                                                                                        ggTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCa                      >  1:2950275/1‑148 (MQ=255)
                                                                                                                         gTGCCGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAg                     <  1:401217/148‑1 (MQ=255)
                                                                                                                            ccGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTaa                  >  2:1013414/1‑148 (MQ=255)
                                                                                                                             cGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTaaa                 <  1:613597/148‑1 (MQ=255)
                                                                                                                             cGTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTaaa                 <  2:2193203/148‑1 (MQ=255)
                                                                                                                              gTCAGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAt                <  1:3411684/148‑1 (MQ=255)
                                                                                                                                 aGCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGa             <  2:2677054/148‑1 (MQ=255)
                                                                                                                                  gCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAg            >  1:1694557/1‑148 (MQ=255)
                                                                                                                                  gCACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAg            >  1:2220692/1‑148 (MQ=255)
                                                                                                                                   cACCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAgg           >  2:2861488/1‑148 (MQ=255)
                                                                                                                                    aCCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCGGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGa          >  2:1000021/1‑148 (MQ=255)
                                                                                                                                    aCCTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGa          <  2:3022867/148‑1 (MQ=255)
                                                                                                                                     ccTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGAt         <  1:1794644/148‑1 (MQ=255)
                                                                                                                                     ccTGCTGCGCGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGAt         <  2:3262555/148‑1 (MQ=255)
                                                                                                                                           gcgcGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTgcgc   <  1:2084322/148‑1 (MQ=255)
                                                                                                                                           gcgcGTAGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTgcgc   <  2:1812946/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATTGAGGATCTGCGCC  >  CP009273/1551475‑1551762

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: ‑