Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 209,608 T→C 7.4% intergenic (+28/‑71) 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,6080TC 7.4% 266.5 / 9.7 94intergenic (+28/‑71)accA/ldcCacetyl‑CoA carboxylase, carboxytransferase, alpha subunit/lysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (44/43);  new base C (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 = 9.85e-01

AAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTG  >  NC_000913/209459‑209756
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                  <  2:878634/150‑1 (MQ=255)
                    aaaGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCg                                                                                                                                  <  1:1342437/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                                                                                                                   <  1:1003704/148‑1 (MQ=255)
                                     ggATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAg                                                                                                                   >  2:1003704/1‑148 (MQ=255)
                                         cTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATg                                                                                                             <  2:302888/150‑1 (MQ=255)
                                           ggCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAt                                                                                                           <  2:84405/150‑1 (MQ=255)
                                           ggCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAt                                                                                                           <  2:1360322/150‑1 (MQ=255)
                                            gCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAtt                                                                                                          <  1:1459703/150‑1 (MQ=255)
                                            gCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGAtt                                                                                                          <  2:1683677/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                                                                                   >  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                                                                                                  >  2:619735/1‑127 (MQ=255)
                                                                           tAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTaaggaagg                                                                                                  <  1:619735/127‑1 (MQ=255)
                                                                           tAAAAAATCGTCGTTATCAGCGCCGGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGTGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCAATGATTAAGGAAGGATTTTCCAGGAGGAACACATGaa                                                                           >  1:1492756/1‑150 (MQ=255)
                                                                            aaaaaaTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCa                                                                                                                    <  1:937984/108‑1 (MQ=255)
                                                                            aaaaaaTCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCa                                                                                                                    >  2:937984/1‑108 (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                                 <  2:876622/150‑1 (MQ=255)
                                                                                                                     cgTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGa                                 >  1:1041469/1‑150 (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)
                                                                                                                                                    cTTCGGTGGCCCTTTTTTATCGCCACGGTTTGTGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCtg  >  2:146678/1‑150 (MQ=255)
                                                                                                                                                    cTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCtg  <  2:418321/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGCGATGGCGGCATCGTTGAAAGCGCAACTGCTGGCGGATCTGGCCGATCTCGACGTGTTAAGCACTGAAGATTTAAAAAATCGTCGTTATCAGCGCCTGATGAGCTACGGTTACGCGTAATTCGCAAAAGTTCTGAAAAAGGGTCACTTCGGTGGCCCTTTTTTATCGCCACGGTTTGAGCAGGCTATGATTAAGGAAGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTG  >  NC_000913/209459‑209756

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: