Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,293 T→C H41R (CAC→CGC)  fecD ← ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,2930TC97.6% 98.9 / ‑4.2 42H41R (CAC→CGC) fecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (22/19);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

AACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGC  >  NC_000913/4512047‑4512525
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
ccccccGCCATGCCGCCCGAAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTcc                                                                                                                                                                                                                                        <  1:72385/249‑1 (MQ=255)
acccccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTcc                                                                                                                                                                                                                                        <  1:365092/249‑1 (MQ=255)
  ccccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCCCGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGTTTAACACCGAGAATATCCGGTGAGGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCgg                                                                                                                                                                                                                                      <  1:339224/251‑1 (MQ=255)
         aTGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCACGCGTCCGCCCTGcc                                                                                                                                                                                                                                >  1:216397/1‑250 (MQ=255)
             cgcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGAAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCg                                                                                                                                                                                                                           >  1:180508/1‑251 (MQ=255)
             cgcacgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCAGCATGGTTAACCCCGAGAATATCCGGTGATGCCAGCGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGCGGGCTGCACCGCCAAACAGTGCCAGCAGCAAACGCGGCCGTCGTAACTCCATCACTACATAATAATGCTCGCGTCCGGCCTGCCCGTCg                                                                                                                                                                                                                           >  1:91383/5‑251 (MQ=255)
                              agcGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGcc                                                                                                                                                                                                           >  1:378014/1‑250 (MQ=255)
                                    agccccaccACGGGCAGTGACGGCATAAGAAGTAGATCCCCCAGAGAGCCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGAGGCCAGAGGGTTGC‑GCACAATCCCCTTTTTCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcac                                                                                                                                                                                                    <  1:154562/243‑1 (MQ=255)
                                          aTCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCAAAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGACCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGAt                                                                                                                                                                                              >  1:111023/1‑251 (MQ=255)
                                              ggggCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACt                                                                                                                                                                                          <  1:260577/250‑1 (MQ=255)
                                                   aGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGAGCGCGCCACGGCACGGGGATCACTCCCat                                                                                                                                                                                     >  1:304043/1‑251 (MQ=255)
                                                       aCGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATgg                                                                                                                                                                                 >  1:349179/1‑251 (MQ=255)
                                                               agaagTAGAGCCCTCACAGAGGCGAGGCGGGGGGAATGTTGAACACCGAAAAGAGCCGGAGATGCCAGAGGGTTGC‑GCACAATCCCCTGTAGCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATaa                                                                                                                                                                         <  1:212197/251‑1 (MQ=255)
                                                                gaataagagCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATACCCGGTGATGCCAGAGGGTTGC‑GCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                        <  1:83911/246‑1 (MQ=255)
                                                                   gTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGGATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGAGCACTCCCATATGGGGTGATATCAgc                                                                                                                                                                     >  1:118431/1‑251 (MQ=255)
                                                                    aagagCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAgcg                                                                                                                                                                    <  1:66396/250‑1 (MQ=255)
                                                                             cacaGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCCCGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCCCTGGGAGCACTCCCATATGGAGTGATAACAGCGCACCGCCCg                                                                                                                                                           >  2:131069/1‑251 (MQ=255)
                                                                                      cAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCAAGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCa                                                                                                                                                   >  1:226313/1‑250 (MQ=255)
                                                                                                 cATGGTTAACACCGAGAATATCCGGTGATGCCAGATGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGACAGCAGCAAGCGCGCCAGTCGTTCCTCCATCAATACTTAATAATGCACGCGACCGGCCTGCCATTCGGTCAGCATCGCGCGCCACGTCACGGGACTCACTCCCATATGAAGAGATAACCGCGCCCAGCCCGCCAGGTCAAGGGtga‑ttaa                                                                                                                                       >  2:168745/1‑247 (MQ=255)
                                                                                                     gTTAACACCGAGAATATCCGGTGCTGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGGAAACAGCGCACAGCCCGCCAgatcggaaggccggagtgg                                                                                                                                     >  2:207731/1‑231 (MQ=255)
                                                                                                     gTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGcc                                                                                                                                                         <  1:207731/229‑1 (MQ=255)
                                                                                                          cacCGAGAATATCCGGTGATGCCAGAGGGTAGC‑GCCCAATCCCCTGTATCAGCACGCCCGACACGGCGAGGGCTGCACCGACAAACCGTGCCAGCAACAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGACGGACAGCAGCGCGCGCCAAGGCACGGGGATCACTCCCATATGGAGTGATACCATCGCACAGCCCGCCAAGGACAGGGTAA‑TGAAAATAACCAgc                                                                                                                             >  2:238516/1‑251 (MQ=255)
                                                                                                             cGAGAATATCCGGGGAGGCTAGAGG‑TTGCCGCACAAGCCTTGGGATCAGCACCCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCTCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCa                                                                                                                          <  1:352109/251‑1 (MQ=255)
                                                                                                              gagaATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCa                                                                                                                          >  1:366680/1‑250 (MQ=255)
                                                                                                              gagaAGATCCTGTGACGCCAGAGGTGTGC‑GCACAATCCCCTGTAGCAGAACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCTCGGCAGTCGGTACGCCATCAATACATAATAATTCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCAGATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCaa                                                                                                                         <  1:158118/251‑1 (MQ=255)
                                                                                                                    agccGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTc                                                                                                                    <  1:112964/248‑1 (MQ=255)
                                                                                                                       cGGTGCTGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGCGGTCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTAGGTACTCCATAAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCTGCGCGCGCCACGGCACGGGGGTCACTCCCATCTGGAGTGATAACAGCGCACAGCACGCCAGGTCACAGGTGA‑TGAAAATAACCATCGCACTTTTAATcc                                                                                                                >  2:102028/1‑251 (MQ=255)
                                                                                                                                  gagGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCAACGACACACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGAACACTCCCATATGGAGAGATAACAGCGCACAGCCCGCCCGGGCAAGGGCGA‑TGCAAATAACCAGCGCAATTTTCATCC‑TCGcctccccc                                                                                                    >  2:381763/1‑248 (MQ=255)
                                                                                                                                      gTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGGGCCAGCAGCAAGCGCGGCAGTCGGTGCTCCATCCATACATAATAATGCTCGCGTCCGGCCTGCCAGTCAGTCAGCAGCGCACGCCACGGCACAGGGATCACTCCCATATGGAGTGATAACAGCGCACACCCCACAAAGGCAAGGG‑GAGCGAATCTAAACAACGGAAGTTTCATCC‑TCGCCTCCTcacac                                                                                                 >  2:129025/1‑249 (MQ=255)
                                                                                                                                                    cccTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCCCGGCACGGGGATCACTCCCATATGGAGTGATAACAACGCACAGCCCACCAGGGCAAGGGGGA‑TGAGAATAACCAACGCCATTTTAATCC‑TCGCCTCCTcaacacccagaaaaaaca                                                                                    >  2:203459/1‑234 (MQ=255)
                                                                                                                                                      cTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCCGGGCAAGAGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                                                 >  1:183624/1‑251 (MQ=255)
                                                                                                                                                        gtaccAGCACGTCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAggg                                                                                <  1:194358/246‑1 (MQ=255)
                                                                                                                                                          aTCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACCAGCCAGACAAAGCAAGGGCTg                                                                             >  1:238006/1‑251 (MQ=255)
                                                                                                                                                              gCACGCCCGCCACGGCGAGGGCTGCACCGACAAAAAGTGCCAGCAACAAGCGCGGACGTCGGTACTCCCTCAAAACATAATAATGATCGCGTCCGGCCTGCCAGACGGCCCGCAGCACGCACCACGGAACGGGCATCACTCCCATAGGGGGCGATAACAGAGCAACGCACTCCAGAGCACGGGGGA‑GGAAAGTAACCAGCGCACATTTCATCCCTC‑CCTCCTCACAAGCCCAACACAGCAAgtgcaacca                                                                          >  2:283088/1‑242 (MQ=255)
                                                                                                                                                              ccccgacctccACGGCTGGGGCTGCCCCGACAACCAGTGCCAGCAGCAAGCGAGGCAGCCGGTACTCCATCAATCCATAATAATGTTCGCGTCCGGCCTGCCATGCGGTCAGCAGCGCGCGCTACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCaa                                                                         <  2:145662/242‑1 (MQ=255)
                                                                                                                                                                           ggCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCCAGGGtga‑tga                                                                                                                                        >  2:365170/1‑176 (MQ=255)
                                                                                                                                                                           ggCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga‑tga                                                                                                                                        <  1:365170/176‑1 (MQ=255)
                                                                                                                                                                                     gCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAAGGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCTCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCg                                                  <  1:81193/251‑1 (MQ=255)
                                                                                                                                                                                             aaaCAGTGCGAGCAGCAAGCTCGGCAGTCGGTACGCCATCAATACCTAATAATGCTCGCTTCCGTCCTGCCCTTCGGGCAGCAGCGCTCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGTGTGA‑TGAAAATAACCAGCGCAATGTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGTGCTGCCAATCAGCTCCAGCACTGAGCCTGCGGGCAGAtc                                          <  2:375964/251‑1 (MQ=255)
                                                                                                                                                                                             aaaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAGAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtc                                          <  1:203459/251‑1 (MQ=255)
                                                                                                                                                                                                  gTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGCAAATAAACAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGAAAAAGCAAGAGCTGCCAATAAGCCCCAGCACTGCGCCTGCGGCCCGATCTCCgg                                     >  2:14299/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCcagc  <  1:176638/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
AACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGC‑GCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGA‑TGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGC  >  NC_000913/4512047‑4512525

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