Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,515,556 G→A S374S (AGC→AGT)  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,5560GA95.2% 163.9 / ‑5.9 66S374S (AGC→AGT) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (32/28);  minor base C (2/1);  total (35/31)
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.00e+00

GGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGG  >  NC_000913/4515313‑4515800
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
ggTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGcc                                                                                                                                                                                                                                                 >  1:265878/1‑250 (MQ=255)
        aaCGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGAGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCGGGAACTGGTAGCCCAGACTCGCCCGCTTGcg                                                                                                                                                                                                                                         >  2:193376/1‑250 (MQ=255)
          cGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGTCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGcgacg                                                                                                                                                                                                                                      >  2:90352/1‑251 (MQ=255)
            cATTTCATGCGTTGATTCATTCAAATAGCGCTAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTCGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCCGCTTGAgacgac                                                                                                                                                                                                                                    >  1:133036/1‑251 (MQ=255)
            cATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGAAGTTAAGCGGCGAGAGGGTGATGCGTTTGCCTTGCTACAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCAAGCTTGCgacgac                                                                                                                                                                                                                                    >  1:52104/1‑251 (MQ=255)
                  aTGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGAcccc                                                                                                                                                                                                                                 >  1:354913/1‑248 (MQ=255)
                      gTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCg                                                                                                                                                                                                                          >  2:202538/1‑251 (MQ=255)
                              attcaAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTAGCTGGAGCGTGGTTCAATACCGCGCACCCAGTAGTTACGCGGCGAGCGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGAGTGTAGAACCCCTGAATGATGAATTTATGCTGGCTGTCTGGACGGAACTGGTAGCCCAGCCTCGCCAGCGTGCGGCGACTCCAGAAGCgggctaacg                                                                                                                                                                                                                  >  2:233939/1‑244 (MQ=255)
                                tcaAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATAcgg                                                                                                                                                                                                                 <  1:339358/250‑1 (MQ=255)
                                 caAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATAcgggc                                                                                                                                                                                                               >  1:383057/1‑251 (MQ=255)
                                  aaaTAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCATGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGCCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATAcgggcg                                                                                                                                                                                                              >  1:325707/1‑251 (MQ=255)
                                      aGCGATAGCCCACGCCCACTGCGTGCGCGGAAGGGCCGAGCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTAGGCTGGCTGGCTGGCTGGAACTGGGAGCCCAGACTGGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTg                                                                                                                                                                                                          <  2:78946/251‑1 (MQ=255)
                                               tcacgcacacTTCGTGCGCGGAAGGGCCGTTCATCAAGATCTGGATGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGAGCTGCGTTTGCATTGCTCCAGGTAGCCGCTGCGCCGGGTTTGGGTGTAGAACCCCTGCATGTTGAATTTATGCTGGCTGTCGGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACGCGAGAAGCGGTCAGAGGGGCGGGCGGAGTGccgg                                                                                                                                                                                                 >  2:113279/2‑249 (MQ=255)
                                                 acgcccacTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCCATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCCGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCGTGCGACGACCCCAGAAGCGGTCATACGGGCGGGAGGAATGCCAg                                                                                                                                                                                                 >  2:193397/1‑249 (MQ=255)
                                                 acgcccacTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCg                                                                                                                                                                                               <  1:154217/251‑1 (MQ=255)
                                                 acgcccacTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCg                                                                                                                                                                                               <  2:362381/251‑1 (MQ=255)
                                                     ccacTTCGTGCTCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATGTATGCGGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACTGCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCg                                                                                                                                                                                           <  1:380724/251‑1 (MQ=255)
                                                           cGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGAGCGGcgtc                                                                                                                                                                                      <  1:341427/250‑1 (MQ=255)
                                                               cgcgGAAGGGCCTATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCTAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGTCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGAGTGCCAGCGCTCGGCGTCGTAAt                                                                                                                                                                                 <  2:180438/251‑1 (MQ=255)
                                                                gcgGAAGGGCCGATCATAAAGATATGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGGGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGGTTGGGTGTAGAAGCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGCACTGGTAGCCCCGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGTGGGATTGCCAGCGGACGGcgttgaaatc                                                                                                                                                                                >  2:110087/1‑244 (MQ=255)
                                                                         ccGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCTCTCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTTTCTGGCTGGAACTGGGAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTGATACGGGCGGGGGGAATTCCAGCGATCGGCGTCGTAATCCGCGCgaga                                                                                                                                                                       <  1:201020/251‑1 (MQ=255)
                                                                            aTCATAAAGATCTGGCTGTAGCGTGGCTCAATTCCGCGCACCCAGTAGTTACGCGGCGAGAGGGGGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTCGAAACCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGCCCCCAGAAGCGGTCATACGGGCGGGGGGATTGCCAGCGCTAGGCGCCGTAATCCGCGCGAGAc                                                                                                                                                                      >  2:172953/1‑249 (MQ=255)
                                                                                  aggATCGGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCAcc                                                                                                                                                              <  2:373298/249‑1 (MQ=255)
                                                                                  aataTCTGGCTGTAGCGTGGCTCAATACCGCGCTCCCAGTAGTTACGCGGCGGGAGGTGGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCAcc                                                                                                                                                              <  2:16251/248‑1 (MQ=255)
                                                                                      tCTGGCTGTAGCGTGGATCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGAGATGCGTTTGCCTTGCTACAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGGGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCAACGGGc                                                                                                                                                          >  2:437026/1‑251 (MQ=255)
                                                                                        tGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGCCGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAAGCCACCGGGCAt                                                                                                                                                        >  1:111781/1‑251 (MQ=255)
                                                                                             gTAGCGTGGCTCAATCCTGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCGGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTc                                                                                                                                                     <  2:311834/249‑1 (MQ=255)
                                                                                             gTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCg                                                                                                                                                    <  1:354964/250‑1 (MQ=255)
                                                                                                gCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCtt                                                                                                                                                <  2:30767/251‑1 (MQ=255)
                                                                                                gCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCAGCGCGAGACAGGCCACCGGGCATGTCGGCtt                                                                                                                                                >  1:249409/1‑251 (MQ=255)
                                                                                                    ggCTTAATACCGCGCTCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACGCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCGGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGTCATGGCGGCTTCAcc                                                                                                                                            <  2:413068/251‑1 (MQ=255)
                                                                                                    ggCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCAc                                                                                                                                             <  1:437026/250‑1 (MQ=255)
                                                                                                     gCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACcg                                                                                                                                           <  1:103445/251‑1 (MQ=255)
                                                                                                        caaaaCCGCGCACCAACTAGATACGCGGAGAAAGGCCGATGCGTTTGCCTTGCTCCAGGAAGACACCGCGCAGGACTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTAGCCCGCTTGCGCAGACCCCAGAAGCGGCCATACGGGCGGGCGGATTGCCAGCGATCGGCCTCGTAATCAGCACGAGACAGGACACCGGGCATCTCGGCTTCCCcgtcg                                                                                                                                        >  2:391042/5‑251 (MQ=255)
                                                                                                               gcgcACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTCGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATAtt                                                                                                                                 <  1:344643/251‑1 (MQ=255)
                                                                                                                      cAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGcagcag                                                                                                                          <  1:369675/251‑1 (MQ=255)
                                                                                                                        gtaggTACGCGGCGAGAGGGTGATTCGTTTGCCTGTCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTTAATTTATGCTGGCTGTCTGGCTGGACCTGGTAGCCCAGACTCGCCAGCTTGCGCCGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTATTCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGc                                                                                                                        <  2:418330/251‑1 (MQ=255)
                                                                                                                         tacttACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGACGCTGCGCAGGGTTTGGGTGTAGAACCACTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCAGCCGGGCGGGCGTATTGCCAGCGATCGGCGTCGTAATCCCCGCGAGACAGGCCACCGGGCAAGTCGGATGCCACGTCGTAAGATTTCAGCAGGCt                                                                                                                       >  2:238945/4‑251 (MQ=255)
                                                                                                                          agtTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCATCAGGCt                                                                                                                       <  1:115701/250‑1 (MQ=255)
                                                                                                                          agtTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCGGTCTGGCTGGAACTGGGAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTTGATTGCCAGCGATCGTCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCt                                                                                                                       <  1:218647/250‑1 (MQ=255)
                                                                                                                          agtTACGCAGCGAGAGGGTGATGAGTTTGCCTTGCTCCAGGAAGCCGCTGCGCAGGGTTTGGGTGAAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGCCGTAATCCCCGCGCGCCAGGCCACCGGGCCTGTCGGCGTCACCGCCGTAATATTGCAGCAGGc                                                                                                                        >  2:381038/1‑249 (MQ=255)
                                                                                                                                 cggcgAGAGGGTGATGCGTTTACCTTGCTCCAGGGAGCCGCTGCGCAGGGTTTGGGAGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGAACCCAGAAGCGGTCATACGGGCGGGTGGGTTGCCAGCGATCGGCGGCGAAATCCGCGCCAGCCAGGCCACCGGGACTGTCGTCGTCACCGTCGTAAGATTgcgagaggctggtgcagg                                                                                                               >  2:31879/1‑235 (MQ=255)
                                                                                                                                  ggcgAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGACCTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGGAATCAGCGCGAAGCAGGCCCCCGGGCATGGCGGCTTCACCGTCGTAAAATTACCGCAGGCTGTTGagcgg                                                                                                              >  2:164081/1‑247 (MQ=255)
                                                                                                                                      agagGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGCTTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTAACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGc                                                                                                         >  1:395906/1‑251 (MQ=255)
                                                                                                                                          ggTGATGCGTTTGCCTTGCTCCAGGGAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTTGCTGGCTGTCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCTTTCATACGGGCGGGTGGATTGCCAGCTAGCGGCGTCGTAATCCTCGCGAGACAGGCGATCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTT‑TGCACCt                                                                                                     <  2:150717/251‑1 (MQ=255)
                                                                                                                                                             tCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCCCCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCCGGCTGCTGAAGGTG‑AGCACCTCATCCGGCGCCTATTTGCt                                                                                  >  2:410794/1‑251 (MQ=255)
                                                                                                                                                              ccAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGc                                                                                   <  1:90422/249‑1 (MQ=255)
                                                                                                                                                                 ggTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTc                                                                              >  1:63884/1‑251 (MQ=255)
                                                                                                                                                                        gCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCCTACGGGCGGGTGGATTGCCAGCGATAGGCGTCGTAAAACGCGCGCGACAAGCCGCCGGGCATGAAGGCTTCACCGTCGTACTACTGCAGCGGGCAGGTGAAGGTGGTG‑ACCGCATCCGGGGCatccttgaatttccgcatc                                                                        >  2:232397/1‑233 (MQ=255)
                                                                                                                                                                           gcgcAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTACCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGAG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgt                                                                    >  1:328776/1‑251 (MQ=255)
                                                                                                                                                                                 ggTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCCCCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGtcgtcg                                                               >  1:439244/1‑250 (MQ=255)
                                                                                                                                                                                  gTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGCCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATGCGGGCGGGTGGATTGCCAGCGATCGTCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGGCGTCGAt                                                             <  2:317765/251‑1 (MQ=255)
                                                                                                                                                                                  gTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCCGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGAt                                                             >  1:371347/1‑251 (MQ=255)
                                                                                                                                                                                    ttGGGTGTCGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAGCTGGTAGGCCAGACCCGCCAGCTTTGGACGACCCCAGAAGCGGGCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATg                                                            <  2:323052/250‑1 (MQ=255)
                                                                                                                                                                                    gtGGGTGTAGAACCCCTGAATGTTGAATTTAGGCGGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTGTTCAGCATCAGGTCGTCGATGc                                                           <  2:341281/250‑1 (MQ=255)
                                                                                                                                                                                     tGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTAGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCCTCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCg                                                          >  1:82722/1‑251 (MQ=255)
                                                                                                                                                                                      gggTGTAGAACCCCTGAATGTTGAATTGATGCTGGCTGTCTGTCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGGCGATGCgg                                                         <  2:383057/251‑1 (MQ=255)
                                                                                                                                                                                           tAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGCCCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGgc                                                    >  2:139543/1‑251 (MQ=255)
                                                                                                                                                                                            aGAACCCCTGACTGTTGAATTTATGCTGGCTGTCTGGATGGAACTGGTAGCCCAGACTCGCCAGCTTGCGCCGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCCCCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGGGgcg                                                   >  2:409722/1‑251 (MQ=255)
                                                                                                                                                                                                   cTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCCGCGCATATTTGCTTTTCCGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCt                                            >  2:317218/1‑251 (MQ=255)
                                                                                                                                                                                                       aTGTTGAATTTATGCTGGCTGTCGGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCGCCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcgcg                                        <  2:195779/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTg                     <  1:238945/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ctggAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGGG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGt              >  1:434954/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggttggCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCTAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGACAGTCACTGCCGCGCGTGCCGGAGTAGagcag       <  2:239165/246‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTAGCCCAGATGCGCCAGCTTGCGACGACCCCATAAGCGTTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGTGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGGCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGTTCCTCATCCGGGTCATATTTGCTTTTCAGCAGCAGGTCGTCGAGGCGGGTGGCGCTGTGCTCGTGCTATTCACTGCCGCGCGTGTCGGAGTAGagcagc      <  1:396399/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCgg  >  1:123665/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTG‑TGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGG  >  NC_000913/4515313‑4515800

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