Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,152 C→T T261I (ACC→ATC)  flgJ → putative peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,1520CT92.3% 111.5 / NA 40T261I (ACC→ATC) flgJputative peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (19/17);  minor base G (0/3);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATA  >  NC_000913/1137922‑1138396
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ggggCAACGGCAAAGCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTATGGTGTCAAAGCCTCTGGCAACTGGACAGGGCCAGTTACTGAAATCACCACGAATGACTATGAAAACGGCGAAGCGCAGAAAGTAAAAGCGTAGTTTCGCGTCTACAGCGCGTTTCTGGAACACTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCTCGCGACGCCACCGTGACGTTCGACGCGCGTGCGGAACAgg                                                                                                                                                                                                                                  >  2:80189/1‑251 (MQ=255)
                                  ccgagcTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCAc                                                                                                                                                                                                >  1:18272/1‑251 (MQ=255)
                                     agcTATAACCTGTTTGGTGTCAAAGCCTCTGGCACCTGGAAAGGGCCAGTAACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGCAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGATCGTATCTGGAAGCCTTGTCGGCTTAAGTTGGGCTGTTAACGCGTAACCCGCGCTACGCAGCCGGGAAGACCGACGCCAATGAGGAAAAGAGGGCGAAGGCCATACAGGCCGCGGGCGAGGCCAcc                                                                                                                                                                                               >  2:272966/1‑249 (MQ=255)
                                         ataACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATc                                                                                                                                                                                           <  1:429901/249‑1 (MQ=255)
                                           accctgtgtggggcgaaaGCCTCTGGCAACGGGAACGGGCCAGTTACTGAAATCACCTCGACGGAATATTAAAACGGCGAAGCGAAGAAAGTAAAAGCTAATTTTCTCGTCGTCAGCTCGTATCTGGAAGCCTTGTCGTATTCCGTTGGGCTGTTAACGCGTACCCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCTGGCTATGCCACCGATCCTCa                                                                                                                                                                                       <  2:499634/236‑1 (MQ=255)
                                                          aaaGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACt                                                                                                                                                                        >  1:271803/1‑251 (MQ=255)
                                                                          ggAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGCCGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACAAGATTcac                                                                                                                                                        >  2:408216/1‑250 (MQ=255)
                                                                              aGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCCAACTCACCAACCAGATTcagc                                                                                                                                                       >  2:333340/1‑248 (MQ=255)
                                                                               gggCCAGTTACTGAAATAACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGGTTCGCTTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGAt                                                                                                                                                   <  2:447601/251‑1 (MQ=255)
                                                                                 gCCAGTTACGGAAATCCCCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGa                                                                                                                                                 <  2:494585/251‑1 (MQ=255)
                                                                                          tGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATa                                                                                                                                         >  1:161799/1‑250 (MQ=255)
                                                                                                caccacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCa                                                                                                                                                                                                 >  1:366210/1‑188 (MQ=255)
                                                                                                caccacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCGGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCa                                                                                                                                                                                                 <  2:366210/188‑1 (MQ=255)
                                                                                                 accacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGCTTCAGCAGATGAAATCGATAAGCGACa                                                                                                                                  >  2:96545/1‑250 (MQ=255)
                                                                                                      gACTGAATATGAAAACGGCGAAGCTAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCGGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACACGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGa                                                                                                                            <  1:139685/251‑1 (MQ=255)
                                                                                                          gAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCaa                                                                                                                        <  1:202629/251‑1 (MQ=255)
                                                                                                              atGAAAACGGCGAAGCGAGTAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTCCGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAAcc                                                                                                                    <  1:394019/251‑1 (MQ=255)
                                                                                                               tGAAAACGGCGACGCGACGAACGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCATTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGTCGCGGGCTATGCCACCGAACCTCACTATGCCCGCAAACCCACCAACATGATTCAGCAGATGAAATCGATCAGCGACAAGGGGAGCAAAACCt                                                                                                                   >  2:244512/1‑251 (MQ=255)
                                                                                                               tGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCt                                                                                                                   >  2:237434/1‑251 (MQ=255)
                                                                                                                gAAAACGGCGAAGCGAAGAAAGTAAAAGCGATGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACATGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCt                                                                                                                   <  1:93344/250‑1 (MQ=255)
                                                                                                                gAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGCTGAAATCGATAAGCGACAAGGTGAGCAAAACCt                                                                                                                   >  2:147666/1‑250 (MQ=255)
                                                                                                                     cGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCACGTCTACAGCTCGAATATGGAAGCCTTATCGGATTAAGTTGGGCTGTTAACGCGAAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCACTACAGGACGCGGGATAGGCCACAGATCCTCAATATGCCCGCAACCTTACCAACCTGATTCAGCAGAGGAAATCGATAAGCGACAAGGAGAGCaagacataaagt                                                                                                              >  2:420025/1‑240 (MQ=255)
                                                                                                                       gcgaagcgaaGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATg                                                                                                           >  1:47168/1‑251 (MQ=255)
                                                                                                                               aagaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTg                                                                                                    >  2:219179/1‑250 (MQ=255)
                                                                                                                                    aGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGGTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGAGCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGAGGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATc                                                                                              <  1:222020/251‑1 (MQ=255)
                                                                                                                                                          tACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCcg                                                                        >  1:489243/1‑251 (MQ=255)
                                                                                                                                                           aCAGCTCGTATCTGGGAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGAGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAACCTCACCAACATGATTCAGCCGATGAAATCAATAAGCGCAAAGGGGAGCAAAACCTACAGTATGACAATTGATAATCTGTTCAGAATAACTCAAATCcgg                                                                       >  2:355095/1‑251 (MQ=255)
                                                                                                                                                           aCAGCTCGTAGCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTcc                                                                         <  1:333340/249‑1 (MQ=255)
                                                                                                                                                            cAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCcggc                                                                      >  1:221283/1‑251 (MQ=255)
                                                                                                                                                                           aGCCTTGTCGGATTACGTTGGGCGGTTAACGCGTAACCCGCGCTACGGCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAGCTGTTCTGAAGAACTCAAGTCCGGCGGGTCGCTGCCGata                                                       <  1:15884/251‑1 (MQ=255)
                                                                                                                                                                                           gTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGAGCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGaa                                       <  1:309634/251‑1 (MQ=255)
                                                                                                                                                                                            ttgggtcgTTAACGCGCAACCCTCGCTACGCCGCCGTGACGAGCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAATGTGAGCAAGACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGaa                                       <  1:278386/243‑1 (MQ=255)
                                                                                                                                                                                             tGGGCTGGTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAg                                      <  1:266050/250‑1 (MQ=255)
                                                                                                                                                                                                  ggttaccccgTAACCCTCGCTACTCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCCGGCCCTACAGGACGCCGTCTATGCCACCGATCCTCACTATGCCCGCAAACGCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCTATAATACTCTGTAATTGAATGCTtat                                <  2:298577/243‑1 (MQ=255)
                                                                                                                                                                                                                ccGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATg                  >  1:488505/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAg             <  1:408284/249‑1 (MQ=255)
                                                                                                                                                                                                                         gccgccGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTg         <  1:221340/251‑1 (MQ=255)
                                                                                                                                                                                                                         gccgccGTGACGAGCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTg         <  2:47168/251‑1 (MQ=255)
                                                                                                                                                                                                                           cgccgTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGAt       <  1:36820/251‑1 (MQ=255)
                                                                                                                                                                                                                                tGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATtaat   >  1:399025/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gacgatCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATACTACTCTGTAATTGAAGGCTTATAAGGAACCTCCATTGCCAGCTTGATtaata  >  2:104547/1‑250 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
GGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATA  >  NC_000913/1137922‑1138396

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