Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,417,032 C→T intergenic (+141/‑49) 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,0320CT97.9% 172.5 / NA 48intergenic (+141/‑49)pta/yfcCphosphate acetyltransferase/putative transporter YfcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/31);  minor base G (0/1);  total (16/32)
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.99e-01

GTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCGCAATCACTGAATCCAAACCAACAAGAAGATGGGCAATGCCCGATACGTTGGTGATTATCTTTTTTGTTGCTATTTTAACCAGCCTTGCCACCTGGGTAGTTCCGGTGGGGATGTTTGACAGTCAGGAAGTGCAGTATCAGGTTGATGGTCAAACAAAAACACGCAAAGTCGTAGATCCACACTCATTTCGC  >  NC_000913/2416791‑2417278
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
gTATGCGCAAGCCGGTTAACGACCTGTCCAGTGGCGCACTGGTGGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCa                                                                                                                                                                                                                                               <  2:19297/251‑1 (MQ=255)
  aTGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCAt                                                                                                                                                                                                                                              <  1:362292/250‑1 (MQ=255)
           ccGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACt                                                                                                                                                                                                                                    <  1:401661/251‑1 (MQ=255)
                tAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACtttgt                                                                                                                                                                                                                                <  1:69213/250‑1 (MQ=255)
                  aCGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTac                                                                                                                                                                                                                             >  2:365241‑M1/1‑249 (MQ=255)
                         gTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataac                                                                                                                                                                                                                        <  1:458650‑M1/249‑8 (MQ=255)
                           cccGTGGCGCACTGGTTGACGATAACGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccct                                                                                                                                                                                                                     <  2:375427‑M1/250‑11 (MQ=255)
                                            gACGATATCTTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCTTCAGCATCCGCAGCTTTGCGCTGCGTATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCACTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTCCCCACTTTGTTacatacccctcaagtgaaaatggctttt                                                                                                                                                                                                   <  1:297215‑M1/251‑29 (MQ=255)
                                                  atCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatg                                                                                                                                                                                              <  1:329879‑M1/250‑34 (MQ=255)
                                                                cgcgCTGACTGCGATTCAGTCTGCACAGCANCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCTGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgtt                                                                                                                                                                                <  2:257457‑M1/250‑48 (MQ=255)
                                                                cgccatgacttCGATTCAGTCTGCACAGCAGCAGTAATCTCGTAATCATCCGCAGCTTTGCGCTGCGGATATCTGAATCGGAAATAATCGCTATTTCCTGTTTTTTATTCTCTTAATGTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTGTGCATACTTCCCGCTTTGTTtcataaccctcaagtgaaaatggcttttttatgttgtgtattgtgtt                                                                                                                                                                                <  2:366685‑M1/245‑48 (MQ=255)
                                                                 gcgcTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgt                                                                                                                                                                              >  2:361024‑M1/1‑202 (MQ=255)
                                                                 gcactgaCTGCGATTCAGTCTGCACAGCAGCAGTAAGCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCACTGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagggaacatggcttttttctgttgtggtttgtgttgg                                                                                                                                                                              >  2:39358‑M1/4‑202 (MQ=255)
                                                                          gCGATTCAGTCTGCCCAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTTATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgc                                                                                                                                                                      <  2:200632‑M1/250‑58 (MQ=255)
                                                                          gCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgca                                                                                                                                                                     <  1:365241‑M1/251‑59 (MQ=255)
                                                                          gCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCACACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgggtgtttgc                                                                                                                                                                      >  1:292463‑M1/1‑193 (MQ=255)
                                                                               tCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggct                                                                                                                                                                 <  1:504317‑M1/250‑63 (MQ=255)
                                                                                      gCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttca                                                                                                                                                           <  1:426340‑M1/249‑69 (MQ=255)
                                                                                       cacaGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgc                                                                                                                                                        <  2:526316‑M1/251‑72 (MQ=255)
                                                                                               cagTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtctt                                                                                                                                                 >  2:331871‑M1/1‑172 (MQ=255)
                                                                                                      ctcGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatc                                                                                                                                         >  1:93269‑M1/1‑165 (MQ=255)
                                                                                                        cGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgcca                                                                                                                                           <  1:177195‑M1/247‑85 (MQ=255)
                                                                                                               atcCGCAGCTTAGCGCTGCGGATATCTGAACCGGAAATAATCACTCTATCCGGTTTTTTATTCGCTTAATTTGCATTAATCCTTTTTGATTATCTTGCTTAACTGCGCTGCATCAATGACTTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataacactcaagtgaaaatggcttttttatgttgtgtattgcgctgtgtgtttgcaggctctttcatgcaggtcttatgacatcttgtcgac                                                                                                                                 >  2:446821‑M1/1‑156 (MQ=255)
                                                                                                                 ccGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaac                                                                                                                              <  1:453366‑M1/251‑98 (MQ=255)
                                                                                                                 ccGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaac                                                                                                                              <  1:433778‑M1/251‑98 (MQ=255)
                                                                                                                 ccGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgac                                                                                                                                    <  1:26951‑M1/245‑92 (MQ=255)
                                                                                                                 ccGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgac                                                                                                                                    >  2:26951‑M1/1‑154 (MQ=255)
                                                                                                                    cAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagt                                                                                                                           <  1:67686‑M1/251‑101 (MQ=255)
                                                                                                                     aGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagta                                                                                                                          <  1:475529‑M1/251‑102 (MQ=255)
                                                                                                                              cTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATGCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATAAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcactcgataatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcagttcttntgccatcttgacgacaaccgtaacattccac                                                                                                                 >  2:270344‑M1/1‑141 (MQ=255)
                                                                                                                                gCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgt                                                                                                               <  1:223980‑M1/251‑113 (MQ=255)
                                                                                                                                 cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgtt                                                                                                              <  1:250994‑M1/251‑114 (MQ=255)
                                                                                                                                                    taatCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttactaagacgt                                                                                           >  2:339491‑M1/1‑119 (MQ=255)
                                                                                                                                                      atCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgttttgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttg                                                                                         >  2:53261‑M1/1‑117 (MQ=255)
                                                                                                                                                        cACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgc                                                                                        <  2:368866‑M1/250‑136 (MQ=255)
                                                                                                                                                            aTTTCCGGTTTTTTATTCTCTTAATTTGCACTAAGCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCGCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcat                                                                                      <  2:103849‑M1/248‑138 (MQ=255)
                                                                                                                                                                             ctctTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcatgaagttga                                                                     >  2:520875‑M1/1‑94 (MQ=255)
                                                                                                                                                                               ctTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaaga                                                                  <  1:532448‑M1/249‑158 (MQ=255)
                                                                                                                                                                                 tAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagat                                                                 <  1:152340‑M1/248‑159 (MQ=255)
                                                                                                                                                                                       gCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTa                                                                                                                                                                                                                              <  2:142036‑M1/85‑2 (MQ=255)
                                                                                                                                                                                                               gCTTAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcg                                 >  1:19402‑M1/1‑60 (MQ=255)
                                                                                                                                                                                                                 ttAACTGCGCTGCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgatt                              >  1:78048‑M1/1‑58 (MQ=255)
                                                                                                                                                                                                                            gCATCAATGAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcaga                    <  2:459038‑M1/250‑204 (MQ=255)
                                                                                                                                                                                                                                    gAATTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaccaaaatg           >  2:22857‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                                      aTTGCGCCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcaggtatagcggttgaatcgcagcgtttcgattgaatggcagaaacaaaatgtc         <  2:78048‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                                           gcCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtattgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgaccacagtaacattcaacgt                                                                                                               >  2:64552‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                           gcCATCTCACTTTGCATACTTCCCACTTTGTTacataaccctcaagtgaaaatggcttttttatgttgtgtaatgtgttgtgtgtttgcaggctctttcatgcaggtcttatgccatcttgacgacaacagtaacattcaacgt                                                                                                               <  1:64552‑M1/144‑113 (MQ=255)
                                                                                                                                                                                                                                             cATCTCACTTTGTATACTTCCCACTTTGGTacataactctcaatggaaaatggattttttgtgttgtgtattgtgttgtgtggttgcaggctctttcatggaggtcttatgccatcttgccgacaacagtaactttcaacgttaaatatgttaataagacgttgcattattgtcctgaagttgaagatagcagttatagcggttgaatcgcagagtgtcgattgaatggcagaaacaaaatgtcgaggatt  <  2:283230‑M1/251‑222 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
GTATGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCGCAATCACTGAATCCAAACCAACAAGAAGATGGGCAATGCCCGATACGTTGGTGATTATCTTTTTTGTTGCTATTTTAACCAGCCTTGCCACCTGGGTAGTTCCGGTGGGGATGTTTGACAGTCAGGAAGTGCAGTATCAGGTTGATGGTCAAACAAAAACACGCAAAGTCGTAGATCCACACTCATTTCGC  >  NC_000913/2416791‑2417278

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