Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,725,582 C→A C44F (TGC→TTC)  ysaB ← putative lipoprotein YsaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,725,5820CA98.6% 257.1 / NA 71C44F (TGC→TTC) ysaBputative lipoprotein YsaB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (31/39);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATTATTCC  >  NC_000913/3725338‑3725827
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
gtgtGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACg                                                                                                                                                                                                                                                 >  1:283407/1‑251 (MQ=255)
    ggATTTAAAGGAAGGGATCAGTATACACGCTGGAGGGAAGATATACAGTACAAAACGGGAAAAAGCATGTCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATCTCATAGCTTCCCTGGAACAGTTCGAAGGCGGTGACGTCTATTTTTTGCTCGCCTGTGTTATAACGGTTTGCCGCCTGATCCTTGAATAACGCatc                                                                                                                                                                                                                                             <  2:44162‑M1/251‑4 (MQ=255)
              gaaggGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAACGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggat                                                                                                                                                                                                                                   <  1:274674‑M1/251‑14 (MQ=255)
                 gggATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcc                                                                                                                                                                                                                                >  2:582733‑M1/1‑235 (MQ=255)
                    aTCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCat                                                                                                                                                                                                                                              >  1:232074‑M1/1‑232 (MQ=255)
                    aTCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCat                                                                                                                                                                                                                                              <  2:232074‑M1/234‑3 (MQ=255)
                          ataCCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatc                                                                                                                                                                                                                        <  1:582733‑M1/250‑25 (MQ=255)
                          ataCCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAAGGGCCATCCGCGTCTAAAGAACATACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTGCATAGCTTCCCTGGAACGGTTCGAGTGCGTTGACGTTGATTTTTGGCTCGCCGGTGTTATAGCGTGTTGCCTACTGATCCTTGAAGAACGCTTcggtctgggagtcacgtttatc                                                                                                                                                                                                                        <  2:255041/250‑22 (MQ=255)
                            aCCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagt                                                                                                                                                                                                                     >  1:363645‑M1/1‑224 (MQ=255)
                                     atggaAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgttttt                                                                                                                                                                                                            >  1:22501‑M1/1‑215 (MQ=255)
                                          aGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCGTGAATAACGCagcggcctcggattcacgtttatcagtgttgtttttgggct                                                                                                                                                                                                       <  1:345111‑M1/251‑42 (MQ=255)
                                                    aacaaaacgggAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAAGAACTCatcggcctcggattcacgtttatcagtgtggtttgtgggctggcagccaga                                                                                                                                                                                             <  2:354369‑M1/247‑52 (MQ=255)
                                                          aCGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaaggg                                                                                                                                                                                         >  1:138095‑M1/1‑194 (MQ=255)
                                                            gggAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtc                                                                                                                                                                                     >  1:577397‑M1/1‑192 (MQ=255)
                                                                 aaaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattaacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggct                                                                                                                                                                                >  2:147743‑M1/1‑187 (MQ=255)
                                                                  aaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctg                                                                                                                                                                               <  1:26540‑M1/251‑66 (MQ=255)
                                                                     cAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGGCGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgata                                                                                                                                                                            <  2:464569‑M1/251‑69 (MQ=255)
                                                                       gggCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatatt                                                                                                                                                                          >  1:138108‑M1/1‑181 (MQ=255)
                                                                                gCATGGGAAGATACAACAAATGGTCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTAAAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcattgttgtttttgggctggcagccagaagggagtacggctgatattttgacaat                                                                                                                                                                  >  2:374025‑M1/1‑172 (MQ=255)
                                                                                gCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATATCTCatcggcctcggattcccgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaata                                                                                                                                                                 <  1:380525‑M1/251‑80 (MQ=255)
                                                                                         gATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggtt                                                                                                                                                        <  1:586088‑M1/251‑89 (MQ=255)
                                                                                            gCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTCCGGAAGCTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGCTTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgc                                                                                                                                                     >  1:243891‑M1/1‑160 (MQ=255)
                                                                                            gCAAAAAAGGGCCATCCGCGTCAAAAGAACAGACAAACGGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatatttttacaataatccgggttcg                                                                                                                                                      <  2:100598‑M1/250‑91 (MQ=255)
                                                                                              aaaacaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTAAGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcg                                                                                                                                                    >  2:22497‑M1/1‑158 (MQ=255)
                                                                                               aaaaaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCACTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGCTCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttctcgat                                                                                                                                                  >  2:148838‑M1/1‑157 (MQ=255)
                                                                                                aaaaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgatt                                                                                                                                                 <  2:235745‑M1/251‑96 (MQ=255)
                                                                                                   aTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctc                                                                                                                                              <  1:266726‑M1/251‑99 (MQ=255)
                                                                                                        cATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaag                                                                                                                                                                                           >  1:444159‑M1/1‑148 (MQ=255)
                                                                                                        cATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaag                                                                                                                                                                                           <  2:444159‑M1/201‑54 (MQ=255)
                                                                                                         aTCCGCGTAAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTGGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccata                                                                                                                                        <  2:249017‑M1/251‑105 (MQ=255)
                                                                                                           ccgagtaaaaagaccagaCAAACTGCGTTTTACGGAAGGTATAACCGCGCATTTCAAAGCTTCCCGGGAACTGTTCTAGGGGGTGGACGTCGATTTTTTGCTCGCCGGTGTTATACCGTTTGGCCGCCTGATCCTTGAAGAACGCatcggccgcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataac                                                                                                                                      <  2:471961‑M1/244‑107 (MQ=255)
                                                                                                           ccGCGTAAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTGGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataac                                                                                                                                      <  2:138108‑M1/251‑107 (MQ=255)
                                                                                                             gcgTCAACAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacac                                                                                                                                    >  1:435257‑M1/1‑143 (MQ=255)
                                                                                                                       aaCAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttttttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataa                                                                                                                          <  1:98975‑M1/251‑119 (MQ=255)
                                                                                                                              aaaCTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataattttt                                                                                                                     >  2:322813‑M1/1‑126 (MQ=255)
                                                                                                                              aaaCTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataattttt                                                                                                                     >  1:261822‑M1/1‑126 (MQ=255)
                                                                                                                               aaCTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagcc                                                                                                                                                                                                >  1:563454‑M1/1‑125 (MQ=255)
                                                                                                                               aaCTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagcc                                                                                                                                                                                                <  2:563454‑M1/173‑49 (MQ=255)
                                                                                                                                 cTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCACAGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCTTGATCCTTGAATAAAGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagagg                                                                                                                >  1:361267‑M1/1‑123 (MQ=255)
                                                                                                                                  tGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttggggtggcagccagaagggagtcaggctgata                                                                                                                                                                            <  2:12814‑M1/190‑69 (MQ=255)
                                                                                                                                  tGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttggggtggcagccagaagggagtcaggctgata                                                                                                                                                                            >  1:12814‑M1/1‑122 (MQ=255)
                                                                                                                                   gCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtg                                                                                                              <  2:417265‑M1/251‑131 (MQ=255)
                                                                                                                                     ccttTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgat                                                                                                            <  2:283407‑M1/250‑133 (MQ=255)
                                                                                                                                         tACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgag                                                                                                         <  1:93706‑M1/250‑136 (MQ=255)
                                                                                                                                            ggAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgc                                                                                                     >  2:175748‑M1/1‑112 (MQ=255)
                                                                                                                                              aaGGTATATCCGCGCATTTCATAGCTTCACTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTTTTATAACGTGTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcgggctgatagtttgacaataatccgggttcgcgattctcgccataacaccaaagaataacttttagaggtgatgagttgctt                                                                                                   <  1:285382‑M1/251‑142 (MQ=255)
                                                                                                                                              aaGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgctt                                                                                                   <  1:175748‑M1/251‑142 (MQ=255)
                                                                                                                                                     aaCCGAGCATTTCATAGCTTCCCTGGACCTGTTCGAAGGCGGTGACGTAGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgcttctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttac                                                                                             >  2:296948‑M1/1‑103 (MQ=255)
                                                                                                                                                         gcgcATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGACTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataac                                                                                        >  1:428414‑M1/1‑99 (MQ=255)
                                                                                                                                                                     gCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCCGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacg                                                                            <  1:265564‑M1/251‑165 (MQ=255)
                                                                                                                                                                        tCCCTGGAACTGTTCGAAGGGGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGGTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacga                                                                           <  2:433887‑M1/249‑166 (MQ=255)
                                                                                                                                                                         cccTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaagg                                                                        >  2:138565‑M1/1‑83 (MQ=255)
                                                                                                                                                                               acctgTTCGAAGGCGGTGACTTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataac                                                                  <  2:261822‑M1/249‑175 (MQ=255)
                                                                                                                                                                                    ttCGACGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggcttcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgct                                                             >  2:343404‑M1/1‑72 (MQ=255)
                                                                                                                                                                                          ggCGGTGCCGTCGATTTTTTGCTCGCCGGTGTTATCACGTTGTGCCGACTGGTCCTTGCATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgcttcgcgccataacacaaaagaataatttttagaggtgatgagttgcttagttacatgacgattgtatgacgacggcataacatgctgt                                                           >  2:389402‑M1/1‑66 (MQ=255)
                                                                                                                                                                                             ggTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcac                                                    >  1:412618‑M1/1‑63 (MQ=255)
                                                                                                                                                                                                  gcgTCGATTTTTTGCTCGCCGGTGGTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatca                                                <  2:284341‑M1/250‑193 (MQ=255)
                                                                                                                                                                                                       gATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtg                                            >  2:34928‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                              tGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgt                                    >  2:237846‑M1/1‑47 (MQ=255)
                                                                                                                                                                                                               gCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgta                                   <  1:34928‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                               gCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgt                                    <  1:337137‑M1/250‑205 (MQ=255)
                                                                                                                                                                                                                  cGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccggggtcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgctgagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaaga                                <  1:343404‑M1/251‑209 (MQ=255)
                                                                                                                                                                                                                     cGGTGTTATAACGTTTTGCCGCCTGATCCTGGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggtgcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaa                             <  2:513588‑M1/251‑212 (MQ=255)
                                                                                                                                                                                                                              aaCGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgc                     <  2:412618‑M1/250‑220 (MQ=255)
                                                                                                                                                                                                                                cGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgat                  >  2:160328‑M1/1‑29 (MQ=255)
                                                                                                                                                                                                                                   tttGCCGCCTGATCCTTGAATAACGCatcggcctcgtattcacgtttatcagtgttgtttttgggctggaagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggagaacgccgtaagaaaatatcttgtgattc                <  2:363645‑M1/250‑225 (MQ=255)
                                                                                                                                                                                                                                    ttGCCGCCTGATCCTTGAATAACGAatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattc                >  2:71633‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                     tGCCGTCTGATCCTTGAATAATGCatcggcctcggatttacgtttatcagtgttgtttttgggctggcagccagaaggtcttcatgctgatattttgacaataatccgggttcgcgattctcgccagaacaccaaagaataatttttagaggtgatgagttgcttagttgcataacgattgtatgacgaaggcataacatgttgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcag              <  1:310961‑M1/250‑227 (MQ=255)
                                                                                                                                                                                                                                       ccgccTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatc           <  1:237846‑M1/251‑230 (MQ=255)
                                                                                                                                                                                                                                            tGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaa       >  2:56005‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                                ccTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaaagatt  <  1:95370‑M1/251‑239 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATTATTCC  >  NC_000913/3725338‑3725827

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