Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,515,972 A→T 5.3% Y236N (TAC→AAC)  fecA ← TonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,9720AT 5.3% 262.0 / 3.6 75Y236N (TAC→AAC) fecATonB‑dependent outer membrane ferric citrate transporter and signal transducer; ferric citrate extracelluar receptor; FecR‑interacting protein
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base T (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCG  >  NC_000913/4515833‑4516118
                                                                                                                                           |                                                                                                                                                  
aTCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcaccgc                                                                                                                                           >  3:410855/1‑149 (MQ=255)
 tCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcaccgca                                                                                                                                          >  3:20755/1‑149 (MQ=255)
 tCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcaccgca                                                                                                                                          >  2:634496/1‑149 (MQ=255)
 tCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAAcgcaccgca                                                                                                                                          >  1:508261/1‑149 (MQ=255)
   aGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAGAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTGCGTTAcggaccgcaga                                                                                                                                        <  3:269274/149‑3 (MQ=255)
               tctTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGGCCTGCGGAATGTCACGGGTAACAAAGTTCACTACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTAc                                                                                                                            >  1:349628/1‑149 (MQ=255)
                 tttCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCa                                                                                                                          >  4:50157/1‑149 (MQ=255)
                 tttCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCa                                                                                                                          >  3:164999/1‑149 (MQ=255)
                  ttCGGGTTGTTGTGTGAAGAGGTTGGCCTGAGCTGACCTTCCACGCCCGCCTCGATCCCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCAc                                                                                                                         <  2:308891/149‑1 (MQ=255)
                      ggTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCa                                                                                                                     >  3:440259/1‑149 (MQ=255)
                        ttgttTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAAt                                                                                                                   >  1:437201/1‑149 (MQ=255)
                                  aaGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATg                                                                                                         >  1:450029/1‑149 (MQ=255)
                                   agagGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGt                                                                                                        <  3:664754/149‑1 (MQ=255)
                                       gTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGcc                                                                                                    <  4:20755/149‑1 (MQ=255)
                                           ggCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGc                                                                                                <  2:218754/149‑1 (MQ=255)
                                           ggCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGc                                                                                                <  2:244484/149‑1 (MQ=255)
                                            gctgagctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCg                                                                                               >  2:253138/1‑149 (MQ=255)
                                                agctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCTGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAAc                                                                                           <  1:624879/149‑1 (MQ=255)
                                                agctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAAc                                                                                           >  1:442024/1‑149 (MQ=255)
                                                agctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAAc                                                                                           <  2:300888/149‑1 (MQ=255)
                                                agctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAAc                                                                                           <  2:430933/149‑1 (MQ=255)
                                                  ctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACgg                                                                                         <  4:440259/149‑1 (MQ=255)
                                                  ctgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACgg                                                                                         <  3:602385/149‑1 (MQ=255)
                                                   tgaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACggg                                                                                        <  3:93495/149‑1 (MQ=255)
                                                    gaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGa                                                                                       <  2:95521/149‑1 (MQ=255)
                                                    gaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGa                                                                                       <  4:164999/149‑1 (MQ=255)
                                                    gaCCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGa                                                                                       <  3:331895/149‑1 (MQ=255)
                                                        ttCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCa                                                                                   <  2:554938/149‑1 (MQ=255)
                                                          ccACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCATCACCGCGTACCACGTCAATGGCATCCAt                                                                                                          >  1:490102/1‑124 (MQ=255)
                                                          ccACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCAt                                                                                                          <  2:490102/124‑1 (MQ=255)
                                                          ccACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGt                                                                                 <  3:50157/149‑1 (MQ=255)
                                                           cACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTg                                                                                <  2:450029/149‑1 (MQ=255)
                                                              gcccgccTCGATACCAAAGTCCTGCGGAATGGTACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGaaa                                                                             <  3:256928/149‑1 (MQ=255)
                                                                ccgccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGGCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGc                                                                           >  1:65377/1‑149 (MQ=255)
                                                                  gccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTg                                                                         >  3:421455/1‑149 (MQ=255)
                                                                     tCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCgg                                                                      >  2:683030/1‑149 (MQ=255)
                                                                             aaaGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGt                                                              >  4:419519/1‑149 (MQ=255)
                                                                                 tCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGa                                                                                                  <  4:2545/109‑1 (MQ=255)
                                                                                 tCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGa                                                                                                  >  3:2545/1‑109 (MQ=255)
                                                                                 tCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAgg                                                          >  2:300630/1‑149 (MQ=255)
                                                                                 tCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAgg                                                          <  1:683030/149‑1 (MQ=255)
                                                                                  ccTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAggg                                                         <  2:437201/149‑1 (MQ=255)
                                                                                   cTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGATACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAgggg                                                        <  1:634496/149‑1 (MQ=255)
                                                                                   cTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGc                                                                        >  4:216003/1‑133 (MQ=255)
                                                                                   cTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGc                                                                        <  3:216003/133‑1 (MQ=255)
                                                                                    tGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGc                                                       >  3:323151/1‑149 (MQ=255)
                                                                                    tGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGc                                                       >  4:496570/1‑149 (MQ=255)
                                                                                          aTGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAggg                                                 >  1:654539/1‑149 (MQ=255)
                                                                                           tGGCACGGGTAACAAAGTTCAACACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACGACCGCGTACCACGTCAATGCCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGCAAGCTGCGGCTGACCGTAAGTGGCGAAgggg                                                <  3:426176/149‑1 (MQ=255)
                                                                                                   gTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACgg                                                                                         <  1:102705/100‑1 (MQ=255)
                                                                                                   gTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACgg                                                                                         >  2:102705/1‑100 (MQ=255)
                                                                                                   gTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGAt                                        <  2:442024/149‑1 (MQ=255)
                                                                                                          aGTTAACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTcc                                 >  2:195627/1‑149 (MQ=255)
                                                                                                            ttCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAgg                                                  <  2:606396/130‑1 (MQ=255)
                                                                                                            ttCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAgg                                                  >  1:606396/1‑130 (MQ=255)
                                                                                                              caccacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCa                             >  3:568521/1‑149 (MQ=255)
                                                                                                               accacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAg                            >  2:628590/1‑149 (MQ=255)
                                                                                                                  acGCCGCCCACGCTCTGCGGTCCCTAACGCACCCCACCACCCCCGCGTACCACGTCAATGGGATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACAGTAAGGGGCGAAGGGGACGGGGATGCCCTCCAACAgggg                         <  2:349628/149‑3 (MQ=255)
                                                                                                                         ccACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAg                  >  4:501901/1‑149 (MQ=255)
                                                                                                                         ccACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAg                  <  1:195627/149‑1 (MQ=255)
                                                                                                                          cACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGc                 <  4:568521/149‑1 (MQ=255)
                                                                                                                          cACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGc                 >  1:558542/1‑149 (MQ=255)
                                                                                                                           aCGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCg                >  2:487571/1‑149 (MQ=255)
                                                                                                                           aCGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCg                <  3:419519/149‑1 (MQ=255)
                                                                                                                            cGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCgg               >  1:332467/1‑149 (MQ=255)
                                                                                                                               tctGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCgg                                                                      >  3:494481/1‑91 (MQ=255)
                                                                                                                               tctGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCgg                                                                      <  4:494481/91‑1 (MQ=255)
                                                                                                                                 tGCGGTGCGTTACGGACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGc          >  1:128658/1‑149 (MQ=255)
                                                                                                                                     gtgcGTTACGGACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGc                                      <  2:660866/114‑1 (MQ=255)
                                                                                                                                     gtgcGTTACGGACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGc                                      >  1:660866/4‑117 (MQ=255)
                                                                                                                                      tCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACggg                                           <  4:406891/111‑1 (MQ=255)
                                                                                                                                      tCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACggg                                           >  3:406891/1‑111 (MQ=255)
                                                                                                                                        cGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGgcgc   <  3:501901/149‑1 (MQ=255)
                                                                                                                                         gTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGcgcg  <  4:410855/149‑1 (MQ=255)
                                                                                                                                         gTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGcgcg  <  4:323151/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
ATCAGGTTGTGCGTCTCTTTCGGGTTGTTTTGTGAAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCG  >  NC_000913/4515833‑4516118

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