Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,417,047 A→C intergenic (+156/‑34) pta → / → yfcC phosphate acetyltransferase/putative transporter YfcC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,417,0470AC98.4% 232.1 / NA 64intergenic (+156/‑34)pta/yfcCphosphate acetyltransferase/putative transporter YfcC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (32/31);  total (33/31)
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.69e-01

ACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCGCAATCACTGAATCCAAACCAACAAGAAGATGGGCAATGCCCGATACGTTGGTGATTATCTTTTTTGTTGCTATTTTAACCAGCCTTGCCACCTGGGTAGTTCCGGTGGGGATGTTTGACAGTCAGGAAGTGCAGTATCAGGTTGATGGTCAAACAAAAACACGCAAAGTCGTAGATCCACACTCATTTCGCAT  >  NC_000913/2416809‑2417280
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
aCGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTac                                                                                                                                                                                                                               >  2:187938‑M1/1‑249 (MQ=255)
       gTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataac                                                                                                                                                                                                                          <  1:547001‑M1/249‑8 (MQ=255)
        tCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccct                                                                                                                                                                                                                       >  1:52442‑M1/1‑241 (MQ=255)
         cccGTTGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCATTGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctc                                                                                                                                                                                                                      >  2:314426‑M1/1‑240 (MQ=255)
         cccGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccct                                                                                                                                                                                                                       <  2:179871‑M1/250‑11 (MQ=255)
         cccGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccct                                                                                                                                                                                                                       <  2:232084‑M1/250‑11 (MQ=255)
          ccGTGGCGCACTGGTTGACGATATCGTCTACNCCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctca                                                                                                                                                                                                                     <  2:317698‑M1/251‑13 (MQ=255)
           cGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGTTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCCCTTTGTTacataaccctcaa                                                                                                                                                                                                                    <  1:284574‑M1/251‑14 (MQ=255)
           cGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaa                                                                                                                                                                                                                    <  1:492295‑M1/251‑14 (MQ=255)
              ggCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagt                                                                                                                                                                                                                  <  1:529481‑M1/250‑16 (MQ=255)
                   aCTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaat                                                                                                                                                                                                            <  1:668950‑M1/251‑22 (MQ=255)
                    cTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatg                                                                                                                                                                                                           >  1:92078‑M1/1‑229 (MQ=255)
                      ggTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatg                                                                                                                                                                                                           <  1:115785‑M1/249‑23 (MQ=255)
                          gacgaaagcgtctacacCATCGCGCTGACTNCGATTCAGGCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttt                                                                                                                                                                                                      <  2:353757‑M1/242‑28 (MQ=255)
                           aCGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACCTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttt                                                                                                                                                                                                    >  1:266194‑M1/1‑222 (MQ=255)
                            cGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggctttttt                                                                                                                                                                                                   >  2:528866‑M1/1‑221 (MQ=255)
                             gATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggctttttta                                                                                                                                                                                                  <  1:388541‑M1/251‑32 (MQ=255)
                             gATATCGTCTACAACATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatgtctttttta                                                                                                                                                                                                  <  1:503983‑M1/251‑32 (MQ=255)
                              atatCGTCTACCCCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCCGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttat                                                                                                                                                                                                 >  2:410147‑M1/1‑219 (MQ=255)
                                         acCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattg                                                                                                                                                                                      >  1:215721‑M1/1‑208 (MQ=255)
                                              cgcgCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgtt                                                                                                                                                                                  <  1:528866‑M1/250‑48 (MQ=255)
                                                        gCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgc                                                                                                                                                                        >  1:608500‑M1/1‑193 (MQ=255)
                                                                  cTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcc                                                                                                                                                             >  2:659870‑M1/1‑183 (MQ=255)
                                                                    gCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcgtttttatgttgtgtattgtgttgtgtgtttgcaggctctttcat                                                                                                                                                            <  1:256123‑M1/250‑70 (MQ=255)
                                                                           agcagTAATCTCGGCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtct                                                                                                                                                    >  2:259103‑M1/1‑174 (MQ=255)
                                                                               gTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatg                                                                                                                                                >  1:258986‑M1/1‑170 (MQ=255)
                                                                                      cGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccagctt                                                                                                                                         <  2:223603‑M1/251‑89 (MQ=255)
                                                                                               ccGCAGCTTTGCGCTGCGGATATCTGAACCGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaac                                                                                                                                <  1:696421‑M1/250‑98 (MQ=255)
                                                                                                         gcgcTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacatt                                                                                                                       <  1:411341‑M1/250‑107 (MQ=255)
                                                                                                               cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgtt                                                                                                                <  1:510034‑M1/251‑114 (MQ=255)
                                                                                                                   tatCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgc                                                                                                                                                          <  1:250435‑M1/205‑72 (MQ=255)
                                                                                                                   tatCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgc                                                                                                                                                          >  2:250435‑M1/1‑134 (MQ=255)
                                                                                                                      cTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatg                                                                                                         >  1:288000‑M1/1‑131 (MQ=255)
                                                                                                                        gAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgtt                                                                                                       <  1:349481‑M1/251‑123 (MQ=255)
                                                                                                                          aCCGGAAATAATCACTATTTCCGGTTTTTTATTCTATTAATTTGCATTAATCCTTTCTTATTATCTTGCTTAACAGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTATGTTccataaccctcaagtgaaaatggcttttttatgttgtgtcttgtgttgtgtgtttgcagcctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgtta                                                                                                      >  2:636043‑M1/1‑127 (MQ=255)
                                                                                                                                      cACTATTTCCGGTTTTTTGTTCTCTTAATTAACATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacatcaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgaaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttactaagacgttgc                                                                                          >  1:559697‑M1/1‑115 (MQ=255)
                                                                                                                                          aTTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattat                                                                                     >  2:19006‑M1/1‑111 (MQ=255)
                                                                                                                                              ccGGTTTTTTATTCTCTTAATTTTCAGTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataacccgcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtc                                                                                 <  2:624409‑M1/251‑145 (MQ=255)
                                                                                                                                                  ttttttATTCCCTTACTTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGATTTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgcgtgcttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaactttcaacgctctatatgttaatacgaagttgcattattgtcctgc                                                                             >  2:578469‑M1/1‑103 (MQ=255)
                                                                                                                                                             ctTAATTTGCCTTAATACTTTCTGAATATCGTGCTTAACTCCGCTGCATCAATGAATAGCGCCATCTCACTTTGCATACTTCACCCTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgccaaccgtaacattcaacgttaaatatgttaataatacgttcaattattgtcctgaagtttaagtt                                                                   >  2:671780‑M1/1‑92 (MQ=255)
                                                                                                                                                              ttAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgggtgtttgcaggctctttcatgaaggtcttatgccatcttgacgacaacagtaacattcaacgttcaatatgttaataagacgttgcattattgtcctgaagttgacgat                                                                   >  1:359817‑M1/1‑91 (MQ=255)
                                                                                                                                                                aaattgcaTTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTCTGTTacataaccctcaagcgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagc                                                                >  1:53361‑M1/4‑89 (MQ=255)
                                                                                                                                                                aaTTTGAATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgacaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttaatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagc                                                                >  2:435833‑M1/1‑89 (MQ=255)
                                                                                                                                                                 aTTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcag                                                              <  1:19006‑M1/251‑164 (MQ=255)
                                                                                                                                                                     gCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtat                                                          >  1:43387‑M1/1‑84 (MQ=255)
                                                                                                                                                                       aTTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtat                                                          >  1:210068‑M1/1‑82 (MQ=255)
                                                                                                                                                                         tAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagc                                                       <  2:272501‑M1/250‑171 (MQ=255)
                                                                                                                                                                          aaTCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcgg                                                     >  1:130725‑M1/1‑79 (MQ=255)
                                                                                                                                                                           aTCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggt                                                    >  2:454698‑M1/1‑78 (MQ=255)
                                                                                                                                                                                ttCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaat                                               >  2:64652‑M1/1‑73 (MQ=255)
                                                                                                                                                                                 tCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatc                                              >  2:37058‑M1/1‑72 (MQ=255)
                                                                                                                                                                                    gATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacatatcccgcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacatcagtaacattcaacgttaaatatgttaataagacgttgcagtagtgtcctgaagttgaagatagcaggtatagcggttgaatcgca                                           <  2:488142‑M1/251‑183 (MQ=255)
                                                                                                                                                                                     aTTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCTCTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcagtattgtccggaagttgaagatagcaggtatagcggttgaatcgcag                                          <  2:565378‑M1/251‑184 (MQ=255)
                                                                                                                                                                                            tGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttc                                    <  1:304070‑M1/250‑190 (MQ=255)
                                                                                                                                                                                             gCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacagcagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgg                                  >  2:671116‑M1/1‑60 (MQ=255)
                                                                                                                                                                                                tAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgatt                                <  2:608500‑M1/250‑194 (MQ=255)
                                                                                                                                                                                                  aCTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattga                              >  1:46370‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                          gCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcaga                      <  2:52442‑M1/250‑204 (MQ=255)
                                                                                                                                                                                                            ctcAATGAATTGCCCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgccaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtactgaagttgaagatagcaggtatagcggttgaatcgcagcgtttagatttaatggcagaacc                   >  2:459679‑M1/2‑45 (MQ=255)
                                                                                                                                                                                                              cAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaaca                  >  1:715790‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                               aaTGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaacaaa                <  1:187938‑M1/251‑210 (MQ=255)
                                                                                                                                                                                                               aaTGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaacaaa                <  1:152867‑M1/251‑210 (MQ=255)
                                                                                                                                                                                                                   aaTTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgccggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaacaaaatgt            >  2:78208‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                                                cACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtggattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaacaaaatgtcgaggatttg  <  2:666010‑M1/248‑224 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
ACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCGCAATCACTGAATCCAAACCAACAAGAAGATGGGCAATGCCCGATACGTTGGTGATTATCTTTTTTGTTGCTATTTTAACCAGCCTTGCCACCTGGGTAGTTCCGGTGGGGATGTTTGACAGTCAGGAAGTGCAGTATCAGGTTGATGGTCAAACAAAAACACGCAAAGTCGTAGATCCACACTCATTTCGCAT  >  NC_000913/2416809‑2417280

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