Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,511,532 T→C 17.3% E295G (GAG→GGG)  fecD ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,5320TC17.3% 108.4 / 2.9 57E295G (GAG→GGG) fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (8/34);  new base C (3/8);  total (13/44)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGG  >  NC_000913/4511305‑4511752
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
aaaGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCt                                                                                                                                                                                                                                 <  2:564850/228‑1 (MQ=255)
          aaatCAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                                                                                                                                                                                                                      <  2:157459/228‑1 (MQ=255)
          aaacaatttttaccggctCGGTTTCCCGCAACCGTTAGGGCCTGATAGGGCGGTTATCGTCCCCGTTGGCAGTGAGAGTGAAACGGCGTTAAGTTCCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                                                                                                                                                                                                                      <  1:342122/211‑1 (MQ=255)
          aaaaCAGTTTAACGGGGTCGAGTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGGTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                                                                                                                                                                                                                      <  1:156802/228‑1 (MQ=255)
          aaaaCAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                                                                                                                                                                                                                      <  1:400642/228‑1 (MQ=255)
                   tAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGaataa                                                                                                                                                                                                              <  2:234689/227‑1 (MQ=255)
                      cAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAACTACCTTGTCTGTCCCGTAACTGACCGTGAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACAAGACAAACCACGGCGCACCGATAATGTCGGTCCGCCAGACAACCAGGAGCGcccacggggggaggaaaactaa                                                                                                                                                                                                           >  1:543051/1‑208 (MQ=255)
                         cGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGGCCGCCAACCGGGAGCCCCAGTgggggggtgataattttcg                                                                                                                                                                                                        >  2:477796/1‑213 (MQ=255)
                          gTCCATTTCCCGCAACCGTTAGGAACGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGGCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGcc                                                                                                                                                                                                      <  2:6934/228‑1 (MQ=255)
                             gATTTCCCGCAACCTTTAGGACGCGTCCGGGGGGGGTCTGTCACCGTTGGCAGTGAGAGAGAAACGTCGTGAAGTACCTTGGCTGTCCCGTAACTGACCGGCAGATTTTCAGTTCGTAAAGTCATTTAACGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATACTGCTCGc                                                                                                                                                                                                       <  1:250440/224‑1 (MQ=255)
                             gATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCcagc                                                                                                                                                                                                   <  1:67475/228‑1 (MQ=255)
                                     gCAACCGGTGGGGCCGATTAGGGCGGGTATTTTCCCCGTTGGCAGTGAGAGTTAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCgg                                                                                                                                                                                            <  1:391017/227‑1 (MQ=255)
                                              aGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCaa                                                                                                                                                                                 <  2:412858/228‑1 (MQ=255)
                                              aGGACCGATCAGGGCGGTGATCTTCCCCCTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCCACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCaa                                                                                                                                                                                 <  2:81010/228‑1 (MQ=255)
                                              aGGACCCATCAGGGTGGTGAATTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATCCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCaa                                                                                                                                                                                 <  2:418920/228‑1 (MQ=255)
                                                    ccagtcgcgggtggtctttCCCAGGTGGCAGTGAGAGTGAAACGTAGTTAAGTTCCTTGTCTGGCCCGTAACTGACCGGCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCCAACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCCCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCa                                                                                                                                                                           <  1:573476/211‑1 (MQ=255)
                                                      aTCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCCGCACGCCAACCCGGAGCTCCAGTGGGGGGTGAATAATTCTCGCCACCAGAACGGCAACCCCCAAACGCAAc                                                                                                                                                                         >  2:400261/1‑228 (MQ=255)
                                                         aGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTCTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCa                                                                                                                                                                      <  2:57067/228‑1 (MQ=255)
                                                               gTGCTCCTCCCCGGTGGCAGGTAGAAGCAAACCTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTc                                                                                                                                                                <  2:543051/228‑1 (MQ=255)
                                                               gTGATCTTCCCCGTTGGCAGTGAGATTGAAACGTCTTTAAGTACCTTGGCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTc                                                                                                                                                                <  2:228917/228‑1 (MQ=255)
                                                                tGATCTTCCCCGTTGGCCGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTACCTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACAACGGCGCACCGATAATGGCGGGCACCACGCCACACCGTACGCCCAGGTGGGGGTGCAAAATTTTCCCCAGCCGATACGCAACCACCAAccaccactccccactc                                                                                                                                                                 >  1:528331/1‑212 (MQ=255)
                                                                tGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAAACGGGAGCTCCCGTGGGGGGTGAATAATTCTCGCCCGCAGATCGGCAACCACCAACAGCAAGGCAACAGTCa                                                                                                                                                               >  2:356466/1‑228 (MQ=255)
                                                                tGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAAACGGGAGCTCCAGGGGGGGAGGAATAATTCTTGCCACCAGATCGGCAACCACCAACCAGAACGCAcactgca                                                                                                                                                               >  2:117607/1‑222 (MQ=255)
                                                                   tCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGGTAAGTCCCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCCTTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAggg                                                                                                                                                            <  2:389247/228‑1 (MQ=255)
                                                                    ttctCCCCGTTGTCTGTGAGACGGAAACTTCGTTAAGTTCCTTTTCTTTCCCGTAACCGCCCTTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTTCCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAggg                                                                                                                                                            <  1:499958/225‑1 (MQ=255)
                                                                     gtcccccGTGTGCAGTGAGAGTGAAACGTCGTTAAGTTCCTTGGCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAACGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGc                                                                                                                                                           <  2:232255/225‑1 (MQ=255)
                                                                             ttGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCCACCGGGAGCTCCCGTGGGGGATGAATAATTCTCGCCAGCAGACCGGCAACCACCAACAGCAAAGCACCTGGCAGGGCCGaaaacg                                                                                                                                                   >  2:501386/1‑225 (MQ=255)
                                                                               ggCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTAACCAACCCCACCGGGAGCCCCAGTGGGGGGTGAATAATTCTCGCCCGCAGAGCGGCGAACCACAACAGCCAACCACCAGTCAgaggggacaaaagc                                                                                                                                                 >  1:590157/1‑215 (MQ=255)
                                                                                    tGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCCACCCGGAGCTCCAGTAGGGGATGGATAATTCTCGCCCGCCGATCCGCAAACACCCACAGCAACACACCTGTCAGGGCCGAAACAGGacgcag                                                                                                                                            >  2:542208/1‑228 (MQ=255)
                                                                                     gagaGTTAACCGCCGTTAAGTTCCTTTTCTTTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCGCCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGt                                                                                                                                           <  1:626546/228‑1 (MQ=255)
                                                                                                   tttAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTCCAGTTCGAAAAGTCATCTATCGCATTCTCACAAGCAACCAGACAAGCCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTc                                                                                                                              <  2:590157/226‑1 (MQ=255)
                                                                                                    ttAAGTTCCTTGTCTGTCCCGTAACTTACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTcc                                                                                                                             <  2:405002/227‑1 (MQ=255)
                                                                                                        gttctttgttttgcccGTAACTGTCCGTCAGATTTTCAGGTCGTAAAGTCATTTATCGCATTCTCACAAGAAACAAGACAAACCCCGGCGCACCGATAAGGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGAGGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGAAGTCTGCGGTGACGTCCACCgg                                                                                                                        <  1:130701/215‑1 (MQ=255)
                                                                                                        gtgccTTGGCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                                                                                                                        <  2:528331/225‑1 (MQ=255)
                                                                                                        atcccTTTTCTTGCCCGTAACGGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                                                                                                                        <  1:477796/225‑1 (MQ=255)
                                                                                                          ctcttttttgtcCCGTAAATGACCGGTAGATTTTCAGTGCGTAAAGTCAGTTACCGCATTCGCCCAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTg                                                                                                                      <  1:72502/220‑1 (MQ=255)
                                                                                                              tgtctgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCAGGAGCCCCCGTGGGGGGTGAAAAATTCTCGCCCGCAGATCGGCAACCACCAACAGCAACGCACCCGTTAAGGCCCAAAAAGGCAGCAGACCGCGGGGaacacacccgcggagag                                                                                                                  >  1:315909/1‑212 (MQ=255)
                                                                                                                   gtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCCGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGGCCGGGCCGAAACAGGCAACAGTCCGCGGTGACGTCCACCCGGGATGCTACg                                                                                                              >  2:85269/1‑227 (MQ=255)
                                                                                                                    tcCCTTTACTGCCCGTCAGATTTTCAGTTCGTAAAGTTATTTATCGCATTCTCACAAGCAACCAGACAACCCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCATATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGca                                                                                                            <  1:501386/228‑1 (MQ=255)
                                                                                                                                         ggagatgtgtaaaagggacagaTCGCATTCCCAAAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACgaga                                                                                        <  2:141504/206‑1 (MQ=255)
                                                                                                                                                 cGTAAATTCTATTATCGCATTCCCACAAGCAACCCGACAAACCCCGGCGCAACGATAATGGCGGGCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAg                                                                               <  1:184375/228‑1 (MQ=255)
                                                                                                                                                                     tccccaagacaacagaaCAAACCACGGCGCCCCGTAAAGGGCGGTCAGCACGCCAACGGGTGGCGCAAGTGGGGGATGAATAAGTCGCGCCAGCAGATCTGCAACCACCAACAGCAACCCACCTGGCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACTTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAggcggc                                                           <  1:642114/212‑1 (MQ=255)
                                                                                                                                                                              caacaaGACAAACCCCGGCGCCCCGTAAATGGCGGCCAGCCCGCCAACCGGGAGCCCCAGGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCCCCAACAGCAACGCACCTGTCAGGGCCGAAACAGGAGGCAGTCTGCGGGGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGt                                                  <  1:225225/228‑1 (MQ=255)
                                                                                                                                                                                 ccAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCCGGAGCTCCAGTGGGGGATGAAAAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCCGTCAGGGCCGAAACAGGCAGCAGGCTGCGGGGACGTCCACCGGTGATGGCGCGCATCAAATGCGGCACCACGAGACCCATAAAGCAAATCGGGCGCCAGGCGGCCACCCCGGtgaa                                               >  1:151737/1‑225 (MQ=255)
                                                                                                                                                                                            cggtgcgcccGATAAGGGTGGACAGACCGCCACCGGGAACCTCCAGGGGGGCATGAATAATGCTCGCCAGCAGGGCGGCAACCACCAACAGCACCGCACCTGTCAGGGCAGAAACGGGCAGCAGTCTGCGGTGACGTCCACCGGTGAGGCTACGCATCATATGCGGCCCCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:229797/220‑1 (MQ=255)
                                                                                                                                                                                            cgggggcACCGATAATGGCCGACAGCGCGCAAACCGGGAGACCCAGGGGGGGAGGAATATTTCCGGGCAGCAGGTCGTGAAGCCCCAACAGCAACGCGCCTGTCAGGGCCGAAACAGGGAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATAGGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:214328/223‑1 (MQ=255)
                                                                                                                                                                                            ccGGGGCACCGAAAATGGGGGTGAGAGCTCCACACGGAAGCGCCAGGGTGGGAGGAATAATTCCGGCCAGAAGAACGGCAACCACGAACAGCAACGCACCGGTCAGGGCCGAAACAGGCGGCAGTCTGCGGTGACGCCCACCGGTGATGCTACGCATCATAGGCGGCACCACGAGACCAATAAGGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:33533/227‑1 (MQ=255)
                                                                                                                                                                                            ccGGCGCACCGATTAGGGCGGTCCACACCCCAACCGGGAGCTCCAGGGGGGGATGAAATATTCCCGCCAGCAAATCGGCAACCCCCAACCGCAACGCACCTGTCAGGGCCGAAACAGGCAGGAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:400261/227‑1 (MQ=255)
                                                                                                                                                                                                 cccccGATAATGGCGGTCAGCACGCCAACCGGGAGCCCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCt                               <  2:147514/225‑1 (MQ=255)
                                                                                                                                                                                                   aCCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAAAAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGGCTGCGGGGCAGCCCACCCGGGGTGCTGCGCATCAAATGCGGCCCCACCGGCCCAATAAACACAATCGGGGCCCCGGCGGCCCCCCCggcgggaggacagcagcaacaaaa                             >  1:193776/1‑206 (MQ=255)
                                                                                                                                                                                                        tAATGGCGGTCAGCCGGCCACCCGGGAGCCCCAGGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACa                        <  2:326815/228‑1 (MQ=255)
                                                                                                                                                                                                             gcgtccaccACTCCCACCGGGAGCCCCCGTGGGGGATGAATAATTCCCCCCAGCATATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGcc                   <  1:356466/220‑1 (MQ=255)
                                                                                                                                                                                                              cGGTCAGCACGCCAACCGGGATCCCCAGGGGGGAATGAAAAATTCTCGCCGGCAGCTCGGCAACCACCAACAGCACCGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGcc                   <  2:560944/227‑1 (MQ=255)
                                                                                                                                                                                                               ggTCAGCACCCCACCGGGGGCCTCCCGGGGGGGATGAATAATGTCCGCCCGCAGACCGGAAACCCACAACACCAAGCCACGTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCCCCACGAGACCAATAAAGCGAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCa                 <  1:372374/228‑1 (MQ=255)
                                                                                                                                                                                                                   gcgccGCCAACCGGGAGCTCCAGGGGGGGTTGATAAATTCCCGCCGGCAGATCGGCAACGACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCGGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGaa              <  1:244773/223‑1 (MQ=255)
                                                                                                                                                                                                                    gcccggcccagcGGGACCCCCAGTGGCGAGATAATAAGTATCGCCAGCACATGGGCAACCACCGACCGCCACCGCACCGTGAGGGCAGAAACAGGCAGCCGTCTGCGGTGACGTCCACGGGTGATGCTACGCAGCATATGCGGCACCCGGAGACCAATAAAGCTAATCGGGCCGCAGTCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAAt             <  1:450243/217‑1 (MQ=255)
                                                                                                                                                                                                                              ccGGAAGCTCCAGGGGGGGAGGAATAAGTCTCGCCAGCAGATCGGCAACCACCACCAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCCGCGGGGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgg    <  2:652132/227‑1 (MQ=255)
                                                                                                                                                                                                                                  gAGCCCCAGTGGGGGATGAATAATTCTCGCCAGCCGATCGGCAACCACCAACAGCCACGCCCCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgggg  <  2:151737/225‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
AAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGG  >  NC_000913/4511305‑4511752

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