Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,367 A→C 16.8% 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,3670AC16.8% 220.1 / 6.9 80T80P (ACC→CCC) proXglycine betaine/proline ABC transporter periplasmic binding protein, cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (54/11);  new base C (2/13);  total (56/24)
Fisher's exact test for biased strand distribution p-value = 5.96e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAACCGGTTGTA  >  NC_000913/2807132‑2807601
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
gCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCt                                                                                                                                                                                                                                         >  1:628787/1‑239 (MQ=255)
 cGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCt                                                                                                                                                                                                                                         >  2:122829/1‑238 (MQ=255)
  gcatgaccttatttttgGCGACACCTTTTGCAACCTTTATTTCTACACAACCTTTGGCTGCCGATCTCCCGGGCAAAGTCTTTACTGAAATTCCAGTTCAGAGCCCCATCACTAAAGAAACCTCCCAGACGCTGCTGGTCAGTGGTGCGCTGGAAAAATTAGGTATACCGTCAAAAAACCCAGCGAAGTAGATTACAACGTTGGCTCCCCCTCGCTTGCTTCCGGCGATGCAACCTGCa                                                                                                                                                                                                                                      <  2:509877/227‑1 (MQ=255)
       aGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCAccg                                                                                                                                                                                                                                   >  2:607118/1‑238 (MQ=255)
         cGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCAACACCTCGCTCGCTTCCGGCGATGCAACCTTCAcc                                                                                                                                                                                                                                    >  1:907067/1‑235 (MQ=255)
           tACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAACCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTg                                                                                                                                                                                                                              <  1:22746/239‑1 (MQ=255)
           tACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTg                                                                                                                                                                                                                              >  2:1075542/1‑239 (MQ=255)
              tttttGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGa                                                                                                                                                                                                                             >  1:338785/1‑237 (MQ=255)
              tttttGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGa                                                                                                                                                                                                                             >  2:882597/1‑237 (MQ=255)
                tttGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTg                                                                                                                                                                                                                         >  1:347391/1‑239 (MQ=255)
                  tGCGACAGCTTTTGCCACGCTTATCTCTACACAAACTTTTGCGGCCGATCTGCCGGGAAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGATTACAACTTTGGCTCCCCCTCGCTTGCTTCCGGCGATGAACCCTTCACCGCCGTGAACTGGa                                                                                                                                                                                                                       <  1:1276681/239‑1 (MQ=255)
                  tGCGACAGCTTTTGCCACGCTTATCTCAACAAAAACTTTTGCGCCCGATCTGCGGGGCAAAGGCATTCCTGTTAATCAATTTCAGAGCCCCATCACTAAAGAACCCTCCCAGACCCTGCTGGTCAGTCGTGCGCTGGAAAAATTGGGTTATACCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGa                                                                                                                                                                                                                       <  2:1016732/239‑1 (MQ=255)
                   gCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGAGCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGACCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGCGATATTAGGCTATACCGTCACCAAACCCAGAGAAGTAGATTACAACGATGGCTACACCTCGCCTTCTTCCGGCGATGCCCCCTTCACCGCCGTGAACTGGa                                                                                                                                                                                                                       >  2:1194570/1‑238 (MQ=255)
                     gACAGCTTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAACCCTTCCAGACGCGGCTGGTCAGGCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTCCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGc                                                                                                                                                                                                                    <  1:1105976/239‑1 (MQ=255)
                       cggCGTTTGCCACGCTTATCCCTACACAAATTTTTGCGGCCGATCTCCCGGGCAAAGGCATTACTGTTAACCCAGTCAAGAGCACCATCACAGAAGAACCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGAAGAAATTAGGTTATACCGTCAACAACCCCAGCGAAGTAGATTCCAACGTTGGCTACCCCTCGCTTGCTTCCGGCGATGCACCCTTCACCGCCGTGAACTGGACGCCa                                                                                                                                                                                                                  <  2:869047/237‑1 (MQ=255)
                                 cccgcttaccttcccacaACCTTTTGCTGCCGATCGGCCGGAAAAAGGATTTACTGTTAATCCAGTTCAGAGCACCATCAGTGAAGAACCCTCCCAGACGCGGCGGGTCAGTCGTGCGCGGGAGAATTTAGGTTATCCCGTCAACAACCCCAGCGAAGTAGATTACACCGTTGGCTACCCCTCGCTTGCTCCGGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGaca                                                                                                                                                                                                        <  2:1012635/225‑1 (MQ=255)
                                 cACGCTTACCTCTACACAAACTTTTGCGGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATCCCGTAAACAACCCCAGCGAAGTAGTTTACAACGTTGGCTACACCTCCCTTGCTTCCGGCGATGAACCCTTCACCGCCGTGAACTGGACGCCACTGCATGac                                                                                                                                                                                                         <  1:117228/238‑1 (MQ=255)
                                                  aaaCTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGcc                                                                                                                                                                                       >  2:1279474/1‑239 (MQ=255)
                                                   aaCTTTCGCCGCCGAGGTGCCGGGCAAAGCCATGATTCAGAATCCAGCTAAGAACACCACCCAGAAAGAACCCTGCCAGACGCTGCGGGTCACTCGCACGATGGAGAAATTAGTTTATCCCGTCAAAAAACCCGGCAAAGTAGATTACAACGTGGGCAGCACCCCCCTTGCTCCCGGCGAGGCACCTTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCg                                                                                                                                                                                      <  2:861009/239‑1 (MQ=255)
                                                    aCTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCgg                                                                                                                                                                                     >  1:378586/1‑239 (MQ=255)
                                                              ccGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAg                                                                                                                                                                           >  2:1163841/1‑239 (MQ=255)
                                                                        gggCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATc                                                                                                                                                                 >  1:497174/1‑239 (MQ=255)
                                                                                gCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAggg                                                                                                                                                         >  2:613228/1‑239 (MQ=255)
                                                                                    tACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTAt                                                                                                                                                     >  1:212850/1‑239 (MQ=255)
                                                                                    tACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTAt                                                                                                                                                     >  1:503027/1‑239 (MQ=255)
                                                                                         ttAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGtt                                                                                                                                                >  1:695857/1‑239 (MQ=255)
                                                                                                   tCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGgcg                                                                                                                                         >  2:245180/1‑236 (MQ=255)
                                                                                                   tCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAAAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGgcg                                                                                                                                         <  1:245180/236‑1 (MQ=255)
                                                                                                    cagacaccaacccctaaaaaaaCCTCCCAGCCTCTGCTGGCCAGTCGTGCCGTGGAAAAATTAGTTTATCCCGTCAACAACCCCAGCGAAGTAGATTACACCTTTGCCTCCCCCCCGCTTGCTCCCGGCGAGGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  2:998436/220‑1 (MQ=255)
                                                                                                    cAGACCACCATCACTGAAAAAACCTCCCAGACGCTGCGGGTCAGTCGTGCGCGGGAGAAATTAGGTTATCCCTTCAAAAACCCCAGCGAAGTAGATTACACCGTTGGCTACACCTCGCTTGCTCCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  1:527072/239‑1 (MQ=255)
                                                                                                     agaccaccaccccaaaaaaaaCCTCCCAGCTGCTGGTGGTCGGCGGTGCGCTGGAAAAATTAGGTTATACCGTCAACAACCCCGGCGAAGTAGATTACACGGTTGGCAACACCCCGCTTGTTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGca                                                                                                                                    <  2:567791/221‑1 (MQ=255)
                                                                                                       agCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGcaca                                                                                                                                  >  2:1308810/1‑239 (MQ=255)
                                                                                                       agCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGcaca                                                                                                                                  >  2:1142047/1‑239 (MQ=255)
                                                                                                            cATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGtt                                                                                                                             >  2:683140/1‑239 (MQ=255)
                                                                                                              tCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTa                                                                                                                            >  1:1046519/1‑238 (MQ=255)
                                                                                                                         aCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTACCGGCGCGGCACAGGGTTACCTGATCGATaa                                                                                                                >  1:104226/1‑239 (MQ=255)
                                                                                                                         aCCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATaa                                                                                                                >  1:438598/1‑239 (MQ=255)
                                                                                                                            ttCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGaa                                                                                                             >  1:1000371/1‑239 (MQ=255)
                                                                                                                              ccAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGaaa                                                                                                            >  1:793431/1‑238 (MQ=255)
                                                                                                                                    gctgctGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAAccgccg                                                                                                     >  2:315120/1‑239 (MQ=255)
                                                                                                                                      tgctgGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCCCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTCACGGCGCGGCACCGGGTTACCTGATCGATAAGCAAACCGCCGa                                                                                                    >  2:366645/1‑238 (MQ=255)
                                                                                                                                       gcgggTCAGTCGTGCGAGGGAGAAATTAGTTTTCACGGACAACACCCACGGCGAAGAAAATACCAGCGTTGCCTACCCCGCGTTTCCGCCGGCCGACGCACCCTACACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGAc                                                                                                   <  2:628787/235‑1 (MQ=255)
                                                                                                                                            cagggcgtgagcggaaaaaaTTATGTTATCCCATCAAAAACCCCGGGGATGAGATTAACATGTTGGCCGCCCCCTCGCTTCTTGCGGGCGGGGAACCTTTCCCCGCGGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTa                                                                                              <  2:719720/221‑1 (MQ=255)
                                                                                                                                                tCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAATGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGAACaaaa                                                                                         >  2:439238/1‑239 (MQ=255)
                                                                                                                                                tCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACaaaa                                                                                         >  2:1338707/1‑239 (MQ=255)
                                                                                                                                                     cgcTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCAcc                                                                                    >  1:989990/1‑239 (MQ=255)
                                                                                                                                                     cgcTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCAcc                                                                                    >  2:199702/1‑239 (MQ=255)
                                                                                                                                                          gagaAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACa                                                                                >  2:555183/1‑238 (MQ=255)
                                                                                                                                                            atggtgtaaaaggaacacaaacaacCCCATCAATGTAGATTACACCGTGGCCCACCCCCGCTTGCCGTCGGCGGAGGAACCCTCCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCTTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCg                                                                             <  2:1004136/220‑1 (MQ=255)
                                                                                                                                                              aaTTAGGTTATACCGTCAAAAACCCCAGGGAAGTAAATTAAACCGTTGCTACCCCCCCCCTTGCTCCCGGCGATCCACCCTTCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGca                                                                           <  1:882597/239‑1 (MQ=255)
                                                                                                                                                                    gTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTa                                                                                                                                                      >  2:987459/1‑158 (MQ=255)
                                                                                                                                                                    gTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTa                                                                                                                                                      <  1:987459/158‑1 (MQ=255)
                                                                                                                                                                        taCCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCa                                                                                                                                                                                                             >  1:309689/1‑99 (MQ=255)
                                                                                                                                                                              caacaaACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGa                                                           >  1:976007/1‑239 (MQ=255)
                                                                                                                                                                              caacaaACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGa                                                           >  2:1130795/1‑239 (MQ=255)
                                                                                                                                                                              caacaaACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGa                                                           >  2:696762/1‑239 (MQ=255)
                                                                                                                                                                               acaacccccAGCAAGGTAAATTACACCTGTGTCTACCCCCCGTTTTCTCCCGGCGATCCACCCTCCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGaa                                                          <  1:122829/233‑1 (MQ=255)
                                                                                                                                                                               aacaacCCCAGCAAAGAGAGCTAAAACGTTGGCTCCACCGCGCTTTCTGCCGGCGAGGAACCCTCCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGaa                                                          <  1:1318864/239‑1 (MQ=255)
                                                                                                                                                                                acaaACCCCGGGAAGTAAATTTACAACTTTGCCCACCCCCCGTTTCCTCCCGGGAATGAACCCTCCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGaa                                                          <  1:1075542/239‑1 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGaa                                                                   >  1:1004136/1‑220 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCaa                                                 >  2:204065/1‑238 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCaa                                                 >  1:217509/1‑238 (MQ=255)
                                                                                                                                                                                         gCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCaa                                                 >  2:779408/1‑238 (MQ=255)
                                                                                                                                                                                               tAGTTTACAACTTGGGCTACCCCCCGCTTGCTTCCGGCGAGGCACCTTTCCCCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  2:786534/148‑1 (MQ=255)
                                                                                                                                                                                               tAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTc                                          >  2:130516/1‑239 (MQ=255)
                                                                                                                                                                                                        aaCGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCa                                   >  1:338041/1‑238 (MQ=255)
                                                                                                                                                                                                          cGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTAc                                                                                                                           >  1:786534/1‑148 (MQ=255)
                                                                                                                                                                                                          cGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATACGAAGATCGCCAAACTGTTCGATAcacacg                                >  1:1297865/1‑234 (MQ=255)
                                                                                                                                                                                                             tGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGACAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGCAAGATCCGAAGATCGCCAAACTGTCCGATACCAACGGCg                             >  2:492371/1‑239 (MQ=255)
                                                                                                                                                                                                             tGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGACAAACTGTTCGATACCAACGGCg                             >  1:152665/1‑239 (MQ=255)
                                                                                                                                                                                                              ggCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCg                             >  1:4581/1‑238 (MQ=255)
                                                                                                                                                                                                                     cacCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGaa                       >  2:449164/1‑239 (MQ=255)
                                                                                                                                                                                                                      acCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCACAGGCGACGGaaa                      >  2:784702/1‑239 (MQ=255)
                                                                                                                                                                                                                       ccTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     >  2:767450/1‑127 (MQ=255)
                                                                                                                                                                                                                       ccTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGc                                                                                                                                     <  1:767450/127‑1 (MQ=255)
                                                                                                                                                                                                                        cTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGaaaa                     >  1:433230/1‑238 (MQ=255)
                                                                                                                                                                                                                           gcttgcttCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAAc                                 >  1:1094193/1‑223 (MQ=255)
                                                                                                                                                                                                                                  tCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCAAAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTaa           >  2:341191/1‑238 (MQ=255)
                                                                                                                                                                                                                                   ccGGCGTATCCAACCTTCCCCGCCGTGACCTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTTTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAAc          <  2:217509/239‑1 (MQ=255)
                                                                                                                                                                                                                                            tGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAACCGGTTGTa  >  2:654265/1‑239 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          
GCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACCATCACTGAAGAAACCTTCCAGACGCTGCTGGTCAGTCGTGCGCTGGAGAAATTAGGTTATACCGTCAACAAACCCAGCGAAGTAGATTACAACGTTGGCTACACCTCGCTTGCTTCCGGCGATGCAACCTTCACCGCCGTGAACTGGACGCCACTGCATGACAACATGTACGAAGCTGCCGGTGGCGATAAGAAATTTTATCGTGAAGGGGTATTTGTTAACGGCGCGGCACAGGGTTACCTGATCGATAAGAAAACCGCCGACCAGTACAAAATCACCAACATCGCACAACTGAAAGATCCGAAGATCGCCAAACTGTTCGATACCAACGGCGACGGAAAAGCGGATTTAACCGGTTGTA  >  NC_000913/2807132‑2807601

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: