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,3320TC100.0% 227.9 / NA 63intergenic (+69/‑20)ydfG/rspR3‑hydroxy acid dehydrogenase YdfG/DNA‑binding transcriptional repressor RspR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/44);  total (19/44)

AAGGCGATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACTCACTCACAAGCAAGAACGAATGACCGTCGAAACGCAACTTAATCCCACACAGCCTGTCAATCAGCAGATTTATCGTATTCTTCGTCGCGACATTGTCCATTGCCTGATTGCTCCAGGCACACCGTTGTCGGAAAAAGAAGTTTCTGTTCGTTTCAATGTGTCACGCCAGCCGGTTCGTGAAGCCTTTATTAAACTGGCGGAAAACGGCCTGATTCAAATTCGTCCGCAACG  >  NC_000913/1628088‑1628563
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
aaGGCGATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAAcccactca                                                                                                                                                                                                                                   <  1:42978/251‑1 (MQ=255)
  agcggggtgaaaatggtggtactcgcctatcaaaataccgttgcattgacgcCAGAAGATCTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACt                                                                                                                                                                                                                                 <  1:311170‑M1/226‑2 (MQ=255)
          cGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaac                                                                                                                                                                                                                         <  2:304504‑M1/251‑10 (MQ=255)
             tAAAGCAGAAAAAACCTATCAAAATCCCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCCGTTAACCCACTCACtacatgaacagc                                                                                                                                                                                                                      >  1:506647‑M1/1‑239 (MQ=255)
                   aaaaaaaaccTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTGACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcagg                                                                                                                                                                                                                <  2:394826‑M1/249‑19 (MQ=255)
                    gAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcagga                                                                                                                                                                                                               <  1:193930‑M1/251‑20 (MQ=255)
                     aaaaaaCCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGCAGCCGTCTGGTGGGTGTCAACCCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCGACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGGTATTGCTTGTCACAAAAAAGTGGTAGACTAATGCAGTTAACCCACTCACtacatgaacagcggcaggac                                                                                                                                                                                                              >  1:545531‑M1/1‑231 (MQ=255)
                       aaaaCCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGTAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacct                                                                                                                                                                                                            <  2:566564‑M1/251‑23 (MQ=255)
                          aCCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggaccttta                                                                                                                                                                                                         <  1:370611‑M1/251‑26 (MQ=255)
                              aacaaAATACCGTTGCATTGAAGCCAGAAGATGTCAGCGAAGCCGGCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAAGTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtc                                                                                                                                                                                                     <  2:443114‑M1/249‑30 (MQ=255)
                              aTCAAAATACCGGTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggt                                                                                                                                                                                                      <  2:169712‑M1/250‑29 (MQ=255)
                                 aaaaTACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggaccgttaggtccgt                                                                                                                                                                                                  <  2:167045‑M1/251‑33 (MQ=255)
                                    aTACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttct                                                                                                                                                                                               >  2:376025‑M1/1‑216 (MQ=255)
                                      aCCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctc                                                                                                                                                                                              >  2:92892‑M1/1‑214 (MQ=255)
                                       ccGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaa                                                                                                                                                                                            <  2:577633‑M1/251‑39 (MQ=255)
                                        cGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaat                                                                                                                                                                                           <  1:249135‑M1/251‑40 (MQ=255)
                                                tGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaatt                                                                                                                                                                                   >  1:110978‑M1/1‑204 (MQ=255)
                                                  aCGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattca                                                                                                                                                                                 <  2:23531‑M1/251‑50 (MQ=255)
                                                      cAGAAGATGTCAGCGAAGACGTCTGGTGGGTGTCAACGCTGCCTGCTCACATCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTACGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacctgcacagaggcaggacctttaggtccgttctcaataccaaattcaactt                                                                                                                                                                             >  2:513708‑M1/1‑198 (MQ=255)
                                                              gTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattc                                                                                                                                                                     >  1:276589‑M1/1‑190 (MQ=255)
                                                                     aaGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagtta                                                                                                                                                              <  1:35075‑M1/251‑69 (MQ=255)
                                                                         cGTCTGGTGGGTGACAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTGAGTAATTTTTAGACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaatacaacttcctgattctcatttaatgt                                                                                                                                                          <  2:534112‑M1/251‑73 (MQ=255)
                                                                              ggtgggtgTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgtttt                                                                                                                                                       <  2:387841‑M1/248‑76 (MQ=255)
                                                                                       cAACGCTGCCTGCTCACGTCAATATCAATTCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTTCCGGGTTATTGCTTGTCACAAAAAAGGGGTAGACTCATGCAGTTAACCCACTCACtacatgatcagcggcaggtcctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttg                                                                                                                                             <  2:493684‑M1/250‑86 (MQ=255)
                                                                                           gccgcctgctcaCGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGGCCCCCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctg                                                                                                                                        <  1:390254‑M1/248‑91 (MQ=255)
                                                                                                ctgcTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattg                                                                                                                                   >  1:152738‑M1/1‑156 (MQ=255)
                                                                                                  gcTCACGTCAATATCAATACCCTGGAAATGATTCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATTCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgc                                                                                                                                  <  2:407232‑M1/250‑97 (MQ=255)
                                                                                                        gTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaa                                                                                                                            >  1:321748‑M1/1‑148 (MQ=255)
                                                                                                                caataCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaa                                                                                                                   <  1:302863‑M1/251‑112 (MQ=255)
                                                                                                                    aCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGACTGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacac                                                                                                                >  1:469801‑M1/1‑136 (MQ=255)
                                                                                                                    aCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacaca                                                                                                               <  2:263929‑M1/251‑116 (MQ=255)
                                                                                                                    aCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacaca                                                                                                               <  1:135558‑M1/251‑116 (MQ=255)
                                                                                                                    aCCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacac                                                                                                                <  1:136660‑M1/250‑115 (MQ=255)
                                                                                                                                atgCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgcca                                                                                                   <  1:603230‑M1/251‑128 (MQ=255)
                                                                                                                                    cGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctt                                                                                               <  1:376025‑M1/251‑132 (MQ=255)
                                                                                                                                     ggTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatcttt                                                                                              <  1:153223‑M1/251‑133 (MQ=255)
                                                                                                                                      gTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTCATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGGAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattccacttcctgattctcagttaatgttttgaaatcgttgctctgaatt                                                                                                                                    >  2:26048‑M1/1‑118 (MQ=255)
                                                                                                                                      gTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaatt                                                                                                                                    <  1:26048‑M1/212‑95 (MQ=255)
                                                                                                                                           ccAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTGTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgt                                                                                        <  2:507704‑M1/251‑139 (MQ=255)
                                                                                                                                                   aTGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCATTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaa                                                                                 <  1:247662‑M1/250‑146 (MQ=255)
                                                                                                                                                   aTGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaa                                                                                 >  1:9286‑M1/1‑105 (MQ=255)
                                                                                                                                                     gCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaacca                                                                              <  2:110978‑M1/251‑149 (MQ=255)
                                                                                                                                                             gAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaatactaaattcaacttcctgattcgcagttaatgttttgaaatcgttgctctgacttgctgagaaatagacaaacacatctttgcttccatctttcgtcgggatgaaaccaaaacct                                                                        <  2:373143‑M1/249‑155 (MQ=255)
                                                                                                                                                             gAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctt                                                                                                           >  1:224881‑M1/1‑95 (MQ=255)
                                                                                                                                                             gAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattcgcagttaatgttttgaaatcgttgctctgaattgctgataaatagacaaacgcatctt                                                                                                           <  2:224881‑M1/214‑120 (MQ=255)
                                                                                                                                                              aaTGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaaaacatctttgctgccatctttcggcgtgatgaaaacaaaacctttt                                                                     >  2:444097‑M1/1‑94 (MQ=255)
                                                                                                                                                              aaTGGCCACCGTCAGTAAATTTTATACCCGGGGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTATCCCACTCACtacatgaacagcggcaggtcctttaggtccgttcgcaataccaaattcatcttccggattctcaggtaatgttttgaaatcgttgctctgaattgctgagaaatatacaaacacatctttgctgccatctttcggcgtgatgaaaccaaatcctttt                                                                     <  2:543487‑M1/251‑158 (MQ=255)
                                                                                                                                                                    cACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcaggg                                                               <  2:9286‑M1/251‑164 (MQ=255)
                                                                                                                                                                    cACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcaggg                                                               <  1:561996‑M1/251‑164 (MQ=255)
                                                                                                                                                                          cAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaac                                                         >  1:603772‑M1/1‑82 (MQ=255)
                                                                                                                                                                                       aCCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacccttttcagggttaaaccatttcactaaac                                            >  1:432439‑M1/1‑69 (MQ=255)
                                                                                                                                                                                                aaCTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattt                                   <  1:504436‑M1/251‑192 (MQ=255)
                                                                                                                                                                                                 aCTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACgacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcatttt                                  <  2:178300‑M1/251‑193 (MQ=255)
                                                                                                                                                                                                    gCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgtt                               >  2:407385‑M1/1‑56 (MQ=255)
                                                                                                                                                                                                    gCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgt                                <  1:537438‑M1/250‑195 (MQ=255)
                                                                                                                                                                                                                  tGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattatt                  <  1:155542‑M1/250‑209 (MQ=255)
                                                                                                                                                                                                                  tGTCACAAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattat                   <  1:584702‑M1/249‑208 (MQ=255)
                                                                                                                                                                                                                     cacaAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattacct              <  2:305717‑M1/251‑213 (MQ=255)
                                                                                                                                                                                                                       caAAAAAGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttt            <  1:41696‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                        aaaaaaGTTGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattacctttt           <  1:334526‑M1/251‑216 (MQ=255)
                                                                                                                                                                                                                          aaaaGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttg          >  2:115773‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                             aggggtagACTCATGCAGTTAACCCACTCACtacatgatcagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgcttagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcagtttgttagacataattattacctgatgaaga      <  2:566224‑M1/248‑221 (MQ=255)
                                                                                                                                                                                                                             aGTGGTAGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaaga      <  1:142056‑M1/251‑221 (MQ=255)
                                                                                                                                                                                                                                   aGACTCATGCAGTTAACCCACTCACtacatgaacagcggcaggacctttaggtccgttctcaataccaaattcaacttcctgattctcagttaatgttttgaaatcgttgctctgaattgctgagaaatagacaaacacatctttgctgccatctttcggcgtgatgaaaccaaaacctttttcagggttaaaccatttcactaaaccagtcattttgttagacataattattaccttttgaagaaatt  >  1:524257‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       
AAGGCGATGACGGTAAAGCAGAAAAAACCTATCAAAATACCGTTGCATTGACGCCAGAAGATGTCAGCGAAGCCGTCTGGTGGGTGTCAACGCTGCCTGCTCACGTCAATATCAATACCCTGGAAATGATGCCGGTTACCCAAAGCTATGCCGGACTGAATGTCCACCGTCAGTAATTTTTATACCCGGCGTAACTGCCGGGTTATTGCTTGTCACAAAAAAGTGGTAGACTCATGCAGTTAACTCACTCACAAGCAAGAACGAATGACCGTCGAAACGCAACTTAATCCCACACAGCCTGTCAATCAGCAGATTTATCGTATTCTTCGTCGCGACATTGTCCATTGCCTGATTGCTCCAGGCACACCGTTGTCGGAAAAAGAAGTTTCTGTTCGTTTCAATGTGTCACGCCAGCCGGTTCGTGAAGCCTTTATTAAACTGGCGGAAAACGGCCTGATTCAAATTCGTCCGCAACG  >  NC_000913/1628088‑1628563

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