Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,561:1 +G 6.3% coding (850/999 nt) fecC ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,5591.G 6.3% 251.3 / 12.0 64coding (852/999 nt)fecCferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

GTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAG  >  NC_000913/4512293‑4512832
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                 
gtgtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCCCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGCACGATCTCCGGGGAATGCCAGCCCGCGTGCCAGTACATCTGCCCGAAGCCTCAACATGGCCCCCAGCACCAAGCCCACTGGCAGTAc                                                                                                                                                                                                                                                                           >  1:220890/1‑275 (MQ=255)
  gtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATATTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCACCACTGCGCCTGCGGACAGATCTCCGGGGAAGGCCAGCGCGCGAGCGAGTACATCTGCCAGCAGCCTCAGCGTCGCCCCCAGCAGCATGCTCACTGGCAGTACg                                                                                                                                                                                                                                                                          >  2:373327/1‑274 (MQ=255)
   tCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGGGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGt                                                                                                                                                                                                                                                                         >  2:96093/1‑274 (MQ=255)
       ncnTGCCGGTGGGGCAGCAGCGCGCGCAACGGCACGGGGATCACTCCCATATGGAGGGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTgcgc                                                                                                                                                                                                                                                                    <  1:303948/274‑1 (MQ=255)
                                  cacggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCCCTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCACCAGCATGCCCACTGGCAGTACGTTCCGCTGATCGAAGCCCGCCCACAAgcccccc                                                                                                                                                                                                                                         >  1:180587/1‑270 (MQ=255)
                                  cacggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCCTCAGCGTGGCCCCCAGCAGCATGCCCACTGGCCGTCCGTTGCGCTGATCGAAGCCTGCCCAGAAgccccc                                                                                                                                                                                                                                          >  1:161431/1‑270 (MQ=255)
                                            ggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGAAAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACCTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTCCGTTGCGCTGATCGAAGCCTGCCCAGATGCGCGCCAAATGTGGCa                                                                                                                                                                                                                               >  2:145826/1‑275 (MQ=255)
                                                   tCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTCCGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACcagca                                                                                                                                                                                                                         >  2:35294/1‑274 (MQ=255)
                                                      cATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAcc                                                                                                                                                                                                                     >  2:257807/1‑275 (MQ=255)
                                                        tatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGTGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGAAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTCATCGAAGCCTGCCCAGAAGCGCACCAGCTGTGGCACCAGCCGACCg                                                                                                                                                                                                                    >  2:382424/1‑274 (MQ=255)
                                                          tGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGCTGTGGCACCAGCAGACCGAt                                                                                                                                                                                                                  >  2:172277/1‑274 (MQ=255)
                                                           ggAGTGATAACAGCGCACAGCCCGCCTGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGGACGTCGCGCTGATCGAACCCTGCCCAGAACCGCCCCAGATGTGGCACCAGCAGACCGa                                                                                                                                                                                                                   >  2:153429/1‑272 (MQ=255)
                                                             aGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATaaa                                                                                                                                                                                                               >  1:277968/1‑274 (MQ=255)
                                                             aGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCAAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGCAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATaca                                                                                                                                                                                                               >  2:355531/1‑272 (MQ=255)
                                                             aGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGACAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATaaa                                                                                                                                                                                                               >  2:245848/1‑274 (MQ=255)
                                                                 aTAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACGGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACACCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCCGCCCAGAAGCGCGCCAGCTGCGGTCCCAGCAGACCGATAAACGcc                                                                                                                                                                                                           >  1:281839/1‑274 (MQ=255)
                                                                   aaCAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCAcc                                                                                                                                                                                                        <  1:340524/275‑1 (MQ=255)
                                                                     cgacggacaGCCCGCCGGGGCAAGGTGGATGAAAATAACCAGCGCAATTTCCATCCTGCCCTCCTCCCAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCgg                                                                                                                                                                                                      <  1:145826/269‑1 (MQ=255)
                                                                             aGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACACGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAACAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTCCACGGAAGCGCGCCCGAAGTGGCCACAGCAGACCGATAAACACCACCGGACCTGCt                                                                                                                                                                                               >  1:214625/1‑274 (MQ=255)
                                                                                 cGCCAGGGAAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTGGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCt                                                                                                                                                                                          <  1:373327/275‑1 (MQ=255)
                                                                                         cAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACCCCACCGGCCCTGCTACGCTGccgaacc                                                                                                                                                                                  >  1:242543/1‑274 (MQ=255)
                                                                                                gatgaAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCCCCAGATGTGGCACCAGCAGACCGAACAACGCCACCGGACCTGCTACGCTGacgacggccccaac                                                                                                                                                                           >  1:77994/1‑264 (MQ=255)
                                                                                                 atgaAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCCTACGCACGCcccacc                                                                                                                                                                           >  1:47155/1‑272 (MQ=255)
                                                                                                     aaaTAACCAGCGCAATTTTCATCCTCGCCTCCTCACAATCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGCTCGAAGCCTGCCCAGAAGCGCGCCAGAt                                                                                                                                                                                                                                     <  2:260152/212‑1 (MQ=255)
                                                                                                       aTAACCAGCGCAATTTTCATCGCCGCCTCCTACCAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAag                                                                                                                                                                     <  2:146797/274‑1 (MQ=255)
                                                                                                        tAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAagc                                                                                                                                                                    <  2:277968/274‑1 (MQ=255)
                                                                                                             aGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGGGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCAc                                                                                                                                                               <  2:281839/274‑1 (MQ=255)
                                                                                                               cgcAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGCTGTGGCACCAGCAGACCGATAAACGCCCCCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTa                                                                                                                                                             >  2:253115/1‑275 (MQ=255)
                                                                                                               cgcAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGACCGCCCCCACCCACAGCAGCACTa                                                                                                                                                             >  2:186991/1‑275 (MQ=255)
                                                                                                               cgcAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTaa                                                                                                                                                            >  2:168751/1‑275 (MQ=255)
                                                                                                                gcAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGACAAACCCCACCGGACCTGCCACGCTGACGCACGCCCCAACCAGAACCAGCCCAAAc                                                                                                                                                           >  1:14947/1‑275 (MQ=255)
                                                                                                                gcAATTTCCTCCCTCGCCTCCTCACAAGACAAGAAAAGCACGGGCAGCTAATCAGCGCAAGCACTGCGCTTGCGGGCAGATCTCCGGGGAAGGCCAGCGGGCGTTCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGAATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAAc                                                                                                                                                           <  1:245848/275‑1 (MQ=255)
                                                                                                                       tCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtga                                                                                                                                                    <  2:77994/275‑1 (MQ=255)
                                                                                                                            cccGCCCCCCCACAAGCCAGACAAAGAAAGGGCGGCCAATCAGCGCCAGCACTGCCCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGTGAGTACCACTGCCAGCAGCGTCAGCGTTGCCCCCAGCAGCAGGCTCACTGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGAc                                                                                                                                                <  1:253115/273‑1 (MQ=255)
                                                                                                                                      cacaAGCCAGTCAAAGCAAGGGCGGCCAAACAGAGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGAGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGTGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGcc                                                                                                                                     <  1:153429/275‑1 (MQ=255)
                                                                                                                                                accaaGCAAGGGCTGCCAATCAGCGCCATCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGtt                                                                                                                           <  2:161431/272‑1 (MQ=255)
                                                                                                                                                aCAAAGCAGGGGCTGCCAAGCAGGGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGtt                                                                                                                           <  2:242543/275‑1 (MQ=255)
                                                                                                                                                aCAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGtt                                                                                                                           <  2:180587/275‑1 (MQ=255)
                                                                                                                                                  aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCAGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCa                                                                                                                         <  2:328821/275‑1 (MQ=255)
                                                                                                                                                      cAAGGGCTGCCAATCATCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                     <  2:47155/275‑1 (MQ=255)
                                                                                                                                                      cAAGGGCTGCCAATCAGAGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCTTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                     <  2:220890/275‑1 (MQ=255)
                                                                                                                                                      cAAGGGCGGCCAGGCAGCGCCAGCGCTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                     <  1:355531/275‑1 (MQ=255)
                                                                                                                                                      cAAGGGCGGCCAATCCGCGCCAGCACTGCGCCTGCGGGCAGGTCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                     <  1:96093/275‑1 (MQ=255)
                                                                                                                                                       aaGGGCTGCCAACCAGCGCCCGCACTGCGCCGGCGGGCAGATCTCCGGGGAAGCCCAGCGCGCGTGCGATTACATCGGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                     <  1:186991/275‑1 (MQ=255)
                                                                                                                                                              gCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGCCAGGTTCACTCCCAGCGTa                                                                                                              >  2:357146/1‑274 (MQ=255)
                                                                                                                                                                 agtCAGCGCCAGCACTGGGCCGGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATggg                                                                                                          <  1:257807/273‑1 (MQ=255)
                                                                                                                                                                         ccAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGACGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCt                                                                                                   <  2:14947/274‑1 (MQ=255)
                                                                                                                                                                                 gcgcCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCACCCAGAACCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCAACCTCAGCGTATGGGCGGTGCTGGCGCTgcc                                                                                          >  1:169575/1‑273 (MQ=255)
                                                                                                                                                                                 gcgcCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCAc                                                                                                                        >  1:366613/1‑245 (MQ=255)
                                                                                                                                                                                     cTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                                                                                                                    <  2:366613/245‑1 (MQ=255)
                                                                                                                                                                                       gCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCaccaacca                                                                                                                                                                         >  1:9023/1‑190 (MQ=255)
                                                                                                                                                                                           gCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAag                                                                                                                                                                     <  2:9023/190‑1 (MQ=255)
                                                                                                                                                                                                         aaGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCg                                                                   >  2:96443/1‑274 (MQ=255)
                                                                                                                                                                                                                gcgcgcGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCcagcag                                                            <  1:96443/274‑1 (MQ=255)
                                                                                                                                                                                                                    gcgTGCGAGTAAATCGGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCCCTGGCAGTACGTTGCGCTGATCGAAGCCGGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAac                                                        <  1:382424/274‑1 (MQ=255)
                                                                                                                                                                                                                        ggcggTACGTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTAGGTTGCGCGGATCGAGGCCTGCCCAGAAGCTGGCCAGATTTGGAACCAGCCAACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGAc                                                   <  1:35294/271‑1 (MQ=255)
                                                                                                                                                                                                                                   tGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATCTTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGCGGTTGCACCGGTTCAGTTGATTCGCCAGCAGCAACCAGACAGGGACTGCa                                         >  2:226074/1‑274 (MQ=255)
                                                                                                                                                                                                                                        ctagcacgaggttGGCCCCAAGCAGCATGCTAACTGGACATACGTTGGGCTAATCGAAGCCTCCCCGGAAGCGCGCCAGATGTGGCACCCGCAGACCAATCAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCCGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCGGCGTATGGGCGGTGCTGTCGCTGAGTTTGAGCAGGTTCAGTTGTTTAGCCAGCAGCAACACGACAGGGACTGCAGTGac                                    <  1:84141/262‑1 (MQ=255)
                                                                                                                                                                                                                                          agcaTCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACata                                                                                                                                                        >  1:113571/1‑156 (MQ=255)
                                                                                                                                                                                                                                            caTCCGCGTGGCCCCCAGCAGAATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGTTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGaccacca                               <  2:365451/275‑1 (MQ=255)
                                                                                                                                                                                                                                              tCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtga                                                                                                                                                    <  2:113571/156‑1 (MQ=255)
                                                                                                                                                                                                                                                    tgacccccAGGAGCAGGCTCACTGGCAGTACGTTGGGCTGATCTAAGCCTGCCCAGAAGCGAGCCAGGTGTGGCACCAGCAGACCGATAAACGCCACCGGGCCTGCTACGCTGACGCACGCACCACCCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAag                       <  1:172277/272‑1 (MQ=255)
                                                                                                                                                                                                                                                            agcagcATGATCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCCACACGACCGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCa               >  2:87536/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                   tGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGTCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATcc        >  2:233067/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          tGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCACGAGCTGCCAGACATCCTGCCAg  >  2:349059/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                 
GTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAG  >  NC_000913/4512293‑4512832

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