Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 209,604 T→C 7.3% intergenic (+24/‑75) accA → / → ldcC acetyl‑CoA carboxylase, carboxytransferase, alpha subunit/lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913209,6040TC 7.3% 253.2 / 10.1 96intergenic (+24/‑75)accA/ldcCacetyl‑CoA carboxylase, carboxytransferase, alpha subunit/lysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (45/43);  new base C (2/5);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

CCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGC  >  NC_000913/209455‑209752
                                                                                                                                                     |                                                                                                                                                    
ccGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGt                                                                                                                                                      >  2:64976/1‑150 (MQ=255)
  ggAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGAACTCGACGTGTTAAGCACTGAAGATTTAAAAAAACGACCTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCa                                                                                                                                                    >  2:152279/1‑150 (MQ=255)
   gAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCAc                                                                                                                                                   <  1:491841/150‑1 (MQ=255)
    aaGCGACGGCGGCATCGTTGAAAGCGCATCTGCTCGCCGATCTGGCAGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTTTCATCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACt                                                                                                                                                  <  1:295108/150‑1 (MQ=255)
     aGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACtt                                                                                                                                                 >  2:952051/1‑150 (MQ=255)
     aGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACtt                                                                                                                                                 >  2:210375/1‑150 (MQ=255)
       cGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCAGAAAAAGGGTCACTTCg                                                                                                                                               >  1:1253149/1‑150 (MQ=255)
        gATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCgg                                                                                                                                              >  2:585678/1‑150 (MQ=255)
            gcggcATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGc                                                                                                                                          <  2:1412276/150‑1 (MQ=255)
            gcggcATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGc                                                                                                                                          >  2:1533679/1‑150 (MQ=255)
                  tCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGCCACCGAAGTGACCCtttt                                                                                                                                    <  2:40799/150‑1 (MQ=255)
                   cGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGCCACCGAAGTGACCCttttt                                                                                                                                   <  1:724431/150‑1 (MQ=255)
                     ttGATAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAATCACTGAAGATTTAAAAAATCGGCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGCCACCGAAGTGACCCtttttca                                                                                                                                 <  1:439271/150‑3 (MQ=255)
                     ttGAAAGCGCAACTGCTGGCGTATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTa                                                                                                                                 <  1:64976/150‑1 (MQ=255)
                       gAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTAta                                                                                                                               >  1:440156/1‑149 (MQ=255)
                       gAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGCCACCGAAGTGACCCtttt                                                                                                                                    >  2:235944/1‑145 (MQ=255)
                       gAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGCCACCGAAGTGACCCtttt                                                                                                                                    <  1:235944/145‑1 (MQ=255)
                        aaaGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCg                                                                                                                              <  1:1342437/150‑1 (MQ=255)
                        aaaGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCg                                                                                                                              <  2:878634/150‑1 (MQ=255)
                            cgcAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCAc                                                                                                                          >  1:1430823/1‑150 (MQ=255)
                            cgcAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGCTGGCCCTATTTTATCGCCAc                                                                                                                          >  2:856282/1‑150 (MQ=255)
                            cgcAACTGCTGGCGGATCAGGCCGATCTCGACGTGTAAAGCACTGAACATTTAAAAAATCGTCGTAATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTACTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCAc                                                                                                                          <  1:94211/150‑1 (MQ=255)
                             gcAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACg                                                                                                                         >  1:1377357/1‑150 (MQ=255)
                                  tgctgGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTg                                                                                                                    >  1:1096583/1‑150 (MQ=255)
                                  tgctgGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGGTTg                                                                                                                    <  2:500248/150‑1 (MQ=255)
                                        cGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAgg                                                                                                              >  1:876622/1‑150 (MQ=255)
                                         ggATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAg                                                                                                               >  2:1003704/1‑148 (MQ=255)
                                         ggATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAg                                                                                                               <  1:1003704/148‑1 (MQ=255)
                                             cTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATg                                                                                                         <  2:302888/150‑1 (MQ=255)
                                               ggCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAt                                                                                                       <  2:1360322/150‑1 (MQ=255)
                                               ggCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAt                                                                                                       <  2:84405/150‑1 (MQ=255)
                                                gCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAtt                                                                                                      <  2:1683677/150‑1 (MQ=255)
                                                gCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAtt                                                                                                      <  1:1459703/150‑1 (MQ=255)
                                                   gATCTCGACGTGTTAAGCTCTGAAGATTTAAAAAATCGTCGTTATCAGCTCCTGATGTGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGGGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaag                                                                                                   <  1:1304649/150‑1 (MQ=255)
                                                   gATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaag                                                                                                   >  2:23162/1‑150 (MQ=255)
                                                   gATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaag                                                                                                   >  1:1576210/1‑150 (MQ=255)
                                                     tctcGACGGGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaagga                                                                                                 <  1:1003872/150‑1 (MQ=255)
                                                      ctcGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaaggaa                                                                                                >  1:418321/1‑150 (MQ=255)
                                                       tcGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaaggaag                                                                                               <  1:1533679/150‑1 (MQ=255)
                                                             tgtTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAATGAtttt                                                                                         <  2:299984/150‑1 (MQ=255)
                                                             tgtTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGAtttt                                                                                         >  2:361088/1‑150 (MQ=255)
                                                             tgtTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGAtttt                                                                                         <  2:854070/150‑1 (MQ=255)
                                                             tgtTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGGGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGAtttt                                                                                         <  2:1309336/150‑1 (MQ=255)
                                                                  aGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCagg                                                                                    <  1:1063562/150‑1 (MQ=255)
                                                                  aGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCagc                                                                                    >  2:1314002/1‑149 (MQ=255)
                                                                      cTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAggagga                                                                                <  1:494947/150‑1 (MQ=255)
                                                                       tGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACTCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGaa                                                                               <  1:994159/150‑1 (MQ=255)
                                                                       tGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGaa                                                                               >  1:491282/1‑150 (MQ=255)
                                                                       tGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGaa                                                                               >  2:372180/1‑150 (MQ=255)
                                                                        gAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCATGAGGAac                                                                              >  2:285110/1‑150 (MQ=255)
                                                                        gAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAac                                                                              <  2:495847/150‑1 (MQ=255)
                                                                           gATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAAcaca                                                                           >  2:517271/1‑150 (MQ=255)
                                                                             tttAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATg                                                                         <  1:210375/150‑1 (MQ=255)
                                                                               tAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaaggaagg                                                                                              >  2:619735/1‑127 (MQ=255)
                                                                               tAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaaggaagg                                                                                              <  1:619735/127‑1 (MQ=255)
                                                                               tAAAAAATCGTCGTTATCAGCGCCGGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGTGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCAATGATTAAGGAAGGATTTTCCAGGAGGAACACATGaa                                                                       >  1:1492756/1‑150 (MQ=255)
                                                                                aaaaaaTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCa                                                                                                                >  2:937984/1‑108 (MQ=255)
                                                                                aaaaaaTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCa                                                                                                                <  1:937984/108‑1 (MQ=255)
                                                                                aaaaaaTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAAc                                                                      <  1:1708539/150‑1 (MQ=255)
                                                                                 aaaaaTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAAca                                                                     >  1:368831/1‑150 (MQ=255)
                                                                                     aTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAAcatcat                                                                 >  2:805785/1‑150 (MQ=255)
                                                                                     aTCGTCGTTATCAGAGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGCCGCCCTTTTTTATCGCCACGGATTCAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAAcatcat                                                                 <  1:896551/150‑1 (MQ=255)
                                                                                         tcgtTATCAGAGCCTGATGAGCTACGGTTACGCTTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGcc                                                             <  2:1279665/150‑1 (MQ=255)
                                                                                          cgtTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCa                                                            <  1:383763/150‑1 (MQ=255)
                                                                                          cgtTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCa                                                            >  2:383763/1‑150 (MQ=255)
                                                                                           gtTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCAt                                                           >  1:1583734/1‑150 (MQ=255)
                                                                                                cAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATgg                                                      >  1:335870/1‑150 (MQ=255)
                                                                                                 aGCGCCTGATGAGCTACGGTAACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATggg                                                     <  2:1171723/150‑1 (MQ=255)
                                                                                                     ccTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCCTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCg                                                 <  2:215048/150‑1 (MQ=255)
                                                                                                      cTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGc                                                >  2:1145610/1‑150 (MQ=255)
                                                                                                          tgaGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATgg                                            >  2:1143333/1‑150 (MQ=255)
                                                                                                            aGCTACGGTTACGCGTAATTCGCAAAAGTTCAGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCg                                          <  2:1253149/150‑1 (MQ=255)
                                                                                                                 cGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCttt                                     <  1:23162/150‑1 (MQ=255)
                                                                                                                  ggTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGACATTATGGGACCGCATGGCCTCtttt                                    >  2:768032/1‑150 (MQ=255)
                                                                                                                     tACGCGTATTTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTtata                                 <  2:801474/150‑1 (MQ=255)
                                                                                                                      aCGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATaa                                >  1:385709/1‑150 (MQ=255)
                                                                                                                         cgTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGa                             >  1:1041469/1‑150 (MQ=255)
                                                                                                                         cgTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGa                             <  2:876622/150‑1 (MQ=255)
                                                                                                                         cgTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGa                             >  1:177816/1‑150 (MQ=255)
                                                                                                                           tAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATg                           >  1:1147196/1‑150 (MQ=255)
                                                                                                                             aTTCGCAAACGATCGGAAAAAGGGTCCCTGCCGTGTCCCTTTTTTATCTCCACGGTTTGAGCAGGCAATGAGTAAGGAAGGATTTTCCTGGAGGAACACATTAACATCATTGCCATTATTGGACCTAATGGCGTCTTTTATAAAGATTAg                         <  2:1038568/150‑1 (MQ=255)
                                                                                                                               tCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGcc                       <  2:359496/150‑1 (MQ=255)
                                                                                                                                 gCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATTAACATCATTGCAATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCa                     >  2:156944/1‑150 (MQ=255)
                                                                                                                                   aaaaGTTCTGAAAAAGGGGCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATc                   <  1:1451432/150‑1 (MQ=255)
                                                                                                                                       gTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAg               <  1:13957/150‑1 (MQ=255)
                                                                                                                                       gTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAg               >  2:725430/1‑150 (MQ=255)
                                                                                                                                        ttCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTATGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGGCTTTTATAAAGATGAGCCCATCAAAGa              >  1:141382/1‑150 (MQ=255)
                                                                                                                                        ttCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGa              >  1:958393/1‑150 (MQ=255)
                                                                                                                                          ataaaaaaGGGCCACCGAAGTGACCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAg               >  2:129645/4‑147 (MQ=255)
                                                                                                                                          ataaaaaaGGGCCACCGAAGTGACCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAg               <  1:129645/144‑1 (MQ=255)
                                                                                                                                             aaaaaGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTgg         <  2:1281352/150‑1 (MQ=255)
                                                                                                                                              aaaaGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACTGCATTGCCATTATGGGACCTCATGGCGTCTTTTATAAATTTGAGACCTTCAATGAACgggg        >  2:840505/1‑149 (MQ=255)
                                                                                                                                               aaaGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAg       >  2:1581277/1‑150 (MQ=255)
                                                                                                                                               aaaGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAg       <  2:1147196/150‑1 (MQ=255)
                                                                                                                                                 aGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTc     <  2:1705064/150‑1 (MQ=255)
                                                                                                                                                 aGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTc     >  2:1182996/1‑150 (MQ=255)
                                                                                                                                                   ggTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCgg   >  2:17271/1‑150 (MQ=255)
                                                                                                                                                    ttCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGCGACCGCAAGGCGTCTTTTATAAAGATGACCCCATCAAAGAACTGGAGACGgc  >  2:405020/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGGAAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGC  >  NC_000913/209455‑209752

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