Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,652,331 T→C intergenic (+2346/+397) ydfD → / ← ynfP lysis protein/protein YnfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3310TC95.7% 156.5 / NA 46intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/24);  minor base A (2/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGG  >  NC_000913/1652086‑1652568
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
tATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                                                                                                                                                                                               <  1:105548/250‑1 (MQ=255)
 aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                                                                                                                                                                                             <  1:553493/251‑1 (MQ=255)
   ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                                                                                                                                                                                           <  1:36363/251‑1 (MQ=255)
      gtgtgtCTGTGTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATATGAAAAATGGAACAATGATGATCAGGAGAAATGTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGACAGCAATCGACGCTCTCCGGATTCAGATGACATTAACGAGA‑CTGATCAAAGAGCATATCATTGATGTTCACCTCTGAGAGTa                                                                                                                                                                                                                                         <  2:497580/250‑1 (MQ=255)
              gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaa‑c                                                                                                                                                                                                                               <  1:580776/251‑1 (MQ=255)
              gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                                                                                                                                                                                   <  2:174841/248‑1 (MQ=255)
              gTTTAGACTGGAGTTCGCCCTGGAGAGTTGTGTGCACTTGGAGGGGAGGACATAGATCTGAAAAATGGATCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTGTCACGGTACTAAAAACACAGGCAGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATAGCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                                                                                                                                                                                                 <  1:365565/250‑1 (MQ=255)
                    aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGaaa                                                                                                                                                                                                                          <  1:278245/250‑1 (MQ=255)
                       ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTC‑ATCTTATTAAGCCAGCAATCCACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGGGTATGGCAGAA‑CAGGAAAA‑c                                                                                                                                                                                                                      >  1:263431/1‑250 (MQ=255)
                                   ggTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCtt                                                                                                                                                                                                         <  2:128894/251‑1 (MQ=255)
                                           gtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTc                                                                                                                                                                                                <  1:428921/251‑1 (MQ=255)
                                                    tGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCAACGGGGCAATC‑TTCTTATTAAGACAGAAATCGACGCCCTCCGGAGTCAGATGACATTTACGAGA‑CTGAGCAAAGAGCATATCATTGCTGTTTACCTCAAAGAGTATGGCAGAA‑CAGCAAAA‑CAAAAATATACATTT‑GTTTTTCAACCttaa                                                                                                                                                                                        >  2:201900/1‑247 (MQ=255)
                                                     ggATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTC‑ATCTTTTTAAACCAGCAATCGGCGCTCTCCGGAGTCAGATGACATTAACGAAA‑ATGAGCACAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCTGAA‑AAGAAAAA‑AAAAAAAGGCCCTTT‑GTTTTTCAACCTGAAgt                                                                                                                                                                                      >  2:115058/1‑251 (MQ=255)
                                                        tGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTC‑ATTTTATTAAGCCAGCAATCGACGCATTCCGGAGTAAGAAGAAATTAACGAGA‑CTGAGGAAAAAGCATAATATTGGTGTTCAACTCAGAGAGTATGGCAGAA‑GAGGAAAA‑AAAAAAAGGACCTTT‑TTTTTTTAAACTGGAGTGTc                                                                                                                                                                                   >  1:181145/1‑251 (MQ=255)
                                                           taggaCATAGATCTGAAAACTGGAACAATGCTGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAAAGGGTCATTC‑ATCTTATTAAACCAACAAACGACGCTCTACGGGGTCAGATGAAATTTACGGGA‑CTGAGCAAAAAGCATATTATTTTTTTTCACCTCAGTGAGGATGGCA‑AA‑CAAAAAAAACAACAATGCACCTTT‑GTTTTTTAACCTGaaaagtca                                                                                                                                                                                  >  2:332264/2‑243 (MQ=255)
                                                                      tCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGGATGGCAGAA‑AAGAAAAA‑AAAAAATGCACCTTT‑GTTTTTCAACCTGAAGAGTCCGCGAG‑AGGaaaa                                                                                                                                                                     >  1:305916/1‑250 (MQ=255)
                                                                                 gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGaaa                                                                                                                                                                                                                          <  1:469893/189‑1 (MQ=255)
                                                                                 gAACAATGATGAACAGGAGAAATTTAGCAAANGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACAGTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGAc                                                                                                                                                         <  2:277681/251‑1 (MQ=255)
                                                                                 gAACAAAGATGAACAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCCCCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGCCAGCGAG‑AGTAAACAATTATGGTGAc                                                                                                                                                         <  2:314615/251‑1 (MQ=255)
                                                                                          tgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGCAACCAATCGGGGCAATC‑ATCTTATTAAGGCAACAATCGACGCAATCCGGAGTCAGATGACATTAACGAGA‑CTGAGAAAAAAGCATATAATTTATGTTTACATCAGGGAGTATGGCAGAA‑AAGCAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGCG‑AGGAAAAAATTATGGTGACCAttttcac                                                                                                                                                >  1:242776/1‑248 (MQ=255)
                                                                                           gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACATAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCATCCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTAcc                                                                                                                                                <  2:32246/250‑1 (MQ=255)
                                                                                                                 gACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAACCGGGCCAATCTAT‑TTATTAAAACAGAAATAGGCGCCCTCCGGAGTCAGATGAACTTAACGAGAACT‑AGCAAAGCACATATCATTGAAGTTCAACTCCGAGAGGATGGCCGAAACA‑AAAAA‑CAAAAATGCACCTTTTGTTTTTCA‑CCTGAACTGTAAGCGAGGAG‑AAAAAATTATGGTGCACATTTTGCCGGTGCATCAATAAGGCAGAtgtg                                                                                                                         >  2:392646/1‑251 (MQ=255)
                                                                                                                    cGTTTCACGGTACCAAAACCACAGGCGGGAACCAACCGGGTCCTTC‑ATCTTATTAAGCCAGCAATCCACGCCCTCCGGAGTCAGATGACATTAAAGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGGGAGTATGGCAGAA‑CCGAAAAA‑CAAAAAAGCACCTTT‑TTTTTTCAAACAGAAGTGTCAGCGAG‑AGAAAAAAATTATGGGGAACATTTTATCCTTGACTCAATAAAGCAGATGTgggg                                                                                                                      >  2:452756/1‑250 (MQ=255)
                                                                                                                    cGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCCTTT‑ATCTTATTAAACCCACAAACGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑AAAAAAAGCACCCTT‑GTTTTTCAAACTTAAATGTCAGGGAG‑AGAAAAAAATTATGGTGACCCTTTTACCGTTGAATCCATAAGGCAGATGggggg                                                                                                                      >  2:152802/1‑250 (MQ=255)
                                                                                                                    cGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGACATTC‑ATCTTCTTAAACCAGCAAACGAAGCACTCCGGGGTCAGATGACCTTAACAGGA‑CCGAGAAAACAGCAAATCATTGATGTTCACCTCAGAAAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCCCCTTT‑GTTTTTCAACCTGAAATGTCCGCGAG‑AGTAAAAAATTTTTGGGAACCTTTTAACGTTTACTCCATAAGGCAGATGggggg                                                                                                                      >  1:582099/1‑250 (MQ=255)
                                                                                                                     cgtgcACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGAAGCTCTCCTGAGTCAGATGACATTAACCAGA‑CTGAGCAAAGAGGATATCAATTATGGTCACCTCAGAGAGGATGGCAGAA‑CAGAAAAA‑CAAAAATGCCACTTT‑GTTTTTCAACCTGAAGTGTCAGCGAA‑AGTAAAAAATTATTGTGACCCTTTTACCGTTGACTCAATAAGGCAGATGTGGGAc                                                                                                                     <  1:339920/247‑1 (MQ=255)
                                                                                                                                ccAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaa                                                                                                          <  1:332264/251‑1 (MQ=255)
                                                                                                                                  aaaaaCACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGGGGGACGCAGCGATaaaa                                                                                                         >  1:208052/1‑250 (MQ=255)
                                                                                                                                    aatcacaGGCGGGATCCAATTGTGGCATTC‑ATCTTAATAAGCCAGCATTCGGCGCTCTCCGGAGTCAGATGACACGAGCGAGA‑CTGAGCAAAGAGCATATCATTGATTTTCCCCTCAGAGAGTATGGCATAA‑CAGAAAAA‑CAAAAATGCACCTGT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGGAAAAAATTATGCTGACCATGTTACCGTTGACTCGATAGGGCAGATGTGGGACGCAGCGCGAAAACgg                                                                                                      <  2:42473/248‑2 (MQ=255)
                                                                                                                                          aGGCGGGAACCCATCGGGTCATTC‑ATCTTATTAAGCCAGAAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAACAATGCACCTTT‑GTTTTTCAACCTGAAGCGTCAGCGAG‑AGTAAAACATTATGGTGACCCTTTTACCGTTGACTCAATAAGGCAGCTGTGGGACGCAGCGATAAAACGTGCCggc                                                                                                >  2:533626/1‑250 (MQ=255)
                                                                                                                                              gggAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                            >  2:152822/1‑251 (MQ=255)
                                                                                                                                               ggAACCAATCGGGAAGTTC‑ATCTTATTAACACAGCAATCGTCGATATCCGGAGTCAGATGACATTAAAGAGG‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑AAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAACGAG‑AGTAAAAAATTATGGTGACAATTTTACCGTTGACTCAATAAGGCAGAGGTGGGACGCAGCGCTCAAAAGAGCCGGAcaccc                                                                                           >  2:376712/1‑247 (MQ=255)
                                                                                                                                                    cAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                      <  2:224884/251‑1 (MQ=255)
                                                                                                                                                    cAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTACGCCATc                                                                                      >  2:73935/1‑251 (MQ=255)
                                                                                                                                                           gTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                               <  1:152802/251‑1 (MQ=255)
                                                                                                                                                                       tattaAGCCAGCAATCGACGCTGTCGGGAGTCAGATGACATTAACGAGA‑CTGTGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGAATGGCAGAA‑CAGACAAA‑CAAACATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCCGATGTGGTACGCAGAGATAAACCGTGCCGGACTCCGCCATCGAACATAATCTCAGt                                                                       >  2:174257/1‑248 (MQ=255)
                                                                                                                                                                            aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACAt                                                               <  1:15517/251‑1 (MQ=255)
                                                                                                                                                                                         cGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCCTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGGCCATTTTACCGTTGACTCAATAAGTCAGATGTGTGACGCAGCGATAAAACGTGCCGGACTCCGCCAGCGAAAATCATATCAGTCGAGACATACTTATGCctgct                                                  <  2:415910/251‑1 (MQ=255)
                                                                                                                                                                                             ctcCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTc                                              >  1:13617/1‑251 (MQ=255)
                                                                                                                                                                                                             cATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACATCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGCTGTGGGACGCAGCGATAAAACGTGCCGGCCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATTCCTGCTGGTCCCTGACCGCTGGTGCt                              >  2:86625/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCg                        <  1:389025/251‑1 (MQ=255)
                                                                                                                                                                                                                     aga‑CTGAGCAAAGAGCTTAGCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAAGGCACCTTT‑GGTGTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACAATTTGACCGTTGTCTCAATAAAGCAGATGTGGGACGCAGCGATAAAACGTGCCGGTCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGGCGCTGTCAGCTGGTGCTAAGCCGGc                      <  2:444509/251‑1 (MQ=255)
                                                                                                                                                                                                                     aga‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGc                      >  1:265186/1‑251 (MQ=255)
                                                                                                                                                                                                                       a‑CTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCCCCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAt                    <  2:98179/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tCATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCCTCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATg   >  1:522348/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cATTGATGTTCACCTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATgg  >  1:510274/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTC‑ATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGA‑CTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAA‑CAGAAAAA‑CAAAAATGCACCTTT‑GTTTTTCAACCTGAAGTGTCAGCGAG‑AGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGG  >  NC_000913/1652086‑1652568

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