Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,879,935 T→G 24.0% T62P (ACC→CCC)  gyrB ← DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,879,9350TG24.0% 131.2 / 8.1 78T62P (ACC→CCC) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (23/33);  new base G (19/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 8.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATTTTACCCTCGCGCTGGATAACCAGCTCCA‑GTTTTTGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACC  >  NC_000913/3879698‑3880164
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
tatatttCCCCGCGCGGTGGATACCCAG‑TCTAGGTTTGTGGGACAGGGC‑TTTGCTCACCGAACCACCACCGCCGTGCAGACCGCGGGATCTTTTAT‑GGGTGTTAGCGCCAAA‑TTTGCCGCCTGCGTGCAGAACGGTCATGATCACTCCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTg                                                                                                                                                                                                                                       <  2:320698/235‑1 (MQ=255)
           tCGCGCTGGATAACCAGCTCCA‑GTTTTTGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGGGCAGAACGGTCATGATCACTTCCGCCGCCGATCCGCCCCCTTCCGGGAGAATACCGGTCGGACTCCCGCGCCCGACATCC‑TGT‑ACAGAGACAGAGGTATAGAGGAGGATGGGGACGAAAAta                                                                                                                                                                                                                              >  2:309357/1‑236 (MQ=255)
           tCGCGCTGGATAACCAGCTCCA‑GTTTTTGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGGGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCTCCTTCCGGGTGAATACCGGTCGGAATCCCGCGCCCGCCATCC‑TGTAACAG‑GCAAGGGTTATCGGCGTGGAGGGGGAGGAAAAtttt                                                                                                                                                                                                                            >  2:698734/1‑238 (MQ=255)
                             tCCA‑GTTTTTGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGCGGGTAATTTTTTTCCACGGACCCCCGAg                                                                                                                                                                                                          >  1:229065/1‑239 (MQ=255)
                             tCCA‑GTTTTTGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGGGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGCCAGAGTTATCGGCGTGCATGTGGACGAGAATTTTTTTACAGGGACCcgcgc                                                                                                                                                                                                           >  1:1107007/1‑237 (MQ=255)
                                       tGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGGGAAAGGGGACGGTAATTTCTTTACAGTGACCCGCGAGCGCTtcgtc                                                                                                                                                                                                 >  2:125018/1‑239 (MQ=255)
                                          gtcggggCGTTTCTG‑ACCCAACCCCCAACGCGTGGCAGACCGCC‑GACAATTTTTAAGGAGTTATCACAAAA‑TTTACCCACTGGGTGAAGAAGGGTCAGGATCATTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGa                                                                                                                                                                                               <  2:870482/234‑1 (MQ=255)
                                          gACAGGGCATTTACT‑ACCGAAACCCAAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGGTTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGa                                                                                                                                                                                               <  1:703526/238‑1 (MQ=255)
                                             aGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGCAAATGGTGACGATAAATTCTTTCAAGTGAcccacacgccgctcgtcggtatc                                                                                                                                                                                           >  1:844070/1‑219 (MQ=255)
                                                                       cGCCGTGCAGCCCGCGGGACACTTTAT‑AGGAGTTACCGTCAAA‑TTTCCCGCCTGCGTGCAGAACGGTAATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCAtg                                                                                                                                                                 <  1:962813/239‑1 (MQ=255)
                                                                            tGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGCCAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCCTCGGCGTTAGAGGTATCT‑ACAACCCCGAATaaccagtgggg                                                                                                                                                            >  2:213861/1‑229 (MQ=255)
                                                                            tGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGGCGATAATTTCTTTTCAGTGACCCGCCAGCGCTTCGTCGATAGAGGTATTA‑ACCACCCCCAAAACCaattgggg                                                                                                                                                            >  2:168495/1‑232 (MQ=255)
                                                                            tGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGACGGTAATTTCTTTaaa                                                                                                                                                                                                                      >  1:181600/1‑180 (MQ=255)
                                                                                         actcTTTAA‑AGGGTTTACCGTAAAG‑TTGCCCGCCGGCGTGCAGAACGGTCATGATCACTCCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGcc                                                                                                                                                <  1:125018/235‑1 (MQ=255)
                                                                                             ttttaagggAGTTTCCGTAAAT‑TTTCCCGCCTGCGTGCAGAACGGTCATGATCCCTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCa                                                                                                                                            <  1:413180/232‑1 (MQ=255)
                                                                                             tttAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAAAGGTGACGATAATTTCTTTTCAGTGACCCGCGAGCCCTTCGTCGGTGGCGGTTTAT‑ACCACCCCCAATACCACGTGGGGGAGACCGGTCCCGGCa                                                                                                                                            >  2:761159/1‑238 (MQ=255)
                                                                                               tat‑aGGAGTTATCGTCAAA‑TTTACCGCCTGCGGGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCATCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTTTCGGCGTGAATGGGGACGATAATTTATTTACAGTGCCACCCGAGAGCTCCTTCTAGAGCGGTATTT‑AACACCCCGAAAACCCTTTGTTGcggcacgtgtccggaatac                                                                                                                                         >  2:325292/1‑221 (MQ=255)
                                                                                                       gTTTTCGTCATT‑TTTCCCGCTGGCGTGCGGAGCGGTCAGTGTCCCTCCCGCCGCCTGTACGCCCTGTTACGGGTGAACGCCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCg                                                                                                                                   <  2:691464/238‑1 (MQ=255)
                                                                                                        ttATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCTCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCCTCGGAGCTAGCGTTAACT‑AACACCTCGAAAAACATGTGGTGCCGACCGGGGGCGGCAACCGTGACGca                                                                                                                                 >  2:793951/1‑238 (MQ=255)
                                                                                                        ttATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGGTAATTTCTTTTCAGTGACCCCCGAGCGCTTCGGCGATAGCGTTATCT‑ACCCCCCCGAATACCCCTTGGGGCAGAccccgtcctttttaaacacttcc                                                                                                                                 >  1:614825/1‑218 (MQ=255)
                                                                                                        ttATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGCCGGTAATTTCTTTACAGTGCCCCCCGGGCGCTTCGCGGTTAGCGGTATTT‑TACACCCCGAATACCCTTTGGTGAACACCGGAGCCGTCACCCggtttcc                                                                                                                                  >  2:64312/1‑232 (MQ=255)
                                                                                                        ttATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGGGCTTCCTCGATACCGGTATCT‑ACCACCCCCAAAACCACGTGGGGCAAACCggaggcgacaccccgtgcgcg                                                                                                                                 >  2:498538/1‑220 (MQ=255)
                                                                                                         taccGTCAAA‑TTTCCCGGCTGTGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                <  1:1012611/236‑1 (MQ=255)
                                                                                                          ttCGTCAAATTTTCCCGCGTGCTGGCAGAACGGTCAGGATCCCTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATGCCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                <  2:46096/238‑1 (MQ=255)
                                                                                                          aTCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATTTGGTGCAGACCGGGGCCGTCATCCGGGTCGcccg                                                                                                                               >  1:92112/1‑237 (MQ=255)
                                                                                                          aTCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGAAAATTTCTTTACAGTGACCCGCGAGGGCTTCGTCGATAGCGTTATCT‑ACCACCTCCAATACCCTGTGGGGCAGACCCGTGCCgcaacccgggttcccgc                                                                                                                               >  1:1057578/1‑223 (MQ=255)
                                                                                                          aTCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCGCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ATAGTGACAGAGTTATCGGCGTGCATGGTGACGATAATTTCTTTACa                                                                                                                                                                                                                      >  1:472115/1‑153 (MQ=255)
                                                                                                          aTCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCACTCTTCCGGGGGAATAACGGTCGGAATGCCGCGCCCGGCATCC‑TGT‑ACCGAGCCAGCGTTATCGTCGCGAATGTTGACGGCTAATTCTTTACAGTGACTCCCGAGCGCGTCTGCGACAGGCTTATCTCACCACC‑CCAAAACCATGTGGGGAAGAGCGGtagcatcccccaggtccccga                                                                                                                               >  2:741585/1‑219 (MQ=255)
                                                                                                          aTCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCCCCGATACGCCCTCTTCGGGGTGAATACCGGTCAGAATGCCGCCCCCGTCATCC‑TGT‑TCAGTGACAGAGTTATCGCCGTACATGCGGACGGTACTTTCTTTACCGTGACCCCCGCGCGGTGCGTAGCTAGCGTTATCT‑CCCGCCTCGCATACCATCCGGTGCAGAcccgctcccaccccagtgtcgccgc                                                                                                                               >  2:67858/1‑216 (MQ=255)
                                                                                                                         cGCCTGCGGGCAGAACGGGCAAGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCggg                                                                                                                 <  1:213861/239‑1 (MQ=255)
                                                                                                                                           tcctGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCggg                                                                                                                 <  1:391191/218‑1 (MQ=255)
                                                                                                                                                aTCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGCTAGCGTTATCT‑ACCACCCCGAATACCAGTGGGTGCCGACCGGGGCCGCCCTCCGTGTCGCCGCAATACATAACCCGGCGCTTACGCACCGCACCCAGc                                                                                            >  2:339936/1‑237 (MQ=255)
                                                                                                                                                    ctccccccgccgATACCCCCTATTCCGGGTGAATACCTGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTTACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTc                                                                                      <  2:184758/233‑1 (MQ=255)
                                                                                                                                                     ttCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGGTATCT‑ACCACCTCGAATACCATGTGGGGCCGACCGGTGCCGTCCTCCGGGTCCCCCATATAAAAACCCTGGCGCCTAAGcccctgtctcctaaaacaa                                                                                      >  2:391191/1‑223 (MQ=255)
                                                                                                                                                     ttCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGGCGCTAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGAATGCGTTATCT‑TCCACCTCGAAAACCATGTGGGGCAGACCGGGGCCGTCACCCCTCTCCCCCGTTTACATACCCGGGCGCTTCCGCCCCCCAACCAAcccatccg                                                                                     >  2:568964/1‑234 (MQ=255)
                                                                                                                                                     ttCCGCCGCCGATACGCCCTCTTCCGGGGGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGGATGGGTACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATTT‑ACCACCTCGAATACCATGAGGTGCAGGCCCGGGCCGTCCTCCGTGTCGCCGATAAACATACCCCGGCCCtttaaaacaactgccgccgccgacc                                                                                     >  2:902434/1‑216 (MQ=255)
                                                                                                                                                           cgccgATTCGCCCTCTTCCGGGGTAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACtt                                                                               <  2:92112/239‑1 (MQ=255)
                                                                                                                                                             ccgATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGAt                                                                                                                                                                                                                                  <  1:244374/91‑1 (MQ=255)
                                                                                                                                                             ccgATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGACGAt                                                                                                                                                                                                                                  >  2:244374/1‑91 (MQ=255)
                                                                                                                                                                 tACGCCCGCTGCCGGGTGAATACGGGTCGGAATGCCGCGCACGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTCACGATAATTTCTTTACAGTGGCCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATGCCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATa                                                                          <  1:64696/238‑1 (MQ=255)
                                                                                                                                                                     cccTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTTCCGTCATCCGTGTCGCCGATATACATACCCGCGCGCTTTCGCACCCCATCCAGCCCTTTTCGGACTTTCATACTtga                                                                     >  1:479352/1‑239 (MQ=255)
                                                                                                                                                                      cctcctCCGGGGTAATCCCGGTCGGAATGCCGCGACCGTCATCC‑TTG‑ACAGGAGCAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGGCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTggag                                                                    <  2:698661/239‑1 (MQ=255)
                                                                                                                                                                          ttCCGGGTTAATCCCCTTCGGAAGTGCGGGCCCGGACTCC‑TTT‑ACAGGGACAGCGTTATCGGCGTTAATGGTGACGATAATCTCTGTACAGTGCCCCTCGAGGCCTTCGTCGATAGCGTTAGCT‑ACCGCCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGCCTTTGATACTGGAGGAGt                                                                <  2:836142/239‑1 (MQ=255)
                                                                                                                                                                          ttCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGCGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTTGTCGATAGCGTTATCT‑ACCACCTTGAATACCATGAGGTGCAGACCGGGGCCGGCCTCCGTTTCGCCGTCATAACTACCCTGGCGCTTACGAACCCCATCCACCCCttttcagcctttggtaacgggggggc                                                                >  2:233442/1‑216 (MQ=255)
                                                                                                                                                                          gccaggttGAATTCCGGTCGGAATGCCCCCGCCGTCAACC‑TGT‑ACAGAGACAGAGTTATCGGCTGGAATGGTGACGATAATTTCTTTACAGGGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTggagga                                                                  <  2:339940/230‑1 (MQ=255)
                                                                                                                                                                             cGGGTTAATCCCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGGTATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCAt                                                             <  1:568964/239‑1 (MQ=255)
                                                                                                                                                                                gTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGat                                                                                                                              <  1:248296/171‑1 (MQ=255)
                                                                                                                                                                                    aTACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGCGACAGCGTTATCGGCGTGAACGGGGACGATAATTTCTTTACAAGGACCCCCCAGCGCTTCGTCGATATCGTTATCT‑ACCACCTCGTATACCATGTGGTGCCGAGCGGGGACGCCTTCCGTGTCGCCGATTTACATACCCGGTCGCCTACTGCCCGCCCCCACCCCTTTCACGAACTTTAAACTCGAGGAGtacaaaaaaat                                                      >  2:370449/1‑229 (MQ=255)
                                                                                                                                                                                    aTACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGGATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCCTCCATAACGTTATCT‑CCCACCTCGAGAAACATGTGGTGCAGACCGGCGCCGTCATCCGGTGCGCCGACATCCAAATCCTGGCGCTTCGCCCCCGCACCCATCCCTTTCGCAACTTTGTAACTGGAGGGGCCACAAGaaat                                                      >  2:266624/1‑237 (MQ=255)
                                                                                                                                                                                    aTACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAAtt                                                      <  1:168495/239‑1 (MQ=255)
                                                                                                                                                                                     tACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTc                                                     <  1:579498/239‑1 (MQ=255)
                                                                                                                                                                                         gggtggAATGCCGCGCCCGTCATCCTTGT‑A‑AGGGAACGAGTTATCGGCGTGAATGGGTACGATAATTTCCTTCCAGTGCCCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACa                                                 <  2:714792/235‑1 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATa                                                                                                                                                                      <  1:168085/121‑1 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCCCATCCAGCCCTTTCCGGACTTTGATCCTGGAGGAGTCCAAAGAATTCGACAt                                                >  1:412792/1‑239 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCCGGCGCTTACGCACCGCATCCCGCCCTTTCAGGACTTTGATATTGGCGGAGCCATAAGaa                                                        >  1:339940/1‑231 (MQ=255)
                                                                                                                                                                                          gTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATa                                                                                                                                                                      >  2:168085/1‑121 (MQ=255)
                                                                                                                                                                                              gAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGGATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCCTTTGGGGCAGACCGGTGCCCCCCTCCCTGTCGCCGCGATACATACCCCGGCGCTTACGCACCCCCTCCCCCCCCTTCCCGACCTTTATACTTGGGGGGTCAAAAGAAtttaaacacaaa                                            >  2:1011147/1‑229 (MQ=255)
                                                                                                                                                                                               aaTGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATaa                                                                                                                                                                                                                                <  2:213618/59‑1 (MQ=255)
                                                                                                                                                                                               aaTGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATaa                                                                                                                                                                                                                                >  1:213618/1‑59 (MQ=255)
                                                                                                                                                                                               aaTGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGGGGAGTCCTAAGAATTCGACATCaaag                                           >  1:714792/1‑237 (MQ=255)
                                                                                                                                                                                                     gcgcCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGGGGCCTAAGAATTCGACATCAACGttttc                                      >  1:908296/1‑236 (MQ=255)
                                                                                                                                                                                                      ccccctttaccc‑TGT‑ACAGAGACAGAGTTTTCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTATGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCg                                    <  2:412792/229‑1 (MQ=255)
                                                                                                                                                                                                             cATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCAtt                              <  1:20228/238‑1 (MQ=255)
                                                                                                                                                                                                                 c‑TGT‑ACAGAGACAGAGTTATCGGTGTGAATGGTGACGATAATTTCTTTACCGTGACCCGCGGGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATAc                         <  1:309357/239‑1 (MQ=255)
                                                                                                                                                                                                                  cttt‑ACAGAGACAGAGTTATCGGTGTTAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCC‑ACCCCCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATAc                         <  1:698734/236‑1 (MQ=255)
                                                                                                                                                                                                                   tGT‑ACAGAGACAGAGTTACCGGCGTTAATGGGGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACt                        <  1:726933/239‑1 (MQ=255)
                                                                                                                                                                                                                     t‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATaa                                                                                                                                                                                                                                <  1:269645/38‑1 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                >  1:500456/1‑128 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCg                                                                                                                                <  2:500456/128‑1 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTAAGCACCGCATCCAGCCCTTTCAGGACCTTGATACTGGAGGAGGCATAAGAATTCGACATCAACGTTTTTTGCTCATTTTTACTTTGGTTa                 >  1:411566/1‑239 (MQ=255)
                                                                                                                                                                                                                           agacagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCCTTATACAAACCCGGGCGCTTACGACCCGCCTCCACCCCTTTCAGGGCTTTGTTACTGGGGGAGGCATAAGAATTCGACCTCAACCTTTCTCCCTCATTTTTCTCTGGGTTa                 >  2:521872/1‑239 (MQ=255)
                                                                                                                                                                                                                             acagagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTTCGCACCGCATCCAGCCCTTTCAGGACTTTGATACTTGAGGAGGCATAAGAATTCGACATCAACGATTCCCGCCCCTTTATACTTGGGTTaa                >  1:535529/1‑238 (MQ=255)
                                                                                                                                                                                                                                gagTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCCGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTTTCGCTCATTTATACTTGGGTTAATCCg            >  1:999033/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGGCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCg            <  1:498538/237‑1 (MQ=255)
                                                                                                                                                                                                                                     gtCGGCGTTAATGGTTACGATAATTTCTTTACAGTGGCCCGCGAGCGGTTCGTCGATAGCGTTATCT‑ACCGCCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTAtt       <  1:691132/238‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGTTACTGGAGGAGGCATAAGAATTCGACATCACCGTTTCCCGCCCATTTTTACTTGGGTTAATCCGTTAtttttc   >  1:586584/1‑235 (MQ=255)
                                                                                                                                                                                                                                           gTGAATGGTGACGATAATTTCTTTACAGGGCCCCGCGAGCGCTTCGGCGATAGCGTTATCT‑ACCCCCTCGAATACCATGTGGTGCAGACCGGGGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTAcc  <  2:999033/238‑1 (MQ=255)
                                                                                                                                                                                                                                           gTGAATGGTGACGATAATTTCTTTACAGGGACCCGCGAGCGCTTCGTCGATCGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGTTTCCGTCATCCGTGTCCCCGATATACATCCCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACATTGTAACTGGGGGTTTATAAAGAATTCGACAGCCCCCTTTTTCTCCCATTTTACCTTGTGGTAATACGTttttttcc   >  2:691132/1‑235 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                      
GAATTTTACCCTCGCGCTGGATAACCAGCTCCA‑GTTTTTGCGACAGGGCGTTTACT‑ACCGAAACACCAACGCCGTGCAGACCGCCGGACACTTTAT‑AGGAGTTATCGTCAAA‑TTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCC‑TGT‑ACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCT‑ACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACC  >  NC_000913/3879698‑3880164

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