Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,067 T→A 6.2% I264N (ATT→AAT)  mhpB → 3‑carboxyethylcatechol 2,3‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,0670TA 6.2% 205.0 / 3.3 64I264N (ATT→AAT) mhpB3‑carboxyethylcatechol 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base A (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAG  >  NC_000913/370919‑371213
                                                                                                                                                    |                                                                                                                                                  
cGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCAtt                                                                                                                                                   >  2:810887/1‑150 (MQ=255)
cGGGTCATTAGCGCCGCTGAGAAGTTTCTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCAtt                                                                                                                                                   >  2:816005/1‑150 (MQ=255)
    tGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                                                               <  1:770228/150‑1 (MQ=255)
    tGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                                                               >  2:905866/1‑150 (MQ=255)
    tGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                                                               >  2:770228/1‑150 (MQ=255)
    tGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                                                               <  2:28063/150‑1 (MQ=255)
      aTTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGa                                                                                                                                             <  2:62834/150‑1 (MQ=255)
       ttAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGaa                                                                                                                                            <  1:1311613/150‑1 (MQ=255)
            gccgcTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgctc                                                                                                                                       <  1:843051/150‑4 (MQ=255)
               gcTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGaca                                                                                                                                    >  2:764892/1‑150 (MQ=255)
                cTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGTAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGacac                                                                                                                                   <  2:58340/150‑1 (MQ=255)
                      aGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGaaa                                                                                                                             >  1:1238495/1‑150 (MQ=255)
                                ggATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                                   >  1:15382/1‑150 (MQ=255)
                                    cAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTc                                                                                                               >  2:995690/1‑150 (MQ=255)
                                      gagaACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgc                                                                                                             >  2:1293001/1‑150 (MQ=255)
                                        gaACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccg                                                                                                           >  2:821571/1‑150 (MQ=255)
                                             cTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCtttt                                                                                                      <  2:1107979/150‑1 (MQ=255)
                                               gCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGc                                                                                                    >  1:455968/1‑150 (MQ=255)
                                                  tCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGc                                                                                                 <  1:407454/150‑1 (MQ=255)
                                                    cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                                   >  1:1285997/1‑130 (MQ=255)
                                                    cGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                                   <  2:1285997/130‑1 (MQ=255)
                                                        cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                                                                           <  1:979524/150‑1 (MQ=255)
                                                          aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGa                                                                                                                                             <  2:319088/98‑1 (MQ=255)
                                                          aCCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGa                                                                                                                                             >  1:319088/1‑98 (MQ=255)
                                                                  tGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTTGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                                                                 <  2:1106642/150‑1 (MQ=255)
                                                                  tGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGc                                                                                 <  1:1052095/150‑1 (MQ=255)
                                                                         aCCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGc                                                                          >  2:1158140/1‑150 (MQ=255)
                                                                                cATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGa                                                                   >  1:1355289/1‑150 (MQ=255)
                                                                                aaTGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGa                                                                   >  2:252412/2‑150 (MQ=255)
                                                                                            ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta                                                       >  1:1349214/1‑150 (MQ=255)
                                                                                               gCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCg                                                    <  1:158436/150‑1 (MQ=255)
                                                                                                 aGGGACGCATACAGGAAGTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGcc                                                  <  2:1150591/150‑1 (MQ=255)
                                                                                                  gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGccc                                                 <  1:995690/150‑1 (MQ=255)
                                                                                                     aCGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAAt                                              >  1:454612/1‑150 (MQ=255)
                                                                                                       gCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGGTAATAACGcccttcgc                                            <  1:781756/150‑6 (MQ=255)
                                                                                                        cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATccc                                           >  1:1379303/1‑150 (MQ=255)
                                                                                                        cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATccc                                           >  1:83382/1‑150 (MQ=255)
                                                                                                         aTACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCg                                          >  2:88219/1‑150 (MQ=255)
                                                                                                          tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCgg                                         <  1:562151/150‑1 (MQ=255)
                                                                                                           aCAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGa                                        >  2:1170068/1‑150 (MQ=255)
                                                                                                           aCAGGAACTGCATGCCGTCAGTAACGAAGAGCTTTCCCCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGa                                        >  2:1327592/1‑150 (MQ=255)
                                                                                                             aGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCCGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGt                                      <  1:905866/150‑1 (MQ=255)
                                                                                                              ggAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTg                                     >  2:589/1‑150 (MQ=255)
                                                                                                                   tGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTg                                >  2:544439/1‑150 (MQ=255)
                                                                                                                   tGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTg                                <  1:764892/150‑1 (MQ=255)
                                                                                                                    ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGc                               >  1:1344044/1‑150 (MQ=255)
                                                                                                                       tGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGGAATTTCTCCGTTTGGCAACTGCCGTAGCGAAGGGCGTTATTACCCCCCAATCCGGGAGTGGATTGCCgg                            >  1:1264084/1‑150 (MQ=255)
                                                                                                                       tGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCgg                            >  2:1114729/1‑150 (MQ=255)
                                                                                                                       tGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCgg                            >  1:361373/1‑150 (MQ=255)
                                                                                                                            tCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTg                       <  1:1293001/150‑1 (MQ=255)
                                                                                                                               gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCt                    >  2:491708/1‑150 (MQ=255)
                                                                                                                                     gtcgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTaa              >  1:1056943/6‑150 (MQ=255)
                                                                                                                                       gagCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgc            >  1:387342/1‑150 (MQ=255)
                                                                                                                                       cgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTa               <  1:1396434/144‑1 (MQ=255)
                                                                                                                                       cgacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTa               >  2:1396434/4‑147 (MQ=255)
                                                                                                                                        agCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcg           <  1:252412/150‑1 (MQ=255)
                                                                                                                                         gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta                                                          <  2:124733/102‑1 (MQ=255)
                                                                                                                                         gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta                                                          >  1:124733/1‑102 (MQ=255)
                                                                                                                                         gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcgc          >  1:649036/1‑150 (MQ=255)
                                                                                                                                         gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcgc          >  1:534935/1‑150 (MQ=255)
                                                                                                                                            ttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAg       <  1:1158140/150‑1 (MQ=255)
                                                                                                                                             tCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGa      >  2:506733/1‑150 (MQ=255)
                                                                                                                                               cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAAc    >  2:362937/1‑150 (MQ=255)
                                                                                                                                                gCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACa   <  2:330040/150‑1 (MQ=255)
                                                                                                                                                  cATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACag  <  2:366900/149‑1 (MQ=255)
                                                                                                                                                  cATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACag  >  1:366900/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAG  >  NC_000913/370919‑371213

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