Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,064 C→G 5.4% A263G (GCC→GGC)  mhpB → 3‑carboxyethylcatechol 2,3‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,0640CG 5.4% 315.1 / 6.8 92A263G (GCC→GGC) mhpB3‑carboxyethylcatechol 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (38/49);  new base G (3/2);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCA  >  NC_000913/370915‑371207
                                                                                                                                                     |                                                                                                                                               
gCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGc                                                                                                                                                 >  1:931679/1‑150 (MQ=255)
  aaCGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCa                                                                                                                                               <  1:515487/150‑1 (MQ=255)
    cGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCAtt                                                                                                                                             >  1:796481/1‑150 (MQ=255)
    cGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCAtt                                                                                                                                             <  1:1518723/150‑1 (MQ=255)
        tGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCg                                                                                                                                         <  1:795990/150‑1 (MQ=255)
         gATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCgg                                                                                                                                        >  1:45004/1‑150 (MQ=255)
          aTTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAGCGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGa                                                                                                                                       >  2:1265067/1‑150 (MQ=255)
           ttAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGaa                                                                                                                                      <  2:523259/150‑1 (MQ=255)
             aGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAg                                                                                                                                    >  1:662620/1‑150 (MQ=255)
               cgccgcTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGGCAATGGCGGAAAgct                                                                                                                                  <  2:610229/150‑3 (MQ=255)
               cgccgcTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTc                                                                                                                                  <  1:147456/150‑1 (MQ=255)
                gccgcTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCg                                                                                                                                 <  1:1511716/150‑1 (MQ=255)
                   gcTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGaca                                                                                                                              <  2:770867/150‑1 (MQ=255)
                    cTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGacac                                                                                                                             <  2:602513/150‑1 (MQ=255)
                    cTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGacac                                                                                                                             <  2:405789/150‑1 (MQ=255)
                     tGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGAcaca                                                                                                                            <  1:1165262/150‑1 (MQ=255)
                      gagaAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACAt                                                                                                                           >  2:1484015/1‑150 (MQ=255)
                       agaAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATg                                                                                                                          <  1:417482/150‑1 (MQ=255)
                        gaAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGa                                                                                                                         >  2:1037060/1‑150 (MQ=255)
                             ttgttgAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCa                                                                                                                    <  1:365511/150‑1 (MQ=255)
                              tgttgAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaa                                                                                                                   >  1:767918/1‑150 (MQ=255)
                               gttgAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaaa                                                                                                                  <  1:983773/150‑1 (MQ=255)
                               gttgAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaaa                                                                                                                  <  1:1312767/150‑1 (MQ=255)
                                ttgAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCaaaa                                                                                                                 <  1:1267917/150‑1 (MQ=255)
                                 tgAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGAATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAAc                                                                                                                <  1:871960/150‑1 (MQ=255)
                                    ggATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                             <  1:1593366/150‑1 (MQ=255)
                                    ggATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                             <  2:934507/150‑1 (MQ=255)
                                    ggATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                             >  1:855596/1‑150 (MQ=255)
                                    ggATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTg                                                                                                             >  2:536039/1‑150 (MQ=255)
                                         agagAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcg                                                                                                        >  1:986883/1‑150 (MQ=255)
                                         agagAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcg                                                                                                        >  1:381219/1‑150 (MQ=255)
                                         agagAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcg                                                                                                        >  2:1091081/1‑150 (MQ=255)
                                         agagAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcg                                                                                                        >  2:236902/1‑150 (MQ=255)
                                            gaACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccg                                                                                                     <  1:1368938/150‑1 (MQ=255)
                                             aaCGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTcgccgc                                                                                                    <  2:931679/150‑1 (MQ=255)
                                                gCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCttt                                                                                                 <  2:987954/150‑1 (MQ=255)
                                                            cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                                                                     >  1:651432/1‑150 (MQ=255)
                                                            cAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTc                                                                                     <  2:243240/150‑1 (MQ=255)
                                                                 cGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGt                                                                                >  2:218159/1‑150 (MQ=255)
                                                                   aTTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGttt                                                                              <  1:536039/150‑1 (MQ=255)
                                                                    tttGGGATAACCAGTTCATGACTTTGCTGGAGCATGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTg                                                                             <  1:408712/150‑1 (MQ=255)
                                                                       gggATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCa                                                                          >  2:503228/1‑150 (MQ=255)
                                                                            aaCCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGc                                                                    >  1:1592372/1‑150 (MQ=255)
                                                                             aCCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGc                                                                    >  2:722522/1‑150 (MQ=255)
                                                                              ccAGTTCATGACTTTGCTTGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCg                                                                   <  1:1664265/150‑1 (MQ=255)
                                                                              ccAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCg                                                                   <  1:928646/150‑1 (MQ=255)
                                                                                  ttCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGc                                                               <  1:1113741/150‑1 (MQ=255)
                                                                                   tCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCg                                                              <  1:1186299/150‑1 (MQ=255)
                                                                                      tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                                                           <  2:45004/150‑1 (MQ=255)
                                                                                      tGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAg                                                           <  2:128699/150‑1 (MQ=255)
                                                                                         cTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGc                                                        >  1:1544464/1‑150 (MQ=255)
                                                                                           ttGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt                                                      <  1:597556/150‑1 (MQ=255)
                                                                                           ttGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGt                                                      <  2:662620/150‑1 (MQ=255)
                                                                                            tGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtt                                                     <  2:986883/150‑1 (MQ=255)
                                                                                             gCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta                                                    <  1:236902/150‑1 (MQ=255)
                                                                                             gCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta                                                    <  1:1462248/150‑1 (MQ=255)
                                                                                             gCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGtta                                                    <  2:123499/150‑1 (MQ=255)
                                                                                               tGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatt                                                  >  1:387851/1‑150 (MQ=255)
                                                                                                ggAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGttatta                                                 >  1:752699/1‑150 (MQ=255)
                                                                                                   gCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCg                                              >  1:672601/1‑150 (MQ=255)
                                                                                                      gggACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGccc                                           <  1:269165/150‑1 (MQ=255)
                                                                                                        gACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCaa                                         >  2:1000051/1‑150 (MQ=255)
                                                                                                          cGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATc                                       <  2:752699/150‑1 (MQ=255)
                                                                                                            cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATccc                                     >  1:1568176/1‑150 (MQ=255)
                                                                                                            cATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATccc                                     >  2:915386/1‑150 (MQ=255)
                                                                                                             aTACAGGAACTGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCgg                                   <  2:1592372/150‑1 (MQ=255)
                                                                                                              tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCgg                                   <  1:1037060/150‑1 (MQ=255)
                                                                                                              tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCgg                                   >  1:506166/1‑150 (MQ=255)
                                                                                                              tACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCgg                                   <  1:784528/150‑1 (MQ=255)
                                                                                                               aCAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGa                                  >  2:519055/1‑150 (MQ=255)
                                                                                                                cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCg                                              <  2:1486225/137‑1 (MQ=255)
                                                                                                                cAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCg                                              >  1:1486225/1‑137 (MQ=255)
                                                                                                                   gAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTgg                              >  2:599026/1‑150 (MQ=255)
                                                                                                                      cTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGAtt                           <  1:1265067/150‑1 (MQ=255)
                                                                                                                        ggATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGc                         <  2:855596/150‑1 (MQ=255)
                                                                                                                             ccGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGAt                    <  1:1484015/150‑1 (MQ=255)
                                                                                                                                 cAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAACCCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTgg                >  1:573580/1‑150 (MQ=255)
                                                                                                                                  aGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGc               >  2:1580984/1‑150 (MQ=255)
                                                                                                                                   gTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCt              >  1:1368482/1‑150 (MQ=255)
                                                                                                                                     aaCGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCg            >  2:1346969/1‑150 (MQ=255)
                                                                                                                                      aCGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGt           >  1:180983/1‑150 (MQ=255)
                                                                                                                                       cGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGtt          <  1:519055/150‑1 (MQ=255)
                                                                                                                                        gAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTa         <  1:503228/150‑1 (MQ=255)
                                                                                                                                         aaGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTaa        >  1:817416/1‑150 (MQ=255)
                                                                                                                                            gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcg     >  1:1429227/3‑150 (MQ=255)
                                                                                                                                            gacTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcg     >  1:77292/3‑150 (MQ=255)
                                                                                                                                            agCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcg     <  2:506166/150‑1 (MQ=255)
                                                                                                                                             gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcgc    >  2:52613/1‑150 (MQ=255)
                                                                                                                                             gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcgc    >  1:1576109/1‑150 (MQ=255)
                                                                                                                                             gCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAgcgc    >  1:1039643/1‑150 (MQ=255)
                                                                                                                                              cTTTCCGGCAATGGCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGcc   >  2:1416552/1‑150 (MQ=255)
                                                                                                                                               tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCa  <  2:381219/150‑1 (MQ=255)
                                                                                                                                               tttCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCa  <  2:123040/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGCTGGAGCAGGGACGCATACAGGAACTGGATGCCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCA  >  NC_000913/370915‑371207

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