Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,515,348 A→C L444V (TTG→GTG)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,3480AC98.6% 228.3 / NA 71L444V (TTG→GTG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (33/37);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCG  >  NC_000913/4515102‑4515585
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
gaagcgtggcggtataGCTCGCTTCTGCGTGCGTGCCTGGGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCAGCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGGATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAg                                                                                                                                                                                                                                            <  2:235290/239‑1 (MQ=255)
 aaGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGc                                                                                                                                                                                                                                           >  1:391210/1‑251 (MQ=255)
      gtgtgTTATAGCTCTCTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGGTCGAAACGCATACCCGGCGTGATGGGCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGGGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg                                                                                                                                                                                                                                      <  2:301516/247‑1 (MQ=255)
      gggCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg                                                                                                                                                                                                                                      <  2:275875/249‑1 (MQ=255)
      gTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAg                                                                                                                                                                                                                                      <  1:470151/251‑1 (MQ=255)
            tataGCTCCCTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATGCACATAGCGATAgcccac                                                                                                                                                                                                                                 <  1:115877/250‑1 (MQ=255)
             ataGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCAGTCACATAGAGTTAgcacacgc                                                                                                                                                                                                                               >  1:83705/1‑251 (MQ=255)
               aGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgccc                                                                                                                                                                                                                             >  1:358920/1‑251 (MQ=255)
               aGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgcc                                                                                                                                                                                                                              <  1:295827/250‑1 (MQ=255)
                 cTCACTTCTTCGTGCGTGCCTGTGCTGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCCACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCCTTTCATGCGTTGATTCATTCACATAGCGATAgcccacgccccc                                                                                                                                                                                                                           >  2:54473/1‑249 (MQ=255)
                 cTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgcccac                                                                                                                                                                                                                           >  1:117027/1‑251 (MQ=255)
                   cACTTCTTCTGGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCTGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTAGCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACt                                                                                                                                                                                                                          <  2:59001/250‑1 (MQ=255)
                    aCTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTc                                                                                                                                                                                                                        >  1:239103/1‑251 (MQ=255)
                         ttcGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGGAGTAACGCAGTTCATGCGTTGATTCATTCACATAGAGATAGCCCACGCCCACTTCGGgcg                                                                                                                                                                                                                   >  1:94263/1‑251 (MQ=255)
                            gtgcgtgcCTGTGATGGCGTTGTTCTGGTCTGACTCGCTATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCCTGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCgg                                                                                                                                                                                                                >  1:400564/1‑251 (MQ=255)
                            gtgcgtgcCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGCTGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGATCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTGGTAACGCCTTTCCTGCGTTGATTCATTCACCTAGCGCTAGCCCACGCCCCCTTCCTGCGCgg                                                                                                                                                                                                                >  2:482232/1‑251 (MQ=255)
                                gtgcCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAgg                                                                                                                                                                                                            <  1:183056/251‑1 (MQ=255)
                                gggcCTGTGATGGCGGTGGTGTTGCATGACTCGATATGTTCGAAACGCATACCCGGCGGGATGGTCCAGTTTCCGATGTCGATTTTGTCATCCAGAGTCCAGGAGTGCGCCTCGGTGCCGGAACGGGTATCGCTGTCGTAAGGGCTTGAGCCGGCCGGCAACTGCCCGCGGCTGGTGGCGGTGTAGTAACGCATTTCATGCGGTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAgg                                                                                                                                                                                                            <  2:2538/249‑1 (MQ=255)
                                      gtgATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCg                                                                                                                                                                                                        <  1:48183/249‑1 (MQ=255)
                                              gttgttCTGGTATGACTCGCTATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGCCGGCAACTGCCCGCTGCTGGTGGCGGTGTAGAAACGCATTTCATGAGTTGATTCTTTCAAATCGCGATAGACCACGCACAATTCGTGCGCGGTCGGGCCGGTCATaaaa                                                                                                                                                                                               >  2:154846/1‑249 (MQ=255)
                                              gttgttCTGGTATGACTCGATATGTTCGAACCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGGGGCGGGGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACCCCAACTTCGTGCGCGGAAGGGCAGATCaaaa                                                                                                                                                                                                 >  2:289515/1‑245 (MQ=255)
                                                tgttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGAt                                                                                                                                                                                             <  1:315839/250‑1 (MQ=255)
                                                 gttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCt                                                                                                                                                                                           <  1:377561/251‑1 (MQ=255)
                                                 gttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCt                                                                                                                                                                                           <  2:117027/251‑1 (MQ=255)
                                                  ttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGACGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGAGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGCTCATAAAGATCtt                                                                                                                                                                                          >  1:453345/1‑250 (MQ=255)
                                                           gACTCGATATGTTCGAAACGCATTCCCGGCGTGATGGTTCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGc                                                                                                                                                                                  <  1:418258/250‑1 (MQ=255)
                                                                   atGTTCGAAACGCATTCGCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCTGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCa                                                                                                                                                                          <  1:477006/250‑1 (MQ=255)
                                                                   atGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCCTGCGGTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGGTCCTAAAGAACTGGCTGGAGCGGGGCTCa                                                                                                                                                                          >  2:183773/1‑250 (MQ=255)
                                                                        cGAAACGCATTCCCGGAGTGATGGTCAAGTGGCCGATGTAGATTTTGTCATCCAGATACCAGGTGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCGTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGAGTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACcg                                                                                                                                                                    <  2:245976/251‑1 (MQ=255)
                                                                           aaCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCgcgc                                                                                                                                                                 >  2:240812/1‑251 (MQ=255)
                                                                                        cGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAc                                                                                                                                                    <  2:218392/251‑1 (MQ=255)
                                                                                         gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGGAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGACCACGCCCACTTCGTGCGAGGCAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTCGTTAcg                                                                                                                                                   >  1:429399/1‑251 (MQ=255)
                                                                                              ggTCCAGGTGCCGATGTCGATTTTGGCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACgcggcg                                                                                                                                              <  2:428386/251‑1 (MQ=255)
                                                                                                  cAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCa                                                                                                                                                                          <  1:34059/219‑1 (MQ=255)
                                                                                                  cAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGAGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGGAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCa                                                                                                                                                                          >  2:34059/1‑219 (MQ=255)
                                                                                                         cGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGGGGCCCAATACCGCGCACCCCGTAGTTACGCGGCGAGAGGGTGATg                                                                                                                                  >  2:76128/1‑251 (MQ=255)
                                                                                                            tGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGCGGGTGATGCGt                                                                                                                               >  1:8273/1‑251 (MQ=255)
                                                                                                               cGATTTTGTCATCCAGATACCAGGCGTGCGTCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCCCTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGtttt                                                                                                                            <  2:158046/251‑2 (MQ=255)
                                                                                                               cGATTTTGTCAGCCAGATACCAGGCGTGCGCCGCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGttt                                                                                                                             <  2:164246/250‑1 (MQ=255)
                                                                                                                gATTTTGTCATCCAGATAACAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCAAGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCCGACGTTACGCGGCGGGGGGTGGATGCGttttc                                                                                                                           >  2:319612/1‑249 (MQ=255)
                                                                                                                                 aCCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCg                                                                                                          >  1:138032/1‑251 (MQ=255)
                                                                                                                                    aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGc                                                                                                                                                                                                          >  1:17340/1‑153 (MQ=255)
                                                                                                                                    aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGc                                                                                                                                                                                                          <  2:17340/153‑1 (MQ=255)
                                                                                                                                    aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTg                                                                                                       >  1:207649/1‑251 (MQ=255)
                                                                                                                                           cgccgcGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAggg                                                                                                <  2:391210/251‑1 (MQ=255)
                                                                                                                                            gccccGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGt                                                                                               <  1:372237/247‑1 (MQ=255)
                                                                                                                                            gcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGt                                                                                               >  1:36060/1‑251 (MQ=255)
                                                                                                                                            gcCTCGGTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCACACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTTGCTTGGCTCAATACCGGGCACCCAGTAGTTACGCGGGGGGAGGGAGATGCGTATGCCTTGCCCCAGGtcgacgatgggccgggg                                                                                               >  2:315808/1‑235 (MQ=255)
                                                                                                                                                  gTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGGAGTAACGCATTTCATGCGGTGATTCATTCACATAGCGATAGCCCACGCCTACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGGAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGgt                                                                                         <  2:476178/251‑1 (MQ=255)
                                                                                                                                                     ccGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTa                                                                                      <  1:240812/251‑1 (MQ=255)
                                                                                                                                                       ggAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGCTCATAAAGATCTGGCTGTAGCGGGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGgtgt                                                                                       >  1:346151/1‑248 (MQ=255)
                                                                                                                                                       ggAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTTGTTACGCGGCGGGAGGGGGATGCGTTTGCCTCGCTCCAGGTAGCCGCTGCGCAGGGTTTgggggt                                                                                       >  2:171324/1‑245 (MQ=255)
                                                                                                                                                        gAACTCGTATCGCGGTCGTAAGGGCTTGCACCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGaa                                                                                   <  1:254569/251‑1 (MQ=255)
                                                                                                                                                        acaCGCGTATAGCGGTCGTACGGGCTTGAGCCGGACGGCACCTGCCCGCTGCTGCTGACGGGGTAGTAAAGCATTTCATGCGTTGATTCATTCACATCGCGATAGCCCACGCCCACTTCGTGCGCGGAACGGCCGAGCATAAAGATCGTGCTGCAGCGTGGGTCAGTACCTCACACCCAGTCGTTTAGCCACGGAGGGGGGAAGCGTTTGCCTCGCTCGAGGTCGCCGCTGCGCCGGGGTTgggggttga                                                                                    >  2:53412/3‑244 (MQ=255)
                                                                                                                                                                aTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGaa                                                                           <  1:319612/251‑1 (MQ=255)
                                                                                                                                                                   gcgGTCGTAAGGGCTTGATCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGAGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTATTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGt                                                                        <  1:46218/251‑1 (MQ=255)
                                                                                                                                                                         gTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGAAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCACTGAATGTTGCAtt                                                                  >  1:3001/1‑251 (MQ=255)
                                                                                                                                                                               gCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCAGGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATgct                                                            <  2:82356/251‑1 (MQ=255)
                                                                                                                                                                                      ccGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCCGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATTCTGGCTGTc                                                     >  2:240996/1‑251 (MQ=255)
                                                                                                                                                                                        ggACGGAAACTGCCCGCTGCTGGGGGCGGTGTAGTAACGCAGTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGGGATGCGTTTGCCTTGCGCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTctg                                                   <  2:453345/251‑1 (MQ=255)
                                                                                                                                                                                          aCGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTctggc                                                 <  1:54473/251‑1 (MQ=255)
                                                                                                                                                                                            ggCAACTGCCCGCTGCTGGGGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTctggctg                                               <  2:400564/251‑1 (MQ=255)
                                                                                                                                                                                                       gctgctGGTGGCGGAGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGGTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGc                                    >  1:8346/1‑251 (MQ=255)
                                                                                                                                                                                                        ctgctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGcc                                   >  1:341890/1‑251 (MQ=255)
                                                                                                                                                                                                             ggtggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAg                                <  1:310360/249‑1 (MQ=255)
                                                                                                                                                                                                              gtggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCCCGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCt                             <  1:334488/251‑1 (MQ=255)
                                                                                                                                                                                                               tggCGGTGTAGTCACGCATTTCCTGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGCGAGGGGGATTCGTTTGCCTTGCTCCAGGTAGCCGCCGCGCAGGGTTTGGGGGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGACGCCCCGGCTc                            >  2:457953/1‑251 (MQ=255)
                                                                                                                                                                                                                          tccctcATTTCCTGCGTTGATTCCTTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGACGCTGCGCAGGGTTTGGGTGTAGAACCCCTGCATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGcg                 >  1:400904/6‑251 (MQ=255)
                                                                                                                                                                                                                           aaCGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGcg                 >  2:390397/1‑250 (MQ=255)
                                                                                                                                                                                                                                aTTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGTCTCGCCAGCTTGCGACGAcc           <  1:404132/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTGATTCATTCACATATCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCg  <  2:239103/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
GAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCG  >  NC_000913/4515102‑4515585

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