Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,767,640 T→C 5.6% N349S (AAC→AGC)  ydiJ ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,767,6400TC 5.6% 234.0 / 3.9 71N349S (AAC→AGC) ydiJputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base C (2/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

CAATATAATCCGCCAGGTGTTCCGGCGGTACGCAGGTATCCTCAGCAAACGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACTTTAGGCAAGCGTGTAATATCCAGCCGCGCTTCGGTAATAAAGGCCAGCGTCCCT  >  NC_000913/1767377‑1767891
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
atcaataCTCCGCCAGGTGTCCCGCGGGTACGAAGGTTCCCTCAGCAACCGGAATGGGCTTGGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCg                                                                                                                                                                                                                                                     <  2:519120/268‑1 (MQ=255)
 aaTATAATCCGCCAGGTGTTCCGGCGGTACGCAGGTAGCCGCAGCAAACGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGccc                                                                                                                                                                                                                                                      <  1:586663/270‑1 (MQ=255)
 aaTATAATCCGCCAGGTGGTCCGGCGGTACGCAGGTATCCTCAGCAAATGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCg                                                                                                                                                                                                                                                     <  1:578510/271‑1 (MQ=255)
                                         tCAGCAAACGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGccc                                                                                                                                                                                                                                                      >  1:587605/1‑230 (MQ=255)
                                         tCAGCAAACGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCa                                                                                                                                                                                                                                                 >  1:828567/1‑235 (MQ=255)
                                          cAGCAAACGGAATTGGCTTAGCGGCACCTTTTGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGccc                                                                                                                                                                                                                                                      <  1:535481/229‑1 (MQ=255)
                                          cAGCAAACGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCAt                                                                                                                                                                                                                                                <  2:828566/235‑1 (MQ=255)
                                              aaaCGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCGTAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAACACCTCTTGATCAGCCACCTCCGCAAACaaaacccccacccc                                                                                                                                                                                                         >  2:1028099/1‑258 (MQ=255)
                                                  ggAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGcc                                                                                                                                                                                                    <  1:638431/271‑1 (MQ=255)
                                                   gAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGATTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACCATGTTCAGCC‑‑AGCATCCCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGa                                                                                                                                                                                                 >  1:115443/1‑271 (MQ=255)
                                                               cGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACAACTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGCGTGCATTTACCCGCCCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAAATCGCTGACGGAATGCAAGACAATAACtac                                                                                                                                                                                       >  1:1101566/1‑269 (MQ=255)
                                                                     ctggagCAGTGCCAACCCGACCACCGGATTTTTTGCGCATTGCATCGATACCTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCC‑‑AGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAg                                                                                                                                                                              <  2:115443/267‑1 (MQ=255)
                                                                      tttGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTACTCAACTCGCt                                                                                                                                                                                                               >  1:141113/1‑240 (MQ=255)
                                                                            aTTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGcc                                                                                                                                                                                                                                                       >  2:467799/1‑194 (MQ=255)
                                                                             ttGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGccc                                                                                                                                                                                                                                                      <  1:467800/194‑1 (MQ=255)
                                                                                  aaGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAAACAACTCGCTGACGGAATGCCAGACAAAATATTCCCCCGCCCGacacaacacc                                                                                                                                                                    >  2:875520/1‑262 (MQ=255)
                                                                                           aaCGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTCCCCGCTCATCAATCAGCGCCTCATCATCACCCGCCAATTCCACAATGTTCACCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCCGACAatat                                                                                                                                                                                           >  1:428415/1‑239 (MQ=255)
                                                                                             cGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGtt                                                                                                                                                                                                                                                            >  2:418871/1‑172 (MQ=255)
                                                                                             cGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGcac                                                                                                                                                                                                                                 >  1:653621/1‑199 (MQ=255)
                                                                                              ggCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTc                                                                                                                                                                                                                                                           <  1:418872/172‑1 (MQ=255)
                                                                                                cTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGcc                                                                                                                                                                                                                                                       >  1:561127/1‑174 (MQ=255)
                                                                                                 ttttttGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGccc                                                                                                                                                                                                                                                      <  2:561126/174‑1 (MQ=255)
                                                                                                                            aaCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGccc                                                                                                                                                                                                                                                      >  2:62787/1‑147 (MQ=255)
                                                                                                                            aaCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAAt                                                                                                                                                                                                       >  1:1031764/1‑194 (MQ=255)
                                                                                                                            aaCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCa                                                                                                                                                                                                   >  2:1113286/1‑198 (MQ=255)
                                                                                                                            aaCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGCCAATATCTTCCCGCGCCAGATTCCGCCCTTTTGAGTCACCCGTCCCCACCGAAACGCCCCGCGCCCCCAccc                                                                                                                         >  2:405223/1‑271 (MQ=255)
                                                                                                                             aCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCg                                                                                                                                                                                                                                                     <  1:62787/147‑1 (MQ=255)
                                                                                                                             aCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATg                                                                                                                                                                                                      <  2:1031763/194‑1 (MQ=255)
                                                                                                                             aCCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAg                                                                                                                                                                                                  <  1:1113287/198‑1 (MQ=255)
                                                                                                                                     cAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACACGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTACACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCTGCGCCAGATTCAGCACTTTTGAGTCAACCGTCACTACCGACAGACCCCACGCCCCAACAATAAAAggag                                                                                                                >  1:926914/1‑270 (MQ=255)
                                                                                                                                                acacCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCt                                                                                                                                                                                                               >  2:1002545/1‑166 (MQ=255)
                                                                                                                                                          gCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGc                                                                                                                                                                                                                                   >  2:594235/1‑136 (MQ=255)
                                                                                                                                                          gCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTc                                                                                                                                                                          >  1:604984/1‑193 (MQ=255)
                                                                                                                                                          gCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGCTCAGccc                                                                                                                                                                                                                                                      >  1:772858/1‑117 (MQ=255)
                                                                                                                                                           ccGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGca                                                                                                                                                                                                                                  <  1:594236/136‑1 (MQ=255)
                                                                                                                                                           ccGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGCTCAGCCCg                                                                                                                                                                                                                                                     <  2:772857/117‑1 (MQ=255)
                                                                                                                                                                 acacCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGTCGCGTTACGCAGCTCAGAATCaaaaaaag                                                                                     >  1:1003316/1‑267 (MQ=255)
                                                                                                                                                                 acacCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGCCGCGt                                                                                                             >  1:893217/1‑247 (MQ=255)
                                                                                                                                                                        cTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCcgc                                                                                                                                   >  1:5163/1‑218 (MQ=255)
                                                                                                                                                                           gatgaCTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGAt                                                                                                                                                                            >  2:583453/1‑174 (MQ=255)
                                                                                                                                                                            atgaCTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGAtt                                                                                                                                                                           <  1:583454/174‑1 (MQ=255)
                                                                                                                                                                                      aTCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGc                                                                                                                                                                                                     >  2:572543/1‑138 (MQ=255)
                                                                                                                                                                                      aTCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAACCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGAc                                                                <  1:1028100/271‑1 (MQ=255)
                                                                                                                                                                                       tCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGcc                                                                                                                                                                                                    <  1:572544/138‑1 (MQ=255)
                                                                                                                                                                                                  aaCCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCcgc                                                                                                                                   >  1:1068553/1‑192 (MQ=255)
                                                                                                                                                                                                   aCCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCcgcg                                                                                                                                  <  2:1068552/192‑1 (MQ=255)
                                                                                                                                                                                                        gcACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCtt                                                                                                                                                                                                                                           >  1:935110/1‑82 (MQ=255)
                                                                                                                                                                                                         cacaGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTg                                                                                                                                                                                                                                          <  2:935109/82‑1 (MQ=255)
                                                                                                                                                                                                            agagTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGAt                                                                                                                                                                            >  1:293812/1‑141 (MQ=255)
                                                                                                                                                                                                             gagTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGAtt                                                                                                                                                                           <  2:293811/141‑1 (MQ=255)
                                                                                                                                                                                                                 gCATTTACCCGCTCATCAATCAGCGCCTTATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGa                                                                                                                                                                                                         >  2:797991/1‑107 (MQ=255)
                                                                                                                                                                                                                 gCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATa                                                                                                                       >  1:856282/1‑189 (MQ=255)
                                                                                                                                                                                                                  cATTTACCCGCTCATCAATCAGCGCCTTATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGaa                                                                                                                                                                                                        <  1:797992/107‑1 (MQ=255)
                                                                                                                                                                                                                  cATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATaa                                                                                                                      <  2:856281/189‑1 (MQ=255)
                                                                                                                                                                                                                               atcaatcaGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGa                                                                                                                                                                                                         >  2:624902/1‑93 (MQ=255)
                                                                                                                                                                                                                                tcaatcaGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGaa                                                                                                                                                                                                        <  1:624903/93‑1 (MQ=255)
                                                                                                                                                                                                                                         gcCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACTTTAGGCAAGCGTGTAATATCCAGCCGCGCTTCGGTAATAAAGGc            <  1:405224/272‑1 (MQ=255)
                                                                                                                                                                                                                                         gcCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACTTTAGGCAAGCGTGTAATATCCAGCCGCGCTTCGGTAATAAAGGc            <  1:1099881/272‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGCGTCAAAAGCGTCCTAATTGACGTTACCCAGACGGAGCCCCTCAGGCAAGCGTGTAAGAGCCAGACGCCCCTCgcaaaacacacac           >  2:838863/1‑260 (MQ=255)
                                                                                                                                                                                                                                          ccTCATCATCACCAGCAAATTCCACAATGCTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCt                                                                                                                                                                                                               >  1:990120/1‑76 (MQ=255)
                                                                                                                                                                                                                                           cTCATCATCACCAGCAAATTCCACAATGCTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTg                                                                                                                                                                                                              <  2:990119/76‑1 (MQ=255)
                                                                                                                                                                                                                                            ctatgccccccAGAAAACTCAACAAGGTTCGGCCCCGGTATCTCGTGGTCAGGCTAACCGGGATTCAACTCGCTGACGGACTGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTGGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACTTTAGGCAAGCGTGTAATATCCAGCCGCGCTTCGGTAATAAAGGCCAg         <  2:926913/263‑1 (MQ=255)
                                                                                                                                                                                                                                                catcaCCAGCAAATTCCACAATGTTCATCCCGAGCATCTCTTGGTAAGGCACCTCGGTAATAAACTCGCTGACGGAAGGACAGACAATATCTTCCCGCGCCAGATTCAGCACTTTCACTTATACAGTCTGTACCAAAAGCGCCagc                                                                                                                                   >  1:562052/1‑146 (MQ=255)
                                                                                                                                                                                                                                                 atcaCCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAg                                                                                                                                                              >  2:964194/1‑118 (MQ=255)
                                                                                                                                                                                                                                                 atcaCCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACttt                                                      >  2:131731/1‑222 (MQ=255)
                                                                                                                                                                                                                                                 atcaCCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACttt                                                      <  1:131731/222‑1 (MQ=255)
                                                                                                                                                                                                                                                  tcaCCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGt                                                                                                                                                             <  1:964195/118‑1 (MQ=255)
                                                                                                                                                                                                                                                   caCCAGCAAATTCCACATTTTTCAGCCCGAGCCTCTCTTGGTCAGGCACATCGGTAATCAACTCGCGGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACTTTAGGCAAGCGTGTAATATCCAGCCGCGCTTCGGTAATAAAGGCCAGCGTCCCt  <  1:609079/272‑1 (MQ=255)
                                                                                                                                                                                                                                                     ccAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCcg                                                                                                                                                                                    >  2:86109/1‑92 (MQ=255)
                                                                                                                                                                                                                                                      cAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCcgc                                                                                                                                                                                   <  1:86109/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acaATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCAc                                                                     >  1:679611/1‑191 (MQ=255)
                                                                                                                                                                                                                                                                  caATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCAcc                                                                    <  2:679610/191‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
CAATATAATCCGCCAGGTGTTCCGGCGGTACGCAGGTATCCTCAGCAAACGGAATTGGCTTAGCGGCACCTTTGGCATTGCCAAGCAGACCAACGGCTTTTTTGCGCATTGCATAGATACGTTCAACCCCCGCCAGCTCACGGCACACCTGCCAGCCGATGACACCTGCTTGATGACTGGCGATCAGCTCATCCAACCGCGCACAGAGTGCATTTACCCGCTCATCAATCAGCGCCTCATCATCACCAGCAAATTCCACAATGTTCAGCCCGAGCATCTCTTGGTCAGGCACATCGGTAATCAACTCGCTGACGGAATGCCAGACAATATCTTCCCGCGCCAGATTCAGCACTTTTGAGTCCACCGTCTCTACCGAAAGCGCCCGCGCCTCAACCATAAACGGCGCGTTACGCAGCGCAGAGTCAAAAGAGTCATATTTGACGTTCACCAGACGGCGCACTTTAGGCAAGCGTGTAATATCCAGCCGCGCTTCGGTAATAAAGGCCAGCGTCCCT  >  NC_000913/1767377‑1767891

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: