Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,935,901 C→G 6.3% intergenic (+153/‑109) PP_2568 → / → PP_2569 ring‑cleaving dioxygenase/MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,935,9010CG 6.3% 246.9 / 5.8 79intergenic (+153/‑109)PP_2568/PP_2569ring‑cleaving dioxygenase/MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base G (2/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGC  >  NC_002947/2935755‑2936046
                                                                                                                                                  |                                                                                                                                                 
ccccGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCg                                                                                                                                                >  1:1787618/1‑150 (MQ=255)
      cAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgt                                                                                                                                          >  2:1610016/1‑150 (MQ=255)
       aaaaTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCGAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGGCCGCCtcgga                                                                                                                                         <  2:1235478/150‑3 (MQ=255)
       aaaaTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                         >  2:1493771/1‑150 (MQ=255)
        aaaTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGGCCGCCtcggat                                                                                                                                        <  1:64669/150‑4 (MQ=255)
          aTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTg                                                                                                                                      <  1:105756/150‑1 (MQ=255)
          aTCGTCACCCAGCCCGCTATCGCCCGACACCCCCGTGCTGCGTTTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGCCGGGTGTATATCCAGCCAAGCCCGATTACGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTg                                                                                                                                      <  1:493968/150‑1 (MQ=255)
              tCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTc                                                                                                                                  >  2:273608/1‑150 (MQ=255)
                 gCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTcgcg                                                                                                                               >  1:126075/1‑150 (MQ=255)
                  caaaCCAGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCCCACGAAAGGCCCGACGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCCTCGTACTTCgcgc                                                                                                                              >  2:462583/3‑150 (MQ=255)
                    aaCCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGTTAACCTGCCAATCTCTATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCaa                                                                                                                            <  1:290585/150‑1 (MQ=255)
                       ccGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCa                                                                                                                         >  2:361732/1‑150 (MQ=255)
                           tATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagc                                                                                                                     >  1:1437209/1‑150 (MQ=255)
                              cGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagcaag                                                                                                                  <  1:2177000/150‑1 (MQ=255)
                               gCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagcaagc                                                                                                                 <  2:2307479/150‑1 (MQ=255)
                                cccGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAgcg                                                                                                                >  2:1432352/1‑150 (MQ=255)
                                   ggCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAATCAAAGCAAGCGCCa                                                                                                             <  1:506179/150‑1 (MQ=255)
                                   ggCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGccc                                                                                                             >  1:2140218/1‑149 (MQ=255)
                                       ccccGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGc                                                                                                         <  2:644079/150‑1 (MQ=255)
                                         ccGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCt                                                                                                       <  2:126075/150‑1 (MQ=255)
                                           ggTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTgg                                                                                                     <  1:1296046/150‑1 (MQ=255)
                                                 gcGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGa                                                                                               <  1:1610016/150‑1 (MQ=255)
                                                     gtcgtcCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGt                                                                                           >  1:227268/1‑150 (MQ=255)
                                                      tcgtcCAGGCCGCCCTCTCGTGCAGGGCCCGGCCGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTg                                                                                          >  1:312511/1‑150 (MQ=255)
                                                              gccgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGtt                                                                                  <  1:463090/150‑1 (MQ=255)
                                                              gccgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGtt                                                                                  <  1:1118707/150‑1 (MQ=255)
                                                              gccgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGtt                                                                                  <  2:588135/150‑1 (MQ=255)
                                                              gccgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATACCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGtt                                                                                  >  1:698429/1‑150 (MQ=255)
                                                                cgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGCGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTtgt                                                                                <  1:1352632/150‑1 (MQ=255)
                                                                cgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTtgt                                                                                >  1:1548513/1‑150 (MQ=255)
                                                                cgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTtgt                                                                                >  1:581581/1‑150 (MQ=255)
                                                                 gccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTtgtg                                                                               <  1:1779412/150‑1 (MQ=255)
                                                                   ccTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAg                                                                             <  2:227268/150‑1 (MQ=255)
                                                                     tctcGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGt                                                                           <  1:273608/150‑1 (MQ=255)
                                                                     tctcGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGt                                                                           >  1:1149669/1‑150 (MQ=255)
                                                                     tctcGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAACAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGt                                                                           >  1:514716/1‑150 (MQ=255)
                                                                             aGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCCGCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCGTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGGGAGGTGGTTGTGAGGTGCTCCGGc                                                                   >  1:555624/1‑150 (MQ=255)
                                                                                  ccGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGcc                                                                  >  2:1315961/1‑146 (MQ=255)
                                                                                  ccGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGcc                                                                  <  1:1315961/146‑1 (MQ=255)
                                                                                  ccGGCGGCTGGATATCCAGCCACTCCCGATGATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGCCGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTa                                                              <  2:1136539/150‑1 (MQ=255)
                                                                                  ccGGAGGGTGGTTATCCAGCCAATCCCGATTATGCTTCAGGAACCTTGCTGCGATCCGCGGCGGCCCGCCTAGTCGTGATTCGCGCAATCAAAGCAAGCGCGAGGGCCGGGGTGCGACGGTGAAGGGGTTGTGAGGTGCTCCGGCCTGTa                                                              <  2:514716/150‑1 (MQ=255)
                                                                                   cggcggGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTAt                                                             >  2:581283/1‑150 (MQ=255)
                                                                                     gcggGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGa                                                           <  2:1587584/150‑1 (MQ=255)
                                                                                         gTGGATATCCAGCTAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTgg                                                                                                     >  2:1850060/1‑104 (MQ=255)
                                                                                         gTGGATATCCAGCTAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTgg                                                                                                     <  1:1850060/104‑1 (MQ=255)
                                                                                          tGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGacaaca                                                      >  2:1987161/1‑150 (MQ=255)
                                                                                            gATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATg                                                    >  2:2098876/1‑150 (MQ=255)
                                                                                              tatCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGAc                                                  >  2:1157863/1‑150 (MQ=255)
                                                                                                tCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGa                                                   >  2:721570/1‑147 (MQ=255)
                                                                                                tCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACgg                                                >  2:1904686/1‑150 (MQ=255)
                                                                                                tCCAGCCAATCCCGACTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGa                                                   <  1:721570/147‑1 (MQ=255)
                                                                                                 ccAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGa                                               <  1:1157863/150‑1 (MQ=255)
                                                                                                   aGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                         <  2:97747/58‑1 (MQ=255)
                                                                                                   aGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCtcgtc                                                                                                                                         >  1:97747/1‑58 (MQ=255)
                                                                                                       aaTCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAc                                         >  1:834444/1‑150 (MQ=255)
                                                                                                        aTCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAcc                                        >  1:860058/1‑150 (MQ=255)
                                                                                                         tCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAccc                                       <  1:1904686/150‑1 (MQ=255)
                                                                                                          cccGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCa                                      <  2:834444/150‑1 (MQ=255)
                                                                                                           ccGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCAt                                     >  1:678988/1‑150 (MQ=255)
                                                                                                             gATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATgg                                   <  1:1493771/150‑1 (MQ=255)
                                                                                                             gATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATgg                                   <  2:581581/150‑1 (MQ=255)
                                                                                                               ttATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCaaa                                                                                                                       >  2:1247612/1‑64 (MQ=255)
                                                                                                               ttATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCaaa                                                                                                                       <  1:1247612/64‑1 (MQ=255)
                                                                                                               ttATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGACTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCg                                 <  2:312511/150‑1 (MQ=255)
                                                                                                                 aTGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGt                               <  1:581283/150‑1 (MQ=255)
                                                                                                                     ttCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGt                                                                                      <  2:1467400/91‑1 (MQ=255)
                                                                                                                     ttCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGt                                                                                      >  1:1467400/1‑91 (MQ=255)
                                                                                                                         tGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCGACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCGCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACa                       >  2:1231517/1‑150 (MQ=255)
                                                                                                                         tGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACa                       >  2:241389/1‑150 (MQ=255)
                                                                                                                         tGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACa                       <  2:1149669/150‑1 (MQ=255)
                                                                                                                            aCCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGca                    <  2:2321437/150‑1 (MQ=255)
                                                                                                                                   tgcGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCGAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCAc             >  1:2389704/1‑150 (MQ=255)
                                                                                                                                   tgcGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGGTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCTAc             >  1:560201/1‑150 (MQ=255)
                                                                                                                                    gcGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACg            <  1:2098876/150‑1 (MQ=255)
                                                                                                                                       aTCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACAAATGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGc         >  1:623918/1‑150 (MQ=255)
                                                                                                                                        gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGt                                                                           <  1:2071/81‑1 (MQ=255)
                                                                                                                                        gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGt                                                                           >  2:2071/3‑83 (MQ=255)
                                                                                                                                        gacGAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCa        >  1:361074/3‑150 (MQ=255)
                                                                                                                                          cGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAgc      >  2:698690/1‑150 (MQ=255)
                                                                                                                                              gcggcCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAgcgcgc  <  1:698690/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCCCGCCAAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGC  >  NC_002947/2935755‑2936046

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