Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,369 A→C 15.6% T80P (ACC→CCC)  proX → glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,807,3690AC15.6% 217.5 / 4.9 74T80P (ACC→CCC) proXglycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (54/7);  new base C (2/11);  total (56/18)
Fisher's exact test for biased strand distribution p-value = 4.78e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAACCGGTTGTA  >  NC_000913/2807134‑2807603
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
gCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCt                                                                                                                                                                                                                                         >  1:628787/1‑239 (MQ=255)
 cGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCt                                                                                                                                                                                                                                         >  2:122829/1‑238 (MQ=255)
       aGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCAccg                                                                                                                                                                                                                                   >  2:607118/1‑238 (MQ=255)
         cGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCAA‑CACCTCGCTCGCTTCCGGCG‑ATGCAACCTTCAcc                                                                                                                                                                                                                                    >  1:907067/1‑235 (MQ=255)
           tACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAACCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTg                                                                                                                                                                                                                              <  1:22746/239‑1 (MQ=255)
           tACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTg                                                                                                                                                                                                                              >  2:1075542/1‑239 (MQ=255)
              tttttGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCACCCTTCACCGCCGTGa                                                                                                                                                                                                                             >  1:338785/1‑237 (MQ=255)
              tttttGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGa                                                                                                                                                                                                                             >  2:882597/1‑237 (MQ=255)
                tttGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTg                                                                                                                                                                                                                         >  1:347391/1‑239 (MQ=255)
                  tGCGACAGCTTTTGCCACGCTTATCTCTACACAAACTTTTGCGGCCGATCTGCCGGGAAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGA‑TTACAACTTTGGCTC‑CCCCTCGCTTGCTTCCGGCG‑ATGAACCCTTCACCGCCGTGAACTGGa                                                                                                                                                                                                                       <  1:1276681/239‑1 (MQ=255)
                  tGCGACAGCTTTTGCCACGCTTATCTCAACAAAAACTTTTGCGCCCGATCTGCGGGGCAAAGGCATTCCTGTTAATCAATTTCAGAGCCCCATCACTAAAGAACCCTCCCAGACCCTGCTGGTCAGTCGTGCGCTGGAAAAATTGGGTTATACCGTCAACAACCCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CCCCTCGCTTGCTTCCGGCG‑ATGCACCCTTCACCGCCGTGAACTGGa                                                                                                                                                                                                                       <  2:1016732/239‑1 (MQ=255)
                   gCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGAGCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGACCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGCGATATTAGGCTATACCGTCACCAAACCCAGAGAAGTAGA‑TTACAACGATGGCTA‑CACCTCGCCTTCTTCCGGCG‑ATGCCCCCTTCACCGCCGTGAACTGGa                                                                                                                                                                                                                       >  2:1194570/1‑238 (MQ=255)
                     gACAGCTTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAACCCTTCCAGACGCGGCTGGTCAGGCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTCCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGc                                                                                                                                                                                                                    <  1:1105976/239‑1 (MQ=255)
                       cggCGTTTGCCACGCTTATCCCTACACAAATTTTTGCGGCCGATCTCCCGGGCAAAGGCATTACTGTTAACCCAGTCAAGAGCACCATCACAGAAGAACCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGAAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGA‑TTCCAACGTTGGCTA‑CCCCTCGCTTGCTTCCGGCG‑ATGCACCCTTCACCGCCGTGAACTGGACGCCa                                                                                                                                                                                                                  <  2:869047/237‑1 (MQ=255)
                                 cACGCTTACCTCTACACAAACTTTTGCGGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATCCCGTAAACAACCCCAGCGAAGTAGT‑TTACAACGTTGGCTA‑CACCTCCCTTGCTTCCGGCG‑ATGAACCCTTCACCGCCGTGAACTGGACGCCACTGCATGac                                                                                                                                                                                                         <  1:117228/238‑1 (MQ=255)
                                                  aaaCTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGcc                                                                                                                                                                                       >  2:1279474/1‑239 (MQ=255)
                                                   aaCTTTCGCCGCCGAGGTGCCGGGCAAAGCCATGATTCAGAATCCAGCTAAGAACACCACCCAGAAAGAACCCTGCCAGACGCTGCGGGTCACTCGCACGATGGAGAAATTAGTTTATCCCGTCAAAAAACCCGGCAAAGTAGA‑TTACAACGTGGGCAG‑CACCCCCCTTGCTCCCGGCG‑AGGCACCTTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCg                                                                                                                                                                                      <  2:861009/239‑1 (MQ=255)
                                                    aCTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCgg                                                                                                                                                                                     >  1:378586/1‑239 (MQ=255)
                                                              ccGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAg                                                                                                                                                                           >  2:1163841/1‑239 (MQ=255)
                                                                        gggCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATc                                                                                                                                                                 >  1:497174/1‑239 (MQ=255)
                                                                                gCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAggg                                                                                                                                                         >  2:613228/1‑239 (MQ=255)
                                                                                    tACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTAt                                                                                                                                                     >  1:212850/1‑239 (MQ=255)
                                                                                    tACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTAt                                                                                                                                                     >  1:503027/1‑239 (MQ=255)
                                                                                         ttAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGtt                                                                                                                                                >  1:695857/1‑239 (MQ=255)
                                                                                                   tCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGgcg                                                                                                                                         >  2:245180/1‑236 (MQ=255)
                                                                                                   tCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAAAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGgcg                                                                                                                                         <  1:245180/236‑1 (MQ=255)
                                                                                                    cAGACCACCATCACTGAAAAAACCTCCCAGACGCTGCGGGTCAGTCGTGCGCGGGAGAAATTAGGTTATCCCTTCAAAAACCCCAGCGAAGTAGA‑TTACACCGTTGGCTA‑CACCTCGCTTGCTCCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  1:527072/239‑1 (MQ=255)
                                                                                                       agCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGcaca                                                                                                                                  >  2:1142047/1‑239 (MQ=255)
                                                                                                       agCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGcaca                                                                                                                                  >  2:1308810/1‑239 (MQ=255)
                                                                                                            cATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGtt                                                                                                                             >  2:683140/1‑239 (MQ=255)
                                                                                                              tCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTa                                                                                                                            >  1:1046519/1‑238 (MQ=255)
                                                                                                                         aCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTACCGGCGCGGCACAGGGTTACCTGATCGATaa                                                                                                                >  1:104226/1‑239 (MQ=255)
                                                                                                                         aCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATaa                                                                                                                >  1:438598/1‑239 (MQ=255)
                                                                                                                            ttCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGaa                                                                                                             >  1:1000371/1‑239 (MQ=255)
                                                                                                                              ccAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGaaa                                                                                                            >  1:793431/1‑238 (MQ=255)
                                                                                                                                    gctgctGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAAccgccg                                                                                                     >  2:315120/1‑239 (MQ=255)
                                                                                                                                      tgctgGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCCCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTCACGGCGCGGCACCGGGTTACCTGATCGATAAGCAAACCGCCGa                                                                                                    >  2:366645/1‑238 (MQ=255)
                                                                                                                                       gcgggTCAGTCGTGCGAGGGAGAAATTAGTTTTCACGGACAACACCCACGGCGAAGAAAA‑TACCAGCGTTGCCTA‑CCCCGCGTTTCCGCCGGCCG‑ACGCACCCTACACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGAc                                                                                                   <  2:628787/235‑1 (MQ=255)
                                                                                                                                                tCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAATGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGAACaaaa                                                                                         >  2:439238/1‑239 (MQ=255)
                                                                                                                                                tCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACaaaa                                                                                         >  2:1338707/1‑239 (MQ=255)
                                                                                                                                                     cgcTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCAcc                                                                                    >  1:989990/1‑239 (MQ=255)
                                                                                                                                                     cgcTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCAcc                                                                                    >  2:199702/1‑239 (MQ=255)
                                                                                                                                                          gagaAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACa                                                                                >  2:555183/1‑238 (MQ=255)
                                                                                                                                                              aaTTAGGTTATACCGTCAAAAACCCCAGGGAAGTAAA‑TTAAACCGTT‑GCTACCCCCCCCCTTGCTCCCGGCG‑ATCCACCCTTCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGca                                                                           <  1:882597/239‑1 (MQ=255)
                                                                                                                                                                    gTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTa                                                                                                                                                      >  2:987459/1‑158 (MQ=255)
                                                                                                                                                                    gTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTa                                                                                                                                                      <  1:987459/158‑1 (MQ=255)
                                                                                                                                                                        taCCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCa                                                                                                                                                                                                             >  1:309689/1‑99 (MQ=255)
                                                                                                                                                                              caacaaACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGa                                                           >  2:1130795/1‑239 (MQ=255)
                                                                                                                                                                              caacaaACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGa                                                           >  1:976007/1‑239 (MQ=255)
                                                                                                                                                                              caacaaACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGa                                                           >  2:696762/1‑239 (MQ=255)
                                                                                                                                                                               acaacccccAGCAAGGTAAA‑TTACACCTGTGTCTA‑CCCCCCGTTTTCTCCCGGCG‑ATCCACCCTCCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGaa                                                          <  1:122829/233‑1 (MQ=255)
                                                                                                                                                                               aacaacCCCAGCAAAG‑AGAGCTAAAACGTTGGCTC‑CACCGCGCTTTCTGCCGGCG‑AGGAACCCTCCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGaa                                                          <  1:1318864/239‑1 (MQ=255)
                                                                                                                                                                                acaaACCCCGGGAAGTAAATTTACAACTTTGCCCA‑CCCCCCGTTTCCTCCCGGGA‑ATGAACCCTCCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGaa                                                          <  1:1075542/239‑1 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGaa                                                                   >  1:1004136/1‑220 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCaa                                                 >  2:204065/1‑238 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCaa                                                 >  1:217509/1‑238 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCaa                                                 >  2:779408/1‑238 (MQ=255)
                                                                                                                                                                                               tAGT‑TTACAACTTGGGCTA‑CCCCCCGCTTGCTTCCGGCG‑AGGCACCTTTCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  2:786534/148‑1 (MQ=255)
                                                                                                                                                                                               tAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTc                                          >  2:130516/1‑239 (MQ=255)
                                                                                                                                                                                                        aaCGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCa                                   >  1:338041/1‑238 (MQ=255)
                                                                                                                                                                                                          cGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTAc                                                                                                                           >  1:786534/1‑148 (MQ=255)
                                                                                                                                                                                                          cGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATACGAAGATCGCCAAACTGTTCGATAcacacg                                >  1:1297865/1‑234 (MQ=255)
                                                                                                                                                                                                             tGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGACAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGCAAGATCCGAAGATCGCCAAACTGTCCGATACCAACGGCg                             >  2:492371/1‑239 (MQ=255)
                                                                                                                                                                                                             tGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGACAAACTGTTCGATACCAACGGCg                             >  1:152665/1‑239 (MQ=255)
                                                                                                                                                                                                              ggCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCg                             >  1:4581/1‑238 (MQ=255)
                                                                                                                                                                                                                    cacCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGaa                       >  2:449164/1‑239 (MQ=255)
                                                                                                                                                                                                                     acCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCACAGGCGACGGaaa                      >  2:784702/1‑239 (MQ=255)
                                                                                                                                                                                                                      ccTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  1:767450/127‑1 (MQ=255)
                                                                                                                                                                                                                      ccTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     >  2:767450/1‑127 (MQ=255)
                                                                                                                                                                                                                       cTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGaaaa                     >  1:433230/1‑238 (MQ=255)
                                                                                                                                                                                                                          gcttgcttCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAAc                                 >  1:1094193/1‑223 (MQ=255)
                                                                                                                                                                                                                                 tCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCAAAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTaa           >  2:341191/1‑238 (MQ=255)
                                                                                                                                                                                                                                  ccGGCGTATCCAACCTTCCCCGCCGTGACCTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTTTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAAc          <  2:217509/239‑1 (MQ=255)
                                                                                                                                                                                                                                          tGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAACCGGTTGTa  >  2:654265/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
GCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGA‑TTACAACGTTGGCTA‑CACCTCGCTTGCTTCCGGCG‑ATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAACCGGTTGTA  >  NC_000913/2807134‑2807603

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