Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,130 C→T P254S (CCG→TCG)  flgJ → putative peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,1300CT96.5% 142.9 / NA 59P254S (CCG→TCG) flgJputative peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (27/29);  minor base G (0/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAG  >  NC_000913/1137885‑1138369
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
atTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCgcgc                                                                                                                                                                                                                                             >  1:461778/1‑251 (MQ=255)
atTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCgcgc                                                                                                                                                                                                                                             >  1:328580/1‑251 (MQ=255)
      tCCTCGCGCAGGCTGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTCTCTGGTAGCCTTGTCGGATTACGTTGGGATGTTAACGCGTAACTCGCGCTCcgcc                                                                                                                                                                                                                                       >  1:2077/1‑251 (MQ=255)
      cccTCGCTCAGGCGGCCCTGGAATTTGGTGGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTAcgcc                                                                                                                                                                                                                                       <  2:511495/250‑1 (MQ=255)
           ccccagccggcACGGGAATTTGGTGGGGGGAAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAACCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCTCGCTACGCCGCCGt                                                                                                                                                                                                                                  <  2:563106/244‑1 (MQ=255)
                  cggcACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGAcc                                                                                                                                                                                                                           >  1:536369/1‑251 (MQ=255)
                         ggattttgtttGGGGGCAAGGGCAAACCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCCCCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTAAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGa                                                                                                                                                                                                                    <  2:177084/242‑1 (MQ=255)
                                  ttGGGGGCAACGGCAACTCCGCCGCGAAAACGGCGAGACGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTAAAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGAGCTACGCCGCCGTGACGCCCGCCGCGCATGCGGAAc                                                                                                                                                                                                           >  2:396538/1‑251 (MQ=255)
                                       ggCACCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGGGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAAGCACCACGACGGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTCCAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAgggg                                                                                                                                                                                                      <  1:354493/251‑1 (MQ=255)
                                       gacaccGGCAAACCCCCCGCAAAAACGGCGAGCCGAGCTATACCCTGTTTGGTGTCAAAGCCTCTGGCAACGGGAAAGGGCCAGTTACGGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAGGCGAAGTGGCGCGTCTACAGCTCGTGTCTGGAAGCCTTGTCGGATTACGTTGGGCTTTTAACGCGTAACTCGCGCTACGCCGCTGTGACGGCCGCCGCTAGTGCGGAAGAgggg                                                                                                                                                                                                      <  2:87233/246‑1 (MQ=255)
                                         cagcgtggagtCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACAGCGTAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTCCAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGCGTGCGGCACAGGTGgc                                                                                                                                                                                                    >  2:182859/11‑251 (MQ=255)
                                                        cgcgAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAgg                                                                                                                                                                                      >  1:358054/1‑250 (MQ=255)
                                                        cgcgAAAACGGCGAGCCGAGCTATAACCTGTGTGGTTTCAAAGCCTCTGGAAACTGGAAAGGGCCAGTTACGGAAATCTCCACGACTGAATATGAAACCGGCGAAGCGAAGAAAGGAAAAGCGAAGTTCCGCGTCTACAGCTCGTGTCTTGAAGCCTTGTCGGATTACGTTGGGCGGTTAACTCGTAACCCGCGCTACTGCGCCGGTACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGa                                                                                                                                                                                     <  1:43076/251‑1 (MQ=255)
                                                         gcgAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAGCTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGCACAGGGGGCGCAGGCCCTACAGGAc                                                                                                                                                                                    >  2:236343/1‑251 (MQ=255)
                                                             aaaCGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGCGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCgg                                                                                                                                                                                >  2:520783/1‑251 (MQ=255)
                                                                                         ggtgacaagctcctgggCAACTGGAAAGGCACTGTTCTGTAATGACCCACGCCTGAAGTTGAAAATGGAGAATCGATGAAATTAAAAGCGAATTTTCGCGTTTCCAGCTCGTATCTGTAGGCCTTGCCGTATTACGTTGGGCTTTAACCCCGAAACTCGCGCTACGCCGCCGGGACGACCGCCGCGAGTGCGGAACAGGGGGCTCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCGCTCTATGCCCg                                                                                                                                                   <  2:17429/236‑1 (MQ=255)
                                                                                            gtCAAAGCCTCTGGCAGCGGGAAAGGGCCATGTACTGAAATCACCCCGTCTGAATGTGAAAAGGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCTTATCTGGAAGCCTTGTCGGATTACGTTGGCCTGTTAACGCGTACCTGGCGCTGCGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGGTATGCCACCGAGCCTCACTATGCCCGCaa                                                                                                                                                <  2:250386/251‑1 (MQ=255)
                                                                                              cAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAAc                                                                                                                                              >  2:436635/1‑251 (MQ=255)
                                                                                                   ccTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCCCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCAcc                                                                                                                                         <  1:345648/251‑1 (MQ=255)
                                                                                                      ctGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAAc                                                                                                                                      >  2:460624/1‑251 (MQ=255)
                                                                                                           aaCTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACTTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGGGCGGAACAGGGGGCGCAGGCACTACAGGGCGCGGGCTCTGCGACCGTACCTCAATATGCCCGCACAATCACACACATGAt                                                                                                                                 >  2:324690/1‑251 (MQ=255)
                                                                                                           aaCTGGAAAGGGACAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAAGGGGCGCAGGCCCTACAGGACGCGGGCTCTGCAACCGATCCTCACTATGCCCGCAAAATAACCAACATGAt                                                                                                                                 >  2:324705/1‑251 (MQ=255)
                                                                                                           aaCTGAAAAGGGCCAGTTACTGACATCAACACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGTCCAGGCCCTACACGACGCGGGCTATGCTCACGCCCCTCCCTATGCCAGCAAACTCACCAACATGAt                                                                                                                                 >  2:587338/1‑251 (MQ=255)
                                                                                                             atGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCTAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGTATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTc                                                                                                                               <  1:156533/250‑1 (MQ=255)
                                                                                                              gggAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCAAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCCCGTATCTGGAAGCCTTGTCGGATTACGTGGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTca                                                                                                                              <  1:520783/250‑1 (MQ=255)
                                                                                                                  aaGGTCCAGTTACTGAACTCACCCCGACTGAATATGAAAACGGCGAAGCGAACAAAGAAAAAGCGAGGTTTCGCGTCAACAGCCCGTATCTGGAAGCTTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTCCGCCGCCGTGACGCCCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTcagcag                                                                                                                          <  2:407098/251‑1 (MQ=255)
                                                                                                                     ggCCAGTTACTGAAATCCCCCCGATAGAATATGAAAGCGGCGAAGCGAAGAAAGTAAAAGCGAAGTTCCCGGTCTACAGCTCGTATCTGGAAGCCTGGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATg                                                                                                                       <  2:488627/251‑1 (MQ=255)
                                                                                                                                gAAATCACCACGACGGAATAGGAAAACGGCGAAGCGAAGAAAGGAAATGCGAAGTTGCGCGTCGACAGCTCGTATCTGTAAGCCTTGTCGGATTCCGTTGGGCTGTGAACGCGTAACTCGCGCTACGCCGCCGGGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGGCCTACAGGACGCGGGGTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGAGTTAACAGATGAAATCGATAAg                                                                                                            <  2:583441/251‑1 (MQ=255)
                                                                                                                                  aaTCACCACGACTAAATATGAAAACGGCAAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCg                                                                                                          <  2:423913/251‑1 (MQ=255)
                                                                                                                                       ccacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGGCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAg                                                                                                     <  2:164422/251‑1 (MQ=255)
                                                                                                                                                gAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGc                                                                                               <  1:517815/249‑1 (MQ=255)
                                                                                                                                                gAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCa                                                                                              >  1:486130/1‑250 (MQ=255)
                                                                                                                                                 aaTATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCaaa                                                                                            <  1:250797/251‑1 (MQ=255)
                                                                                                                                                            ggagaagcgaggAAAGTAAAAGGGAAGTTTCGCGACTACAGCTTGTTGCTGGAAGCTTTGCGGGATTACGTGGGGCTGTTACCGCGAAACTGGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAACCCTACAGTAt                                                                                 <  2:576056/248‑1 (MQ=255)
                                                                                                                                                             gcgaagcgaaGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATg                                                                                <  1:217431/251‑1 (MQ=255)
                                                                                                                                                                aagcgaagAAAGGAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAAc                                                                             >  2:462144/1‑251 (MQ=255)
                                                                                                                                                                          aGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATc                                                                   >  2:226345/1‑251 (MQ=255)
                                                                                                                                                                                agggAAGTTTCGCGTCTACAGCTCTTATCTGGAAGCCTTGGCGGATTACGTGGGGCTGTTAACGCGTAACTCGCGCCACCCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATGCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTTATTATCTGTTCt                                                             <  2:350949/248‑1 (MQ=255)
                                                                                                                                                                                 gCGAAGTTTCGCGTCTACATCTCGTATCGGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTg                                                            <  1:380061/251‑1 (MQ=255)
                                                                                                                                                                                            cgTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCa                                                   >  1:77352/1‑249 (MQ=255)
                                                                                                                                                                                             gTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGt                                                >  2:94658/1‑251 (MQ=255)
                                                                                                                                                                                             gTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGt                                                >  1:152638/1‑251 (MQ=255)
                                                                                                                                                                                                        gTATCTGGAAGCCTTGTCGGATTACGATGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGAGGGTCg                                     >  1:333180/1‑251 (MQ=255)
                                                                                                                                                                                                         tATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGc                                    >  1:288560/1‑251 (MQ=255)
                                                                                                                                                                                                          aTCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCt                                   >  2:114749/1‑251 (MQ=255)
                                                                                                                                                                                                           tCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCACGCGAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTACGGCGGGACGCtc                                  >  1:566324/1‑250 (MQ=255)
                                                                                                                                                                                                              ggAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCACGTCCGGCGGGTCGCTGCCg                               >  2:469803/1‑251 (MQ=255)
                                                                                                                                                                                                              gaaagcttggccGGATTACGTTGGTCTTTTAACGCGGAACGCGCGCTACCCCGCCGTGCCGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGGCCGGCGTGTCGCTGcc                                <  2:265493/239‑1 (MQ=255)
                                                                                                                                                                                                                       gTCGGATTACGTGGGGCTGTAACCGCGTAACTCGCGCACCCCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATAct                       <  2:486130/250‑1 (MQ=255)
                                                                                                                                                                                                                         cGGATTCCGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATActct                     <  1:253343/250‑1 (MQ=255)
                                                                                                                                                                                                                         cGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGTTAATACTCTg                    >  1:136763/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGGATTACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGGCGCTGCCGATAATACTCTg                    <  1:213801/251‑1 (MQ=255)
                                                                                                                                                                                                                            atttcgtggggctttgaccgcgTAACGCGCGCTACCCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTa                  <  2:328580/232‑1 (MQ=255)
                                                                                                                                                                                                                              tACGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAAt                <  1:460624/250‑1 (MQ=255)
                                                                                                                                                                                                                               aCGTTGGGCTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAAtt               <  1:392259/250‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTGGGTTGTTAACGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGaa            <  1:574344/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCTGTTAACGCGTAACTCGAGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCACGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGACTAACTCAAGTCCGGGGGGTCGCTGACGATAATACTCTGTAATTGAAGGCt        >  2:189355/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCTATAATACTCTGGAATTGAAGGCTTATAAg  <  1:165192/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCGCGTAACTCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAg  <  1:270005/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAG  >  NC_000913/1137885‑1138369

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