Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,628,332 T→C intergenic (+69/‑20) ydfG → / → rspR 3‑hydroxy acid dehydrogenase YdfG/DNA‑binding transcriptional repressor RspR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,628,3320TC95.0% 204.3 / NA 60intergenic (+69/‑20)ydfG/rspR3‑hydroxy acid dehydrogenase YdfG/DNA‑binding transcriptional repressor RspR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/34);  minor base A (2/1);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACTCACTCACAAGCAAGAACGAATGACCGTCGAAACGCAACTTAATCCCACACAGCCTGTCAATCAGCAGATTTATCGTATTCTTCGTCGCGACATTGTCCATTGCCTGATTGCTCCAGGCACACCGTTGTCGGAAAAAGAAGTTTCTGTTCGTTTCAATGTGTCACGCCAGCCGGTTCGTGAAGCCTTTATTAAACTGGCGGAAAACGGCCTGATTCAAATTCGTCCGCAACGTGGCAGCTAC  >  NC_000913/1628099‑1628573
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
ggTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaac                                                                                                                                                                                                                                   >  1:573595‑M1/1‑241 (MQ=255)
               ccctatcaaaagcactttgcattgacgccagaagATGGCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGTGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACACACTCACtacatgaacagcggcaggaccttta                                                                                                                                                                                                                   <  2:152041‑M1/235‑26 (MQ=255)
                 cTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttagg                                                                                                                                                                                                                 >  1:400257‑M1/1‑224 (MQ=255)
                      aaaatccagttgcattgACGGCAGAAGATGTCAGCGAAGCCGTCTGGTGGGAGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCAATGCCGGACTGAATGTCAACCGTCAGGAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTTTCACAAAACAGTGGAAGACTCATGCAGTTAACACACACACtatatgaacagcggcaggacctttagggcagt                                                                                                                                                                                                            >  2:603284‑M1/9‑219 (MQ=255)
                          tACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAACGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggaactttaggtccgttctc                                                                                                                                                                                                        >  2:420086‑M1/1‑215 (MQ=255)
                              gTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaat                                                                                                                                                                                                     >  1:298716‑M1/1‑211 (MQ=255)
                                tGCATTGACGCCAGAAGATGTCAGCGCAGCCGTCTGGTGGGTGTCAACGCTTCCTGCTCACGTCAATAACAATTCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacttgaacagcggcaggacctttaggtccgttctcaatcc                                                                                                                                                                                                   <  2:220037‑M1/250‑42 (MQ=255)
                                  cATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagc                                                                                                                                                                                                                                <  2:523103‑M1/219‑13 (MQ=255)
                                  cATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagc                                                                                                                                                                                                                                >  1:523103‑M1/1‑207 (MQ=255)
                                   atttacgccAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaa                                                                                                                                                                                               <  2:448839‑M1/247‑46 (MQ=255)
                                      gACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTCTTGCTTGTCACAACAAAGTGGTAGCCTAAGGCAGTTCACCCACTCACtacatgaacaacggtaggcacattaggtacggtctcaatacaaaattc                                                                                                                                                                                            >  2:741389‑M1/1‑203 (MQ=255)
                                          ccAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaact                                                                                                                                                                                        <  2:210245‑M1/251‑53 (MQ=255)
                                           cAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaactt                                                                                                                                                                                       <  1:51589‑M1/251‑54 (MQ=255)
                                                  tGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGCCTCATGCAGTTAACCCACTCACtacatgaacagcggcaggccctttaggtccgttctcaataccaaattcaacttcctgctt                                                                                                                                                                                >  2:106858‑M1/1‑191 (MQ=255)
                                                   gTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattc                                                                                                                                                                               <  1:340909‑M1/251‑62 (MQ=255)
                                                      aGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcc                                                                                                                                                                            >  2:161963‑M1/1‑187 (MQ=255)
                                                         gtagcCGTCTGGTGGGTGTCAACACTGCCTGCTCACGTCAATATCAACACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCATCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaattccaaattcaacttcctgattctcagtt                                                                                                                                                                         <  1:25291‑M1/249‑68 (MQ=255)
                                                         gAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctnaataccaaattccacttcctgattctcagtt                                                                                                                                                                         >  2:360309‑M1/1‑184 (MQ=255)
                                                         gAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagtt                                                                                                                                                                         >  1:163290‑M1/1‑184 (MQ=255)
                                                           aGCCTTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCTGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCTCTCACtacatgaacagcgacaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaa                                                                                                                                                                       <  1:449277‑M1/251‑70 (MQ=255)
                                                            gCCGTATGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctccataccaaattcaacttcctgattctcagttaat                                                                                                                                                                      >  2:445866‑M1/1‑181 (MQ=255)
                                                               gTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACACTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGAAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACACAAAAGTGGTAGACTCATGCAGTTAACACACTCACtacatgaacagcggcaggacctttcggtacgttctcaataccaaattcaacttcctgaatctccgttaatgtt                                                                                                                                                                   >  2:62913‑M1/1‑178 (MQ=255)
                                                                         tgtCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctccgttaatgttttgaaatcgt                                                                                                                                                         >  1:518604‑M1/1‑168 (MQ=255)
                                                                             aaCGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgct                                                                                                                                                     >  2:250470‑M1/1‑164 (MQ=255)
                                                                                    cctgcTCACGTCAATATCAATCCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAAGTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAATCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaat                                                                                                                                               <  2:486501‑M1/250‑94 (MQ=255)
                                                                                       gcTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattg                                                                                                                                             <  2:298716‑M1/249‑96 (MQ=255)
                                                                                       gcTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGGTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaagtcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgct                                                                                                                                           <  1:39935‑M1/251‑98 (MQ=255)
                                                                                          cACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttcccaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctg                                                                                                                                          >  2:486671‑M1/1‑151 (MQ=255)
                                                                                                       ataCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaaca                                                                                                                           <  1:161963‑M1/251‑114 (MQ=255)
                                                                                                              ggAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatcttt                                                                                                                    <  1:36626‑M1/251‑121 (MQ=255)
                                                                                                               gAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttg                                                                                                                   >  1:191952‑M1/1‑130 (MQ=255)
                                                                                                                    gatgCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgcc                                                                                                              <  1:748370‑M1/251‑127 (MQ=255)
                                                                                                                         cGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctt                                                                                                         <  1:500920‑M1/251‑132 (MQ=255)
                                                                                                                            ttACCCAAAGCGATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTATACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttc                                                                                                       <  2:437047‑M1/250‑134 (MQ=255)
                                                                                                                                 cAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtg                                                                                                 >  1:176062‑M1/1‑112 (MQ=255)
                                                                                                                                  aaaGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtg                                                                                                 >  2:149269‑M1/1‑111 (MQ=255)
                                                                                                                                     gCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatga                                                                                             <  2:191952‑M1/251‑144 (MQ=255)
                                                                                                                                          gCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgtgaaatagacaaacacatctttgctgccatctttcggcgtgatgaaacca                                                                                        <  1:420086‑M1/251‑149 (MQ=255)
                                                                                                                                          gCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaacc                                                                                         <  1:381849‑M1/250‑148 (MQ=255)
                                                                                                                                          gCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaacc                                                                                         >  1:4666‑M1/1‑103 (MQ=255)
                                                                                                                                            cGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGATATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatggtttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaa                                                                                       >  1:73201‑M1/1‑101 (MQ=255)
                                                                                                                                                          aCCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggt                                                                        <  2:4666‑M1/251‑165 (MQ=255)
                                                                                                                                                                    atgtttatgCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGCCTCATGCAGTTAACCCACTCACtgcatgaacagcggcatgtcctttaggtccgttctcaataccaaattcaccttcctgatgctcagttaatgttgtgaaatcgtggctctgaactgctgagaaatagacagacacatcttttctgccatctttcggcgtgatgaatccaaaacctttttcagggttaaaccattt                                                              <  1:434476‑M1/248‑175 (MQ=255)
                                                                                                                                                                    aTTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccattt                                                              <  1:268816‑M1/251‑175 (MQ=255)
                                                                                                                                                                      ttttATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttc                                                             <  1:445866‑M1/250‑176 (MQ=255)
                                                                                                                                                                      ttttATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttc                                                             <  2:715501‑M1/250‑176 (MQ=255)
                                                                                                                                                                         tataCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcact                                                          <  1:62913‑M1/250‑179 (MQ=255)
                                                                                                                                                                            aCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaac                                                      >  2:495162‑M1/1‑69 (MQ=255)
                                                                                                                                                                                ggcttaactGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagataatgttttgaaatagttgctctgaattgatgagaaatagacaaacacatctttgctgccctctttcggcgttatgaaaccaaaacctttttcaggtttaaacaatttccataaaccagt                                                  >  2:176861‑M1/5‑65 (MQ=255)
                                                                                                                                                                                     aaCTGCCGGGTTGTTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattt                                             <  1:293711‑M1/251‑192 (MQ=255)
                                                                                                                                                                                         gCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgtt                                         <  1:592978‑M1/251‑196 (MQ=255)
                                                                                                                                                                                            gggTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttaga                                      <  1:501734‑M1/251‑199 (MQ=255)
                                                                                                                                                                                                 aTTGCTTTTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataa                                 <  2:674407‑M1/251‑204 (MQ=255)
                                                                                                                                                                                                 aTTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataa                                 <  1:486671‑M1/251‑204 (MQ=255)
                                                                                                                                                                                                   tGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgtggctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataat                                <  1:603284‑M1/250‑205 (MQ=255)
                                                                                                                                                                                                                  aGTGGTAGACTCATGCAGTTAACCCACTCACnacaggaacatcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaaga                <  2:340027‑M1/251‑221 (MQ=255)
                                                                                                                                                                                                                        aGACTCATGCAGTTAATCCACTCACtacatgaacagcggcaggacctttagggccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctttgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgaggaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacttaattatttccttttgaagaaattag          <  2:487457‑M1/251‑227 (MQ=255)
                                                                                                                                                                                                                        aGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaagaaattag          <  1:360309‑M1/251‑227 (MQ=255)
                                                                                                                                                                                                                         gACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaagaaattagc         >  1:237460‑M1/1‑24 (MQ=255)
                                                                                                                                                                                                                                gCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaagaaattagcccttggg  >  1:541481‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
GGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACTCACTCACAAGCAAGAACGAATGACCGTCGAAACGCAACTTAATCCCACACAGCCTGTCAATCAGCAGATTTATCGTATTCTTCGTCGCGACATTGTCCATTGCCTGATTGCTCCAGGCACACCGTTGTCGGAAAAAGAAGTTTCTGTTCGTTTCAATGTGTCACGCCAGCCGGTTCGTGAAGCCTTTATTAAACTGGCGGAAAACGGCCTGATTCAAATTCGTCCGCAACGTGGCAGCTAC  >  NC_000913/1628099‑1628573

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