Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 353,610 C→G 12.5% A302G (GCG→GGG)  prpE → propionate‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913353,6100CG12.5% 227.2 / 2.1 85A302G (GCG→GGG) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (30/43);  new base G (8/4);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACG  >  NC_000913/353386‑353836
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
gtggcTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGa                                                                                                                                                                                                                                         <  1:1028198/227‑1 (MQ=255)
   gcTGGAATCCAACGAACCCTCCTGCATTCCTTCCCCTTCCGGCACGACCGGCAAACCTAAAGGCGGGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTa                                                                                                                                                                                                                                      <  2:1533885/227‑1 (MQ=255)
       gAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTa                                                                                                                                                                                                                                      <  2:1824467/223‑1 (MQ=255)
         aTCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGGGTATCGctt                                                                                                                                                                                                                                >  1:1560899/1‑227 (MQ=255)
             aaCGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACTATCGTTTTCGa                                                                                                                                                                                                                            >  1:1176809/1‑227 (MQ=255)
                       cctGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGGCATCGGTTACGAACGAGTGTCGa                                                                                                                                                                                                                 >  1:337926/1‑228 (MQ=255)
                         tGCATTCTCTACACTGCCGGCACGACCGGCAAACCTAACGGCGTGAAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGAcc                                                                                                                                                                                                              <  2:2028126/228‑1 (MQ=255)
                             ttccacacccTTCCGGGACGACCGGCAAACCTAAAGGCGTGCAGCGTGGCGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                           <  1:849098/218‑1 (MQ=255)
                             ttCTCTACCCTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGGCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                           <  1:1418330/227‑1 (MQ=255)
                                tatACATTTCCGGCACGACCGGCCAACCTAAAGGCGTGCTGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGCCCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGcc                                                                                                                                                                                                         <  2:1297303/224‑1 (MQ=255)
                                        tCCGGGACGGCCGGCCAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcggc                                                                                                                                                                                               <  1:310776/228‑1 (MQ=255)
                                        tCCGGCACGAACGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTgcggc                                                                                                                                                                                               <  2:1807705/228‑1 (MQ=255)
                                         ccGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACTATGGGTTACGAAGGGTTGCCGACCTGGCCGGACTgcgggg                                                                                                                                                                                              >  1:152539/1‑226 (MQ=255)
                                         ccGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGGATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCTCATGTTTTTGAAGGAATGCCGCCCTGGGCgggcgtgcgcg                                                                                                                                                                                              >  1:1164217/1‑219 (MQ=255)
                                          cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAGTGGGACTCTCGTTTACGAAGGATTGCCGACACTGCCCGGCTgcggcgg                                                                                                                                                                                            >  1:1826922/1‑227 (MQ=255)
                                          cGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGGGGTGGGGCATTCGGATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGACTCTCGGTTACGAAGGATTGCCGACCTGGCCgggagtcgggt                                                                                                                                                                                              >  1:1559638/1‑218 (MQ=255)
                                           ggCACGACGGGCAAACCTATAGGCGTTCAGCGTGCCGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                           <  2:1452807/228‑1 (MQ=255)
                                           ggCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGGCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGa                                                                                                                                                                                                                                         <  1:1965797/184‑1 (MQ=255)
                                            gCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACGATCGGTTACGAAGGAATTCCGAACTGGCCGGGCTgcggcgctg                                                                                                                                                                                          >  1:305344/1‑225 (MQ=255)
                                            gCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCg                                                                                                                                                                                                                                   >  1:1021286/1‑189 (MQ=255)
                                             cACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGt                                                                                                                                                                                                                                  <  2:1019426/189‑1 (MQ=255)
                                                 aCCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGGCTATCGTTTACGAAGGATTGCCGACCTGGCCCGACTGCGGCGTGtggtgg                                                                                                                                                                                    >  1:1561244/1‑228 (MQ=255)
                                                   cGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTgtggt                                                                                                                                                                                      <  1:906157/224‑1 (MQ=255)
                                                          ccTAAAGGCGTGCAGCGTGACGTCGGCGGATAGGCGGTGGCGCTGGCGCCCTCGATGGACACCATTTTGGGCGGCAAAGCGGGCAGCGTTTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAAtcgt                                                                                                                                                                            <  1:1730130/227‑1 (MQ=255)
                                                              aaGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCgaga                                                                                                                                                                       <  2:1048525/228‑1 (MQ=255)
                                                              aaGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCgag                                                                                                                                                                        <  1:1994258/227‑1 (MQ=255)
                                                                    tGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGGGGTGGGGCATTCGGATATCGTTTACGCGCCGCTGCTGGGGGGGGTGGCGGACAACGGTTACGAAGGAATGGCGACCTGGGCGGGGCGCTGCGTGTGGGTGGGCATAGTCGAGCAataat                                                                                                                                                                 >  2:1180844/1‑226 (MQ=255)
                                                                     gCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTgg                                                                                                                                                                                                           <  1:43006/187‑1 (MQ=255)
                                                                      catcGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAg                                                                                                                                                               <  2:1185060/225‑1 (MQ=255)
                                                                                   gcggATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGGGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGAGGGCGGCTATCGTTTACGAAAGATTGCCGACCTGGCCGGACTTCGGCGTGTGGTGGACAATCGTCGCGAAAAATCAGGTTAGCCggcgg                                                                                                                                                   >  1:1526978/1‑224 (MQ=255)
                                                                                         atGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGGGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACTATCGTTTACGAAGGATGTCCGACCTGGCCGGGCTGCGGCGTGTGGGGGGCGATCGTCGAGAAAAATCAGGGTAGCCggtggtttgac                                                                                                                                             >  2:144183/1‑218 (MQ=255)
                                                                                               tGGCGCTGGCGACCTCGATGGACCCCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCg                                                                                                                                                                                                        <  1:716147/164‑1 (MQ=255)
                                                                                                              ctgtGGACCCCTTTTTTGGCTGCGAAGGGGGCGGCGTGTTCTTTGCCGCATGGGATATCGGCTGGGGGGTGGGTCGTTCGATTAGGGTTGCCGCGCCGCAGCGGGCGGGGAGGGCGACTATCGTTTACGAAGGATTGCCGCCATGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATTTTCTCAGCGCCGACCGCCATTCGCGTg                                                                                                                        <  2:1607553/224‑1 (MQ=255)
                                                                                                                            tgtGGCGGCAAAGCGGGCAGAGTGTTCTTTTGCGCATCGGAGAGCGGCTGGGTGGTGGGGCATTCGTATATCGTTTCCGCGCCGCGGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                         <  2:304823/225‑1 (MQ=255)
                                                                                                                               ggcggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACt                                                                                                                                                                                                                                       >  1:1573347/1‑102 (MQ=255)
                                                                                                                                gcggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTa                                                                                                                                                                                                                                      <  2:1570556/102‑1 (MQ=255)
                                                                                                                                 cggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAAtcgtcg                                                                                                                                                                          <  1:43660/158‑1 (MQ=255)
                                                                                                                                  ggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCg                                                                                                                                                                                                                                   >  1:338266/1‑103 (MQ=255)
                                                                                                                                  ggcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTAc                                                                                                                                                                                                                              >  1:35361/1‑108 (MQ=255)
                                                                                                                                   gcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGt                                                                                                                                                                                                                                  <  2:337696/103‑1 (MQ=255)
                                                                                                                                   gcAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACg                                                                                                                                                                                                                             <  2:35295/108‑1 (MQ=255)
                                                                                                                                             ctgCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCa                                                                                       <  1:345749/225‑1 (MQ=255)
                                                                                                                                                      cTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTTCGAAGGATTGCCGACCTGGCCGGGCTGCGGCGGGTGGTGGGCAATCGTCGGGAAATATCAGGTTAGCCCGATGTTTTCAGCGCCGACCCCCCTTACAGTGGTGGAAAAAATTCCTTACGCTGAAATTTGCAAAAACGAtac                                                                             >  2:1991906/1‑226 (MQ=255)
                                                                                                                                                          tGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGCGACTTTCGTTTTCGAAGGATTGCCGACCTGGCCGGACTGCGGCGGGTGGGGGGCAATCGGCGAGAAATATCAGGGTAGCCGGATGTTTTCAGCGCCGACCGCAATTTGCGGGCTGGTAAAATTCCCTACCGCTGTAAATCGCCAACACCAtcccccg                                                                         >  2:1688086/1‑223 (MQ=255)
                                                                                                                                                            cgcATCGGATATCGGCTGGGGGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGGGGCTATCGTTTACGAAGGATTGCCGACCTGGCCgg                                                                                                                                                                                                       >  2:1427583/1‑104 (MQ=255)
                                                                                                                                                             gcATCGGGTGTGGGCTGGGTGGTGGGGCGTTCGTATATCGTTTGAGCGGCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCGGGACTGCGGCGTGTTGTGGACAATCGTCGAGAAATATCAGGTTAGc                                                                                                                                                         <  1:991974/147‑1 (MQ=255)
                                                                                                                                                             gcATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGa                                                                                                                                                                                                      <  1:1430140/104‑1 (MQ=255)
                                                                                                                                                                 cGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATGGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCGCGTCGCTgg                                                                 <  1:938102/228‑1 (MQ=255)
                                                                                                                                                                   gttttCGGCTGGGTGGTGGGCCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                               <  2:1823685/224‑1 (MQ=255)
                                                                                                                                                                   gaattCGGCTGGTTGGTGGGGCATTCGTATATCGGTTCCGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                               <  1:1318892/224‑1 (MQ=255)
                                                                                                                                                                      atggcgtGGGTGGTGGGGCATGCGTATTTCGTTTCCCCGCCCATGCGGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTg                                                           <  1:1961146/222‑1 (MQ=255)
                                                                                                                                                                      atCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGGGGGGACTATCGGTTTCGAAGGATTGCCGCCCCGGCCGGACTGCGGCGGGGGGGGGGACAACGTCGGGGAAAATCAGGGTAGCCCGGGGGTTTAGGGGCCGACCGCCATTGCCGGGCTGGAAAAATTCCCAACCCCGTGAATTCGCAAAAACGAGTGCTCCTCGCTCGGAAgg                                                            >  2:1777974/1‑227 (MQ=255)
                                                                                                                                                                          gCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGctc                                                        <  1:384547/227‑1 (MQ=255)
                                                                                                                                                                            tGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGGGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTAt                                                     <  2:1174702/228‑1 (MQ=255)
                                                                                                                                                                             gggTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTTCGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCCTTTGCGTGCTGAAAAAATTCCCTACCTCTGAAAATCGCCAACACGATCTCTCGTCTCTGGAAGAGCTCTAt                                                     >  2:199090/1‑227 (MQ=255)
                                                                                                                                                                              ggtggtGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgt                                                                                                                                                                                            >  1:394365/1‑96 (MQ=255)
                                                                                                                                                                              gggggtggGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCCGGCCGGGCTGCGGCGGGGGGGGGGCAAACGGCGGGGAAAATCAGGGTTGGCGGGTGGTTTCCGGCCCGGCCCCCCATCCCGTGCTGGAAAAAATCTCTACCGCCGGAAATTCCCACACCGCAACCCTCGGCCCTGGAAGGGCCCTac                                                     >  2:1961212/4‑226 (MQ=255)
                                                                                                                                                                               gtggtgGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                           <  2:393716/96‑1 (MQ=255)
                                                                                                                                                                                       gCATTCGTAGATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGgag                                           <  1:1995455/227‑1 (MQ=255)
                                                                                                                                                                                        cATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCCTTTCCGTGCTGAAAAAATTACCTACCGCTGAAATTCGCAAACACGAACTCTCGTCGCTGGac                                                               >  2:992278/1‑206 (MQ=255)
                                                                                                                                                                                        cATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGc                                                                                        >  1:38430/1‑183 (MQ=255)
                                                                                                                                                                                         aTTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCa                                                                                       <  2:38360/183‑1 (MQ=255)
                                                                                                                                                                                         aTTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTGACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                               <  1:994066/206‑1 (MQ=255)
                                                                                                                                                                                          ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGGTGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTTCCTACCGCTGAAATTACCGAACACGAGCTTCCGTCGCTGGAAGTGgtcccgctccttc                                              >  2:2017899/1‑211 (MQ=255)
                                                                                                                                                                                          ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGGCTATCGTTTACGAAGGATTGCCGACCCGGCCGGGCGTCGGCGGGGGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGGTGTTCTCAGCGCCGACCGCCATTCGCGGGCTGAAAAAATTTCCTACCGCTGAAATTCGCAAACACGAACTCTCGTCGCTGGGACTGATCTATc                                                    >  2:449122/1‑215 (MQ=255)
                                                                                                                                                                                           tcgttgatCTTTTACGGGCCGCTGCTGGCGGGGATGGGGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCAAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATct                                                   <  1:449857/209‑1 (MQ=255)
                                                                                                                                                                                           tCGTATATCGGTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGc                                                                     <  1:2021508/199‑1 (MQ=255)
                                                                                                                                                                                                  tCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCCGAAATTCGCAAACACGAACTCTCGTAGCTGGAAGTGCTCCATCTGGCCGGGGaaagcctggagt                               >  2:2025790/1‑218 (MQ=255)
                                                                                                                                                                                                     tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGa                                                                                                                                                     >  1:429469/1‑111 (MQ=255)
                                                                                                                                                                                                     tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGACCTCTCGTCGCTGGAAGGGCTCTATCCGGCCGGTGAAACGCTGGACGAgg                            >  2:864054/1‑227 (MQ=255)
                                                                                                                                                                                                     tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGGGTGGGGGGCAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTTGCCAAAACGAACTTCTGTCCCTGGCAGTGCTCTATCTGGCTGGAGGACCGCTGGACGAGc                            >  2:1354928/1‑228 (MQ=255)
                                                                                                                                                                                                     tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCCGGCCGGGCTGCGGCGTGTGGGGGACAATCGTCGAGAAATATCAGGTTAGCCGGGTGGTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTTCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCCCTGGAAAGGCCCTAACTGGCTGGGGAACCGCTGGACGAGc                            >  2:713135/1‑228 (MQ=255)
                                                                                                                                                                                                     tttACGCGCCGCTGCTGGCGGGGATGGCGACCATCGTTTACGAAGGATTGCCGACCTGGCCGGACTggggcgg                                                                                                                                                                                            >  2:1455317/1‑72 (MQ=255)
                                                                                                                                                                                                      ttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCgtg                                                                                                                                                                                           <  1:1457921/73‑1 (MQ=255)
                                                                                                                                                                                                      ttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGAt                                                                                                                                                    <  2:428753/111‑1 (MQ=255)
                                                                                                                                                                                                       tACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTTCCGACCTGGCCGGACTGCGGCGGGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCCGTAATTCGCAAAAACGATCTCTCGGCGCTGGAAGGGGTCCATCTGGCTTGGGAACCCCTGGACGAGcca                          >  2:1357587/1‑227 (MQ=255)
                                                                                                                                                                                                       tACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCCTTTGCGTGCTGAAAAAATTTCCTACCGCTGAAATTTGCAAAAACGATCTCCCCTCGCTTGAAGGTGCCTATCTTGCTGGGGGACCTCCGGGCGatgcc                          >  2:1136793/1‑224 (MQ=255)
                                                                                                                                                                                                               gctgctGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGc                                                          >  1:58737/1‑188 (MQ=255)
                                                                                                                                                                                                               gctgctGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGGGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGCCCGCCATTCCCGTGCTGAAAAAATTCCCTACCGCTGCAATTTGCCAACACGAACTCCCGGCGCCGGAAGTGGCCTATCTTGCTGGGGACACGCCGGCAGAGCGGAGCGcccg                  >  2:1425311/1‑226 (MQ=255)
                                                                                                                                                                                                                ctgcggGGGGGAATGGCGCTTTTCGTTTACGAAGGATTCCCGCCCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGct                                                         <  2:58631/188‑1 (MQ=255)
                                                                                                                                                                                                                   cgggggggaATGGCGGCTATCGTTTACGAAGGATTCCCGACCTGGCCGGACTGGGGCGTGTGGTGGACAATCGGCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg              <  1:1360030/223‑1 (MQ=255)
                                                                                                                                                                                                                     gggggggATGGGGACTATAGTATACGGAGGTGTGCCGCCCTGGCCGGACTGCGGGGTGTGGTGGACAATCGTCGAAAAATATGAGGTTAGCCGGATTTTCCCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGCCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt            <  2:337357/225‑1 (MQ=255)
                                                                                                                                                                                                                        ggggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGGGGGGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAgcgc                                                                                                                                         >  2:1890471/1‑104 (MQ=255)
                                                                                                                                                                                                                         gggATGGGGACTTTCGTTTACGAAAGATTGCCGACCTGGCCGGCCTGGGGCGTTTGGTGGACAATCGTCGAGAAATATCAGTTTAGCCGGTTGTTCTCAGCGcc                                                                                                                                        <  1:1893819/104‑1 (MQ=255)
                                                                                                                                                                                                                          gaaTGGCGACTACCGTTTACGAAGGATGGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCa       <  1:1357367/226‑1 (MQ=255)
                                                                                                                                                                                                                                cGACTATCGTTTACGAAGGTTTGCCGACCTGGCCGGACGGCGGCGTGTGGTGGACAATCGGCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACg  <  1:509388/227‑1 (MQ=255)
                                                                                                                                                                                                                                cGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGa                                                                >  2:46245/1‑166 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         
GTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCATTTTTGGCGGCAAAGCGGGCAGCGTGTTCTTTTGCGCATCGGATATCGGCTGGGTGGTGGGGCATTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACG  >  NC_000913/353386‑353836

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