Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 627,428 C→A 100% I693I (ATC→ATA)  cstA → carbon starvation protein involved in peptide utilization; APC peptide transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273627,4280CA98.6% 682.3 / 0.4 218I693I (ATC→ATA) cstAcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (112/102);  total (115/103)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 218 total aligned reads displayed.

TGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTT  >  CP009273/627284‑627570
                                                                                                                                                |                                                                                                                                              
 gCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTgg                                                                                                                                            >  2:2185291/1‑148 (MQ=255)
     ggTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCa                                                                                                                                        >  1:2717013/1‑148 (MQ=255)
       ttAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAgg                                                                                                                                      >  2:1563999/1‑148 (MQ=255)
        tAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGc                                                                                                                                     >  2:179540/1‑148 (MQ=255)
        tAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGc                                                                                                                                     >  2:1956947/1‑148 (MQ=255)
        tAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGc                                                                                                                                     >  2:391058/1‑148 (MQ=255)
             aTCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAg                                                                                                                                <  1:110682/148‑1 (MQ=255)
              tcttctTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAgg                                                                                                                               >  1:1278637/1‑148 (MQ=255)
              tcttctTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAgg                                                                                                                               >  1:1881274/1‑148 (MQ=255)
                            gtggtggtgGTTCTGGCACTTTTCTCTATTAAGACGGCACTTGCGGCATTGAAAGATCCGACGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTc                                                                                                                 <  1:1200249/145‑1 (MQ=255)
                                tggtggTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCagag                                                                                                             <  2:90993/148‑1 (MQ=255)
                                 ggtggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGt                                                                                                            <  2:2699533/148‑1 (MQ=255)
                                 ggtggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGt                                                                                                            >  2:1972607/1‑148 (MQ=255)
                                 ggtggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGt                                                                                                            >  2:2402858/1‑148 (MQ=255)
                                  gtggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTg                                                                                                           <  2:3347514/148‑1 (MQ=255)
                                  gtggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTg                                                                                                           >  2:2958207/1‑148 (MQ=255)
                                  gtggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTg                                                                                                           <  1:1249795/148‑1 (MQ=255)
                                   tggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGa                                                                                                          <  1:1437819/148‑1 (MQ=255)
                                      ttCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggg                                                                                                       <  2:711373/148‑1 (MQ=255)
                                         tGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggggg                                                                                                    <  1:2151768/148‑1 (MQ=255)
                                           gCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGa                                                                                                  >  1:2584527/1‑148 (MQ=255)
                                            cACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGAt                                                                                                 >  1:1049947/1‑148 (MQ=255)
                                             aCTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATg                                                                                                <  1:2462993/148‑1 (MQ=255)
                                              cTGTTCTCGTTTATGACGTCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAATATGTCGTGGAGATAGTGTCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGt                                                                                               <  2:579565/148‑1 (MQ=255)
                                               tGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGtt                                                                                              <  2:1849538/148‑1 (MQ=255)
                                                 ttCTCTATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCTAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTgg                                                                                            <  1:335165/148‑1 (MQ=255)
                                                  tctcGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGc                                                                                           >  1:3215801/1‑148 (MQ=255)
                                                    tcGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGa                                                                                         >  2:3082917/1‑148 (MQ=255)
                                                     cGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGaa                                                                                        <  2:2182002/148‑1 (MQ=255)
                                                     cGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGaa                                                                                        <  1:1625418/148‑1 (MQ=255)
                                                     cGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGaa                                                                                        >  1:1638767/1‑148 (MQ=255)
                                                     cGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGaa                                                                                        <  2:3304673/148‑1 (MQ=255)
                                                      gATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAAt                                                                                       <  2:203339/148‑1 (MQ=255)
                                                       aTTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTCAAGTCAGAGTGAGGGGGGGATGTTGGCGACTg                                                                                      >  1:2949186/1‑148 (MQ=255)
                                                        ttAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGt                                                                                     >  2:2976356/1‑148 (MQ=255)
                                                        ttAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGt                                                                                     <  1:3399756/148‑1 (MQ=255)
                                                         tAAGACGGCTCTTGCGGCATTTAAAGATCAGAAGCAAACGTCGAAATAAAAGCCGTATGAACCAATGCCTGAAAATGTCGAGGAGATTGTGGCGCTGTCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGtt                                                                                    <  2:1252146/148‑1 (MQ=255)
                                                          aaGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTg                                                                                   <  2:117061/148‑1 (MQ=255)
                                                           aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTgg                                                                                  >  1:2869200/1‑148 (MQ=255)
                                                                gCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAg                                                                             <  1:3019649/148‑1 (MQ=255)
                                                                gCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAg                                                                             <  2:730972/148‑1 (MQ=255)
                                                                 cACTTGCGGCATTGAAATATCCGAAGCC‑ACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTg                                                                           <  2:618559/148‑1 (MQ=255)
                                                                 cACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGt                                                                            >  2:1498463/1‑148 (MQ=255)
                                                                    ttGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGcc                                                                         >  2:2499161/1‑148 (MQ=255)
                                                                    ttGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGcc                                                                         >  1:988875/1‑148 (MQ=255)
                                                                     tGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGccc                                                                        >  2:210281/1‑148 (MQ=255)
                                                                       cGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAg                                                                      <  1:3042056/148‑1 (MQ=255)
                                                                          cATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGAACACTAAAGTCAGAGGGTGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                                                                   >  2:1983567/1‑148 (MQ=255)
                                                                              gAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                               >  2:1711115/1‑148 (MQ=255)
                                                                               aaaGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGACTGTTGGCTTAGTGCCCAGGGTTCCct                                                              >  2:1613366/1‑148 (MQ=255)
                                                                               aaaGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCct                                                              >  1:820539/1‑148 (MQ=255)
                                                                                aaGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGATGAGATAGTGGCGCAGGCAAAAGTCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                                             <  2:1287328/148‑1 (MQ=255)
                                                                                aaGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                                             <  1:2188990/148‑1 (MQ=255)
                                                                                 aGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctct                                                            >  1:1676828/1‑148 (MQ=255)
                                                                                  gATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                           >  2:110785/1‑148 (MQ=255)
                                                                                    tCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAc                                                         >  1:1071699/1‑148 (MQ=255)
                                                                                    tCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAc                                                         <  2:238659/148‑1 (MQ=255)
                                                                                     ccGAAGCCAACGGCGAAAGACACGCCGTATGAACCNATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCNCACTAAAGTCAGAGTGAGGGGGGGATGTTNGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcn                                                        <  1:137436/148‑2 (MQ=255)
                                                                                        aaGCCAACGGCTAAATAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                                     <  1:349630/148‑1 (MQ=255)
                                                                                        aaGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                                     >  2:3269195/1‑148 (MQ=255)
                                                                                          gCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc                                                   >  1:945414/1‑148 (MQ=255)
                                                                                          gCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc                                                   >  1:2292762/1‑148 (MQ=255)
                                                                                           ccAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAcc                                                  <  2:1349494/148‑1 (MQ=255)
                                                                                             aaCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                                <  1:2055185/148‑1 (MQ=255)
                                                                                              aCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCctc                                               <  1:643750/148‑1 (MQ=255)
                                                                                                ggCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCtctc                                             >  2:2550050/1‑148 (MQ=255)
                                                                                                 gCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcc                                            <  2:297696/148‑1 (MQ=255)
                                                                                                     accacaCGCCGTATGAACCAATGCCGGACAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGAGAGGGGGGGATGTTGGCGAAGGTTGGCTTAGTGCCCAGGGTTCCCTCTCAGCCTAACCCTCTCCCCgg                                        >  2:3029216/6‑148 (MQ=255)
                                                                                                     aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCgg                                        >  1:199266/1‑148 (MQ=255)
                                                                                                      agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGt                                       >  1:1728629/1‑148 (MQ=255)
                                                                                                      agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGt                                       >  2:2445933/1‑148 (MQ=255)
                                                                                                         aaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTggg                                    >  2:325728/1‑148 (MQ=255)
                                                                                                          aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgggg                                   <  1:555688/148‑1 (MQ=255)
                                                                                                           cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGc                                  >  1:2831470/1‑148 (MQ=255)
                                                                                                            gCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                 <  1:2656895/148‑1 (MQ=255)
                                                                                                             ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGa                                >  2:2598492/1‑148 (MQ=255)
                                                                                                               gTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                              >  2:2586542/1‑148 (MQ=255)
                                                                                                               gTATGAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                              >  1:1020228/1‑148 (MQ=255)
                                                                                                                 aTGAACCAATGCCGGAATATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgggg                            <  1:2028446/148‑1 (MQ=255)
                                                                                                                   gAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                          <  1:1498463/148‑1 (MQ=255)
                                                                                                                   gAACCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                          <  1:2652489/148‑1 (MQ=255)
                                                                                                                    aaCCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACt                         >  2:1720843/1‑148 (MQ=255)
                                                                                                                     aCCAATGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg                        >  2:605745/1‑148 (MQ=255)
                                                                                                                      cccaTGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                       <  2:3405809/145‑1 (MQ=255)
                                                                                                                        aaTGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAACGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAcc                     <  2:2868188/148‑1 (MQ=255)
                                                                                                                         aTGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg                    >  1:3074777/1‑148 (MQ=255)
                                                                                                                         aTGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg                    >  1:1043673/1‑148 (MQ=255)
                                                                                                                          tGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                   >  1:2454670/1‑148 (MQ=255)
                                                                                                                          tGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                   >  2:1388496/1‑148 (MQ=255)
                                                                                                                          tGCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                   >  2:1650342/1‑148 (MQ=255)
                                                                                                                           gCCGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg                  >  1:889994/1‑148 (MQ=255)
                                                                                                                            ccGGAAAATGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgc                 >  1:2492152/1‑148 (MQ=255)
                                                                                                                                aaaaTGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt             >  2:424516/1‑148 (MQ=255)
                                                                                                                                 aaaTGTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt            >  1:1760202/1‑148 (MQ=255)
                                                                                                                                     gTCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATa        >  1:2326554/1‑148 (MQ=255)
                                                                                                                                      tCGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg       >  2:2997437/1‑148 (MQ=255)
                                                                                                                                       cGAGGAGATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGc      >  2:863406/1‑148 (MQ=255)
                                                                                                                                         aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGTGGGGACTGACCGAGCGCGTTGATAGCAt    >  2:845577/1‑148 (MQ=255)
                                                                                                                                         aggagATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAt    <  1:1682445/148‑1 (MQ=255)
                                                                                                                                          ggagATAGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAtt   >  1:1954480/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTT  >  CP009273/627284‑627570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑