Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 209,610 C→A 7.0% intergenic (+30/‑69) 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,6100CA 7.0% 281.8 / 9.1 94intergenic (+30/‑69)accA/ldcCacetyl‑CoA carboxylase, carboxytransferase, alpha subunit/lysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base A (2/5);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTG  >  NC_000913/209462‑209759
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                                           <  1:235944/145‑1 (MQ=255)
                gAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGCCACCGAAGTGACCCtttt                                                                                                                                           >  2:235944/1‑145 (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:84405/150‑1 (MQ=255)
                                        ggCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAt                                                                                                              <  2:1360322/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                                                                                                          >  1:1576210/1‑150 (MQ=255)
                                            gATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaag                                                                                                          >  2:23162/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:854070/150‑1 (MQ=255)
                                                      tgtTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGAtttt                                                                                                >  2:361088/1‑150 (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                                                                                      >  2:372180/1‑150 (MQ=255)
                                                                tGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGaa                                                                                      >  1:491282/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                                                                                                     <  1:619735/127‑1 (MQ=255)
                                                                        tAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaaggaagg                                                                                                     >  2:619735/1‑127 (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                                                                   >  2:383763/1‑150 (MQ=255)
                                                                                   cgtTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCa                                                                   <  1:383763/150‑1 (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:177816/1‑150 (MQ=255)
                                                                                                                  cgTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGa                                    <  2:876622/150‑1 (MQ=255)
                                                                                                                  cgTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGa                                    >  1:1041469/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                      >  2:725430/1‑150 (MQ=255)
                                                                                                                                gTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAg                      <  1:13957/150‑1 (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:1182996/1‑150 (MQ=255)
                                                                                                                                          aGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTc            <  2:1705064/150‑1 (MQ=255)
                                                                                                                                            ggTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCgg          >  2:17271/1‑150 (MQ=255)
                                                                                                                                             ttCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGCGACCGCAAGGCGTCTTTTATAAAGATGACCCCATCAAAGAACTGGAGACGgc         >  2:405020/2‑150 (MQ=255)
                                                                                                                                                 cTTCGGTGGCCCTTTTTTATCGCCACGGTTTGTGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCtg     >  2:146678/1‑150 (MQ=255)
                                                                                                                                                 cTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCtg     <  2:418321/150‑1 (MQ=255)
                                                                                                                                                    cGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCtggtg  <  1:422017/150‑1 (MQ=255)
                                                                                                                                                    cGGTGGCCCGTTTTTATCGGCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGCGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATACAGATGAGCCCATCAAAGAACTGGAGTCGGCGCtggtg  <  1:1029783/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTG  >  NC_000913/209462‑209759

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: