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

AAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGG  >  NC_000913/4511305‑4511752
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
aaaGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCt                                                                                                                                                                                                                                 <  2:564850/228‑1 (MQ=255)
          aaatCAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTgg‑ggg                                                                                                                                                                                                                      <  2:157459/228‑1 (MQ=255)
          aaacaatttttaccggctCGGTTTCCCGCAACCGTTAGGGCCTGAT‑AGGGCGGTTATCG‑TCCCCGTTGGCAGTGAGAGTGAAACGGCGTTAAGTTCCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTgg‑ggg                                                                                                                                                                                                                      <  1:342122/211‑1 (MQ=255)
          aaaaCAGTTTAACGGGGTCGAGTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGGTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTgg‑ggg                                                                                                                                                                                                                      <  1:156802/228‑1 (MQ=255)
          aaaaCAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTgg‑ggg                                                                                                                                                                                                                      <  1:400642/228‑1 (MQ=255)
                   tAACAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGaataa                                                                                                                                                                                                              <  2:234689/227‑1 (MQ=255)
                      cAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAACTACCTTGTCTGTCCCGTAACTGACCGTGAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACAAGACAAACCACGGCGCACCGATAATGTCGGTCCGCCAGACA‑ACCAGGAGCGcccacggggggaggaaaactaa                                                                                                                                                                                                           >  1:543051/1‑208 (MQ=255)
                         cGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGGCCGCCA‑ACCGGGAGCCCCAGTgg‑gggggtgataattttcg                                                                                                                                                                                                        >  2:477796/1‑213 (MQ=255)
                          gTCCATTTCCCGCAACCGTTAGGAAC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGGCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGcc                                                                                                                                                                                                      <  2:6934/228‑1 (MQ=255)
                             gATTTCCCGCAACCTTTAGGACG‑CGTCCGGG‑GGGGGTCTGTCACCGTTGGCAGTGAGAGAGAAACGTCGTGAAGTACCTTGGCTGTCCCGTAACTGACCGGCAGATTTTCAGTTCGTAAAGTCATTTAACGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATACTGCTCGc                                                                                                                                                                                                       <  1:250440/224‑1 (MQ=255)
                             gATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCcagc                                                                                                                                                                                                   <  1:67475/228‑1 (MQ=255)
                                     gCAACCGGTGGGGCC‑GATTAGGGCGGGTATTT‑TCCCCGTTGGCAGTGAGAGTTAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCgg                                                                                                                                                                                            <  1:391017/227‑1 (MQ=255)
                                              aGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCaa                                                                                                                                                                                 <  2:412858/228‑1 (MQ=255)
                                              aGGACC‑GATCAGGGCGGTGATCT‑TCCCCCTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCCACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCaa                                                                                                                                                                                 <  2:81010/228‑1 (MQ=255)
                                              aGGACC‑CATCAGGGTGGTGAATT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATCCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCaa                                                                                                                                                                                 <  2:418920/228‑1 (MQ=255)
                                                    ccagtcgcgggtggtctt‑tCCCAGGTGGCAGTGAGAGTGAAACGTAGTTAAGTTCCTTGTCTGGCCCGTAACTGACCGGCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCCAACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCCCGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑Ca                                                                                                                                                                           <  1:573476/211‑1 (MQ=255)
                                                      aTCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCCGCACGCCA‑ACCCGGAGCTCCAGTGG‑GGGGTGAATAATTCTCGCCACCAGAACGGC‑AACCCCCAAACG‑CAAc                                                                                                                                                                         >  2:400261/1‑228 (MQ=255)
                                                         aGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTCTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCa                                                                                                                                                                      <  2:57067/228‑1 (MQ=255)
                                                               gTGCTCC‑TCCCCGGTGGCAGGTAGAAGCAAACCTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTc                                                                                                                                                                <  2:543051/228‑1 (MQ=255)
                                                               gTGATCT‑TCCCCGTTGGCAGTGAGATTGAAACGTCTTTAAGTACCTTGGCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTc                                                                                                                                                                <  2:228917/228‑1 (MQ=255)
                                                                tGATCT‑TCCCCGTTGGCCGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTACCTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACAACGGCGCACCGATAATGGCGGGCACCACGCCACACCGTACGC‑CCAGGTG‑GGGGTGCAAAATTTTCCCCAGCCGATACGC‑AACCACCAAccaccactccccactc                                                                                                                                                                 >  1:528331/1‑212 (MQ=255)
                                                                tGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑AACGGGAGCTCCCGTGG‑GGGGTGAATAATTCTCGCCCGCAGATCGGC‑AACCACCAACAG‑CAAGGCAACAGTCa                                                                                                                                                               >  2:356466/1‑228 (MQ=255)
                                                                tGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑AACGGGAGCTCCAGGGG‑GGGAGGAATAATTCTTGCCACCAGATCGGC‑AACCACCAACCA‑GAACGCAcactgca                                                                                                                                                               >  2:117607/1‑222 (MQ=255)
                                                                   tCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGGTAAGTCCCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCCTTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAggg                                                                                                                                                            <  2:389247/228‑1 (MQ=255)
                                                                    ttctCCCCGTTGTCTGTGAGACGGAAACTTCGTTAAGTTCCTTTTCTTTCCCGTAACCGCCCTTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTTCCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAggg                                                                                                                                                            <  1:499958/225‑1 (MQ=255)
                                                                     gtcccccGTGTGCAGTGAGAGTGAAACGTCGTTAAGTTCCTTGGCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAACGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGc                                                                                                                                                           <  2:232255/225‑1 (MQ=255)
                                                                             ttGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCC‑ACCGGGAGCTCCCGTGG‑GGGATGAATAATTCTCGCCAGCAGACCGGC‑AACCACCAACAG‑CAAAGCACCTGGCAGGGCCGaaaacg                                                                                                                                                   >  2:501386/1‑225 (MQ=255)
                                                                               ggCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTAACCAACCCC‑ACCGGGAGCCCCAGTGG‑GGGGTGAATAATTCTCGCCCGCAGAGCGGCGAACCA‑CAACAG‑CCAACCACCAGTCAgaggggacaaaagc                                                                                                                                                 >  1:590157/1‑215 (MQ=255)
                                                                                    tGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCC‑ACCCGGAGCTCCAGTAG‑GGGATGGATAATTCTCGCCCGCCGATCCGC‑AAACACCCACAG‑CAACACACCTGTCAGGGCCGAAACAGGacgcag                                                                                                                                            >  2:542208/1‑228 (MQ=255)
                                                                                     gagaGTTAACCGCCGTTAAGTTCCTTTTCTTTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCGCCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGt                                                                                                                                           <  1:626546/228‑1 (MQ=255)
                                                                                                   tttAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTCCAGTTCGAAAAGTCATCTATCGCATTCTCACAAGCAACCAGACAAGCCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTc                                                                                                                              <  2:590157/226‑1 (MQ=255)
                                                                                                    ttAAGTTCCTTGTCTGTCCCGTAACTTACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTcc                                                                                                                             <  2:405002/227‑1 (MQ=255)
                                                                                                        gttctttgttttgcccGTAACTGTCCGTCAGATTTTCAGGTCGTAAAGTCATTTATCGCATTCTCACAAGAAACAAGACAAACCCCGGCGCACCGATAAGGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGAGGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGAAGTCTGCGGTGACGTCCACCgg                                                                                                                        <  1:130701/215‑1 (MQ=255)
                                                                                                        gtgccTTGGCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAGCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                                                                                                                        <  2:528331/225‑1 (MQ=255)
                                                                                                        atcccTTTTCTTGCCCGTAACGGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                                                                                                                        <  1:477796/225‑1 (MQ=255)
                                                                                                          ctcttttttgtcCCGTAAATGACCGGTAGATTTTCAGTGCGTAAAGTCAGTTACCGCATTCGCCCAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTg                                                                                                                      <  1:72502/220‑1 (MQ=255)
                                                                                                              tgtctgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCAGGAGCCCCCGTGG‑GGGGTGAAAAATTCTCGCCCGCAGATCGGC‑AACCACCAACAG‑CAACGCACCCGTTAAGGCCCAAAAAGGCAGCAGACCGCGGGGaacacacccgcggagag                                                                                                                  >  1:315909/1‑212 (MQ=255)
                                                                                                                   gtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCCGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGGCCGGGCCGAAACAGGCAACAGTCCGCGGTGACGTCCACCCGGGATGCTACg                                                                                                              >  2:85269/1‑227 (MQ=255)
                                                                                                                    tcCCTTTACTGCCCGTCAGATTTTCAGTTCGTAAAGTTATTTATCGCATTCTCACAAGCAACCAGACAACCCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCATATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGca                                                                                                            <  1:501386/228‑1 (MQ=255)
                                                                                                                                         ggagatgtgtaaaagggacagaTCGCATTCCCAAAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACgaga                                                                                        <  2:141504/206‑1 (MQ=255)
                                                                                                                                                 cGTAAATTCTATTATCGCATTCCCACAAGCAACCCGACAAACCCCGGCGCAACGATAATGGCGGGCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAg                                                                               <  1:184375/228‑1 (MQ=255)
                                                                                                                                                                     tccccaagacaacagaaCAAACCACGGCGCCCCGTAAAGGGCGGTCAGCACGCCA‑ACGGGTGGCGCAAGTGG‑GGGATGAATAAGTCGCGCCAGCAGATCTGC‑AACCACCAACAG‑CAACCCACCTGGCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACTTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAggcggc                                                           <  1:642114/212‑1 (MQ=255)
                                                                                                                                                                              caacaaGACAAACCCCGGCGCCCCGTAAATGGCGGCCAGCCCGCCA‑ACCGGGAGCCCCAGGGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCCCCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGAGGCAGTCTGCGGGGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGt                                                  <  1:225225/228‑1 (MQ=255)
                                                                                                                                                                                 ccAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCCGGAGCTCCAGTGG‑GGGATGAAAAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCCGTCAGGGCCGAAACAGGCAGCAGGCTGCGGGGACGTCCACCGGTGATGGCGCGCATCAAATGCGGCACCACGAGACCCATAAAGCAAATCGGGCGCCAGGCGGCCACCCCGGtgaa                                               >  1:151737/1‑225 (MQ=255)
                                                                                                                                                                                            cggtgcgcccGATAAGGGTGGACAGACCGCCA‑CCGGGAACCTCCAGGGG‑GGCATGAATAATGCTCGCCAGCAGGGCGGC‑AACCACCAACAG‑CACCGCACCTGTCAGGGCAGAAACGGGCAGCAGTCTGCGGTGACGTCCACCGGTGAGGCTACGCATCATATGCGGCCCCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:229797/220‑1 (MQ=255)
                                                                                                                                                                                            cgggggcACCGATAATGGCCGACAGCGCGCAA‑ACCGGGAGACCCAGGGG‑GGGAGGAATATTTCCGGGCAGCAGGTCGTG‑AAGCCCCAACAG‑CAACGCGCCTGTCAGGGCCGAAACAGGGAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATAGGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:214328/223‑1 (MQ=255)
                                                                                                                                                                                            ccGGGGCACCGAAAATGGGGGTGAGAGCTCCA‑CACGGAAGCGCCAGGGT‑GGGAGGAATAATTCCGGCCAGAAGAACGGC‑AACCACGAACAG‑CAACGCACCGGTCAGGGCCGAAACAGGCGGCAGTCTGCGGTGACGCCCACCGGTGATGCTACGCATCATAGGCGGCACCACGAGACCAATAAGGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:33533/227‑1 (MQ=255)
                                                                                                                                                                                            ccGGCGCACCGATTAGGGCGGTCCACACCCCA‑ACCGGGAGCTCCAGGGG‑GGGATGAAATATTCCCGCCAGCAAATCGGC‑AACCCCCAACCG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGGAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGa                                    <  1:400261/227‑1 (MQ=255)
                                                                                                                                                                                                 cccccGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCCCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCt                               <  2:147514/225‑1 (MQ=255)
                                                                                                                                                                                                   aCCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAAAAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGGCTGCGGGGCAGCCCACCCGGGGTGCTGCGCATCAAATGCGGCCCCACCGGCCCAATAAACACAATCGGGGCCCCGGCGGCCCCCCCggcgggaggacagcagcaacaaaa                             >  1:193776/1‑206 (MQ=255)
                                                                                                                                                                                                        tAATGGCGGTCAGCCGGCCA‑CCCGGGAGCCCCAGGGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACa                        <  2:326815/228‑1 (MQ=255)
                                                                                                                                                                                                             gcgtccaccACTCCC‑ACCGGGAGCCCCCGTGG‑GGGATGAATAATTCCCCCCAGCATATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGcc                   <  1:356466/220‑1 (MQ=255)
                                                                                                                                                                                                              cGGTCAGCACGCCA‑ACCGGGATCCCCAGGGG‑GGAATGAAAAATTCTCGCCGGCAGCTCGGC‑AACCACCAACAG‑CACCGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGcc                   <  2:560944/227‑1 (MQ=255)
                                                                                                                                                                                                               ggTCAGCACCCCA‑CCGGGGGCCTCCCGGGG‑GGGATGAATAATGTCCGCCCGCAGACCGGA‑AACCCACAACAC‑CAAGCCACGTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCCCCACGAGACCAATAAAGCGAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCa                 <  1:372374/228‑1 (MQ=255)
                                                                                                                                                                                                                   gcgccGCCA‑ACCGGGAGCTCCAGGGG‑GGGTTGATAAATTCCCGCCGGCAGATCGGC‑AACGACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCGGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGaa              <  1:244773/223‑1 (MQ=255)
                                                                                                                                                                                                                    gcccggcccagcGGGACCCCCAGTGGCGAGAT‑AATAAGTATCGCCAGCACATGGGC‑AACCACCGACCGCCACCGCACC‑GTGAGGGCAGAAACAGGCAGCCGTCTGCGGTGACGTCCACGGGTGATGCTACGCAGCATATGCGGCACCCGGAGACCAATAAAGCTAATCGGGCCGCAGTCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAAt             <  1:450243/217‑1 (MQ=255)
                                                                                                                                                                                                                              ccGGAAGCTCCAGGGG‑GGGAGGAATAAGTCTCGCCAGCAGATCGGC‑AACCACCACCAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCCGCGGGGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgg    <  2:652132/227‑1 (MQ=255)
                                                                                                                                                                                                                                  gAGCCCCAGTGG‑GGGATGAATAATTCTCGCCAGCCGATCGGC‑AACCACCAACAG‑CCACGCCCCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgggg  <  2:151737/225‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                               
AAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACC‑GATCAGGGCGGTGATCT‑TCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCA‑ACCGGGAGCTCCAGTGG‑GGGATGAATAATTCTCGCCAGCAGATCGGC‑AACCACCAACAG‑CAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGG  >  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: ‑