Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,985,816 T→C intergenic (‑224/‑31) yqeF ← / → yqeG putative acyltransferase/putative transporter YqeG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,985,8160TC98.4% 216.9 / ‑7.3 61intergenic (‑224/‑31)yqeF/yqeGputative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/27);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGG  >  NC_000913/2985573‑2986062
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
aCAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGAGCAATAGTACATGATTGTTAATTACAGTTATATTAATTGTTTAAAAAAAGTGATTTTTATAACAAAGGAAATATGCCTGAGAAGCAGTCAGAGACATAACTGGCACGTAGGATTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGAGTTAATCACCCCTACTTAACCTGGTTTATGACAAACTGCACTCCATTTTGTTTa                                                                                                                                                                                                                                                  >  2:334843/1‑250 (MQ=255)
      acaACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGt                                                                                                                                                                                                                                            >  1:531575/1‑250 (MQ=255)
               ttCATAACAACATCGGAACAATGAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTGAATTGTTTAAAAAAAGTGATTTATATCAAAAAGGAAATACGCATGAGCAGAAGTCAGAGACATAGCGGGCACGTAGGGTTTGCAACCACTACCCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATGCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCCGTGTCATTTGTAgatgat                                                                                                                                                                                                                                     <  2:46236/248‑1 (MQ=255)
                  atcccatccTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGTCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTGTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGTTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTc                                                                                                                                                                                                                               <  1:531112/244‑1 (MQ=255)
                                  aaTTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                                >  1:148063‑M1/1‑235 (MQ=255)
                                     tAACAGGGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGGTTAAAAAAAGAGATTTTTATCACAAAGGAAATATGCATGAGCAGCAGTCAGAGACATAACTGGCACGTGGGGTTTGCAACCACTAACACACAAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTAGGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacattt                                                                                                                                                                                                            >  1:196285‑M1/1‑232 (MQ=255)
                                      aaCAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttg                                                                                                                                                                                                           <  1:481607‑M1/251‑21 (MQ=255)
                                         aGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccg                                                                                                                                                                                                        >  1:295217‑M1/1‑228 (MQ=255)
                                         aGTGCAATAGTAAATGATTGTTAACTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCGCCCATACTTAACCTGGTTTATGGTAAATTGCCCGCCATTCTGTTGAATTTGTAGATGATACGTTCtaaaaaacccccgacatgtgccg                                                                                                                                                                                                        <  1:102102‑M1/251‑24 (MQ=255)
                                          gTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgg                                                                                                                                                                                                       >  1:184875‑M1/1‑227 (MQ=255)
                                                    aaaTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgag                                                                                                                                                                                              <  1:55749‑M1/250‑34 (MQ=255)
                                                    aaaTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCGGAGCAGCAGTCAGAGACATAACTTGCACGTAGGGTTTTCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTTGTTTATGGTAAATTGCCCTCCATTCTGTTTTATTTGTAGATGATACGTTCtagaaaaaccccgacatatgccgggggtgtgagc                                                                                                                                                                                             <  2:35805‑M1/251‑35 (MQ=255)
                                                    aaaTGATTGTTAATGACAGTTATCTGAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCGCCAGTAGAGGGGTAGATAGGGTGTTAATCTCCCATACTTAACCTTGTTTATGGTAAATTGCCCTCCATTCTTTTTAATTTGTAGATGATACGTTCgaaaaaaaccccgacatttgccggggttgttagc                                                                                                                                                                                             <  2:176434‑M1/251‑35 (MQ=255)
                                                        gATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAgat                                                                                                                                                                                                                                        <  1:458879/204‑1 (MQ=255)
                                                        gATTGTTAATTACAGTTAACTTAATTGTTTAACAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGAAAATTGCCCTCCATTCTGTTTAATTTGTAgat                                                                                                                                                                                                                                        >  2:458879/1‑204 (MQ=255)
                                                                    cAGTTATCTTAATTGTTTAAAAAAAGAGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttat                                                                                                                                                                             >  2:225683‑M1/1‑201 (MQ=255)
                                                                      gTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattt                                                                                                                                                                           <  2:157397‑M1/251‑53 (MQ=255)
                                                                                      aaaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgcc                                                                                                                                                           >  2:138835‑M1/1‑183 (MQ=255)
                                                                                          aaaGTGGTTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACACACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggc                                                                                                                                                       >  2:120631‑M1/1‑179 (MQ=255)
                                                                                          aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggg                                                                                                                                                        >  1:398439‑M1/1‑179 (MQ=255)
                                                                                           aaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggga                                                                                                                                                       >  1:116076‑M1/1‑178 (MQ=255)
                                                                                              tGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaac                                                                                                                                                     >  1:405064‑M1/1‑175 (MQ=255)
                                                                                               gattgtcaccacaaaggaaataTGCCTGAGCAGCAGGCAGAGACATAACGGGCTCGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaacccccgacatttgccggggttgtgagcataacgtaatgcggattttaccggacgcatcgccgggaacaga                                                                                                                                                  <  2:408052‑M1/242‑78 (MQ=255)
                                                                                                aTTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacaga                                                                                                                                                  >  1:367615‑M1/1‑173 (MQ=255)
                                                                                                     tATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGCCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGCAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattntaccggacgcaccgccgggaacagaataacg                                                                                                                                            >  2:238785‑M1/1‑168 (MQ=255)
                                                                                                          caAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGTGCGTTAATCTCCCATACTTATCCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcg                                                                                                                                       <  1:126428‑M1/251‑89 (MQ=255)
                                                                                                          caAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcg                                                                                                                                       >  1:181412‑M1/1‑163 (MQ=255)
                                                                                                               gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                                                                                                                                  >  1:156790‑M1/1‑158 (MQ=255)
                                                                                                                  atatGCCTGAGCAGCAGTCAGAGACAGAACTGGCACGTAGTGTTTGCAACCACTAACCCACCAAGAGAGGGGTAGATAGGGCGTTAAGCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatttgaccggacgcatcgccgggaacagaataacggcgcggatggtat                                                                                                                               <  2:49080‑M1/251‑97 (MQ=255)
                                                                                                                    atGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatgg                                                                                                                             >  1:320243‑M1/1‑153 (MQ=255)
                                                                                                                              agcagTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtttgtg                                                                                                                    >  2:127821‑M1/1‑143 (MQ=255)
                                                                                                                                cagTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaa                                                                                                                  >  1:238833‑M1/1‑141 (MQ=255)
                                                                                                                                      gagaCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcatt                                                                                                           <  2:184875‑M1/251‑117 (MQ=255)
                                                                                                                                           aTAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccat                                                                                                      >  1:295218‑M1/1‑130 (MQ=255)
                                                                                                                                             aaCTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatac                                                                                                    <  2:156790‑M1/251‑124 (MQ=255)
                                                                                                                                              aCTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCCTTCTGTTTAATTGGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacatcattaccatacg                                                                                                   <  2:368390‑M1/251‑125 (MQ=255)
                                                                                                                                                      gtagggtttgCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgatac                                                                                           <  1:349847‑M1/247‑133 (MQ=255)
                                                                                                                                                        agtgtttGCAAACACTAACCCACCAATAGATGGTTAGATAGGGCTGTAATCTCCCATACTTATCCTTGTTTATGGTAAATTGCCCTCCAGTCTGTTGAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacggaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcatgtccatccggtcgatacct                                                                                         <  2:156656‑M1/248‑135 (MQ=255)
                                                                                                                                                        agggttTGCAACCACTAACCCACCACTAGAGGGGTAGATAGGGCGTTAATCTCCCATAATTAACCTGGTTTATGGTAAATTGCCCTACATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattcacaaaaggtcgataccc                                                                                         >  2:491914‑M1/3‑117 (MQ=255)
                                                                                                                                                                     aCTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttgg                                                                                                                      <  2:497192‑M1/209‑106 (MQ=255)
                                                                                                                                                                     aCTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttgg                                                                                                                      >  1:497192‑M1/1‑104 (MQ=255)
                                                                                                                                                                                aaTAGTGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggt                                                                 <  2:192156‑M1/251‑159 (MQ=255)
                                                                                                                                                                                        ggTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacccgctgtcggcggtaaggcctt                                                         >  2:213245‑M1/1‑85 (MQ=255)
                                                                                                                                                                                            gATAGGGCGTTAATCGCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttcc                                                     <  2:295218‑M1/251‑171 (MQ=255)
                                                                                                                                                                                                       aaTCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtga                                           <  1:485628‑M1/250‑181 (MQ=255)
                                                                                                                                                                                                                 aCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttc                                 >  1:491915‑M1/1‑60 (MQ=255)
                                                                                                                                                                                                                    tAACCTGGTTTATGGTAAATTGCCCTCCATTCTGCTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggaatggctgttggtgaacagccttaccatacgggtgataccaattaccagacccgctgtccgcggtaagccaagtttccgtgaggggaagtaatcttcctcg                             >  2:462404‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                                     aaCCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTATATGATACGTTCtaaaaaaaccccgacattagccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaagcttcatcgt                            <  2:89632‑M1/251‑196 (MQ=255)
                                                                                                                                                                                                                          ggTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTTTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggatcagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaa                        <  1:343594‑M1/250‑200 (MQ=255)
                                                                                                                                                                                                                           gTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagcac                       >  1:319563‑M1/1‑50 (MQ=255)
                                                                                                                                                                                                                             ttATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacat                     <  1:138835‑M1/250‑203 (MQ=255)
                                                                                                                                                                                                                             ttATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTGCtaaaaaaacccctacatttgccggggtggtgagcataccgtaatgcttatttttccggacgcatcgccgggaacagaataacgtcgctgattgtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacat                     <  2:356985‑M1/250‑203 (MQ=255)
                                                                                                                                                                                                                               aTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgccatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaatcacgtcgcggatggtatggatgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcttagaacatcg                   <  1:481197‑M1/250‑205 (MQ=255)
                                                                                                                                                                                                                                  gTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctg                                                                         >  1:288710‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                                                  gTAAATTGCCCTCCATGCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgccatttgccggggttgtgagcataacgtaatgcttattttaccggacgcagcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccagacggtcgataccaattcccagacctgctg                                                                         <  2:288710‑M1/193‑151 (MQ=255)
                                                                                                                                                                                                                                   tAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaacccagacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgaaaccaattcccagacctgcggtcggcggtaagccaggttccagtctggngaacgaaacttcatctgagaacatcgctttc              >  2:240004‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                                                      aTTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattcccatacggtcggtaccaattcccagacctgatgtcggcggtaagcccagttcccgtgcggtgcnggcattttcatcgtagaaaaccgcttcgtcc           >  2:233637‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                                       ttgccatacattctgttTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtttgtgaacagcattaccatacggtcgataccaatttccagaccggcaggcggcggtaagccatgtttcagtgcggtgtagaaaacttcatcgcagaacaacggctcgtcgc          >  1:255455‑M1/9‑38 (MQ=255)
                                                                                                                                                                                                                                        tGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggggccgtaatcttcctcgtagaacatcgcttcgtcggc         >  2:389561‑M1/1‑37 (MQ=255)
                                                                                                                                                                                                                                          cccTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacattttccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctggtggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgtttcagtgcggtgacgtaatcttcctcgtagaacatctcttcctcggccc       >  2:15253‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                                            cTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtcgtcacct     >  1:255012‑M1/1‑33 (MQ=255)
                                                                                                                                                                                                                                               cATTCTGTTTACTTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggattgaatggctgttggtgaacagcattaccatacggtcgataccaattaccagacctgctgtcggcggtaagacatgttacagtgcggtgccggaatcttcctcgtagaacaacgcgtaggtatcaccttcg  >  2:160793‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
ACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAGGG  >  NC_000913/2985573‑2986062

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