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→AGTfecAferric 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

GGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGG  >  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)
                                                                                                                                      agagGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGCTTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTAACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGc                                                                                                         >  1:395906/1‑251 (MQ=255)
                                                                                                                                          ggTGATGCGTTTGCCTTGCTCCAGGGAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTTGCTGGCTGTCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCTTTCATACGGGCGGGTGGATTGCCAGCTAGCGGCGTCGTAATCCTCGCGAGACAGGCGATCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTTTGCACCt                                                                                                     <  2:150717/251‑1 (MQ=255)
                                                                                                                                                             tCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCCCCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCCGGCTGCTGAAGGTGAGCACCTCATCCGGCGCCTATTTGCt                                                                                  >  2:410794/1‑251 (MQ=255)
                                                                                                                                                              ccAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGc                                                                                   <  1:90422/249‑1 (MQ=255)
                                                                                                                                                                 ggTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTc                                                                              >  1:63884/1‑251 (MQ=255)
                                                                                                                                                                        gCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCCTACGGGCGGGTGGATTGCCAGCGATAGGCGTCGTAAAACGCGCGCGACAAGCCGCCGGGCATGAAGGCTTCACCGTCGTACTACTGCAGCGGGCAGGTGAAGGTGGTGACCGCATCCGGGGCatccttgaatttccgcatc                                                                        >  2:232397/1‑233 (MQ=255)
                                                                                                                                                                           gcgcAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTACCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGAGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGgt                                                                    >  1:328776/1‑251 (MQ=255)
                                                                                                                                                                                 ggTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCCCCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGtcgtcg                                                               >  1:439244/1‑250 (MQ=255)
                                                                                                                                                                                  gTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGCCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATGCGGGCGGGTGGATTGCCAGCGATCGTCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGGCGTCGAt                                                             <  2:317765/251‑1 (MQ=255)
                                                                                                                                                                                  gTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCCGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGAt                                                             >  1:371347/1‑251 (MQ=255)
                                                                                                                                                                                    ttGGGTGTCGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAGCTGGTAGGCCAGACCCGCCAGCTTTGGACGACCCCAGAAGCGGGCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATg                                                            <  2:323052/250‑1 (MQ=255)
                                                                                                                                                                                    gtGGGTGTAGAACCCCTGAATGTTGAATTTAGGCGGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTGTTCAGCATCAGGTCGTCGATGc                                                           <  2:341281/250‑1 (MQ=255)
                                                                                                                                                                                     tGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTAGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCCTCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCg                                                          >  1:82722/1‑251 (MQ=255)
                                                                                                                                                                                      gggTGTAGAACCCCTGAATGTTGAATTGATGCTGGCTGTCTGTCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGGCGATGCgg                                                         <  2:383057/251‑1 (MQ=255)
                                                                                                                                                                                           tAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGCCCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGgc                                                    >  2:139543/1‑251 (MQ=255)
                                                                                                                                                                                            aGAACCCCTGACTGTTGAATTTATGCTGGCTGTCTGGATGGAACTGGTAGCCCAGACTCGCCAGCTTGCGCCGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCCCCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGGGgcg                                                   >  2:409722/1‑251 (MQ=255)
                                                                                                                                                                                                   cTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCCGCGCATATTTGCTTTTCCGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCt                                            >  2:317218/1‑251 (MQ=255)
                                                                                                                                                                                                       aTGTTGAATTTATGCTGGCTGTCGGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCGCCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTcgcg                                        <  2:195779/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTCTGGCTGGAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTg                     <  1:238945/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ctggAACTGGTAGCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGGGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGt              >  1:434954/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggttggCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCTAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGACAGTCACTGCCGCGCGTGCCGGAGTAGagcag       <  2:239165/246‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTAGCCCAGATGCGCCAGCTTGCGACGACCCCATAAGCGTTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGTGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGGCGTAATATTGCAGCAGGCTGTTGAAGGTGTGTTCCTCATCCGGGTCATATTTGCTTTTCAGCAGCAGGTCGTCGAGGCGGGTGGCGCTGTGCTCGTGCTATTCACTGCCGCGCGTGTCGGAGTAGagcagc      <  1:396399/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gCCCAGACTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCgg  >  1:123665/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTGTTGAAGGTGTGCACCTCATCCGGCGCATATTTGCTTTTCAGCATCAGGTCGTCGATGCGGGTGGCGCTGTGCTCGCGCCAGTCACTGCCGCGCGTGCCGGAGTAGAGCAGCGCGG  >  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: