Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,739,066 A→T 7.2% intergenic (‑87/+70) proY ← / ← hutH proline (histidine) APC transporter/histidine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,739,0660AT 7.2% 632.7 / 2.0 240intergenic (‑87/+70)proY/hutHproline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (101/105);  new base T (3/18);  total (117/123)
Fisher's exact test for biased strand distribution p-value = 2.32e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

GAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATAT  >  NC_002947/5738917‑5739215
                                                                                                                                                     |                                                                                                                                                      
gAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGa                                                                                                                                                        >  2:674911/1‑150 (MQ=255)
gAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGa                                                                                                                                                        <  1:556589/150‑1 (MQ=255)
 aGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                                                                       <  1:5502078/150‑1 (MQ=255)
 aGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                                                                       <  1:3619368/150‑1 (MQ=255)
   ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATgg                                                                                                                                                     <  1:5386419/150‑1 (MQ=255)
   ggCCATCAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTCTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACGGGGGATgg                                                                                                                                                     >  2:5988739/1‑150 (MQ=255)
     ccATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGGGAGATTCAGGGGGCCGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgtg                                                                                                                                                  >  2:1813810/1‑150 (MQ=255)
       aTGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtt                                                                                                                                                <  1:3586674/150‑1 (MQ=255)
         gAAACGGTTGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtttt                                                                                                                                              <  1:595264/150‑1 (MQ=255)
            aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                                                           >  2:5588765/1‑150 (MQ=255)
            aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                                                           <  2:378649/150‑1 (MQ=255)
            aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCGCCCACCGGGGCTGGTGTTTTGcc                                                                                                                                           >  1:6630538/1‑149 (MQ=255)
            aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCCTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                                                           <  2:4930734/150‑1 (MQ=255)
             cGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCAc                                                                                                                                          <  1:1280900/150‑1 (MQ=255)
                aTGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTa                                                                                                                                       <  1:869847/150‑1 (MQ=255)
                     gCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCt                                                                                                                                  <  2:4307361/150‑1 (MQ=255)
                       gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                                >  2:3799896/1‑150 (MQ=255)
                       gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                                <  2:5451909/150‑1 (MQ=255)
                        ggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgtg                                                                                                                               <  1:6706037/150‑1 (MQ=255)
                         gCACTTAGCCCGCGCTTGAGACCTTGAGGTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTgg                                                                                                                              <  1:1476891/150‑1 (MQ=255)
                           aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTTAATTGTTTTTGTGGTCGTGAGATTCTGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                            <  1:3558780/150‑1 (MQ=255)
                           aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                            <  2:1526952/150‑1 (MQ=255)
                           aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                                                            <  1:3552033/150‑1 (MQ=255)
                             ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGc                                                                                                                          <  1:3854041/150‑1 (MQ=255)
                             ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCGCCCACACGGGGTGGTGTTTTTCACCGACGGGAGTGGGAGc                                                                                                                          >  1:2398635/1‑150 (MQ=255)
                              tcgCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCg                                                                                                                         <  2:6016864/148‑1 (MQ=255)
                               aGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                                                                        <  2:583251/150‑1 (MQ=255)
                               aGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                                                                        <  1:1711679/150‑1 (MQ=255)
                               aGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCGGTGGGAGCgg                                                                                                                        >  1:6166280/1‑150 (MQ=255)
                                 cccGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGgc                                                                                                                      >  2:506057/1‑150 (MQ=255)
                                   cgcgCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCGGTGGGTGCGGgcgc                                                                                                                    >  1:2712635/1‑150 (MQ=255)
                                       gttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTGGTGTGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGGGGGCGCGccc                                                                                                                <  1:5036615/149‑1 (MQ=255)
                                          gagaCCTTGAGCTTGTTGCATGTCACGTCCTTGATTTGTTTTTGTGGTGGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCcgcg                                                                                                             <  2:4035682/150‑1 (MQ=255)
                                           agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTTCACCTACTTCTGTGGGAGCGGGCGCGCCcgcca                                                                                                            >  2:2025000/1‑148 (MQ=255)
                                                  gAGCTTGTTGCATGTCCCGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCAGCTCCCACAGGGTATGTTTTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGgccg                                                                                                     <  1:4722092/150‑1 (MQ=255)
                                                    gCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCCCCCCCCCACCCAAACGGGATGGTGTTTTTGACCAACATCGGTGGGAGCGGGGGCGCCCCCCAACGgccccc                                                                                                   >  1:2962671/1‑150 (MQ=255)
                                                    gCTTGTTGCATGTCACGTCCTTGAATTGTGTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGGACGCCCCCACCCCCCGGGGCTGGTGGTTTTGCCCTACGTCGGGGGGTGCGGGCGCGCCCGCGCAGGgccgcc                                                                                                   >  1:2962688/1‑150 (MQ=255)
                                                             aTGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGccc                                                                                          >  2:2246701/1‑150 (MQ=255)
                                                                tCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAg                                                                                       >  2:5886133/1‑150 (MQ=255)
                                                                    gTCCTTGAATTGTTTTTGGGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCCGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGc                                                                                   <  1:6771439/150‑1 (MQ=255)
                                                                        ttGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCACCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTa                                                                               >  2:1196045/1‑150 (MQ=255)
                                                                          taaTTTTTTTTGTGGTCGTGTGATTCTGGGGGCGGCTGCGCAGCCCCTTCGGGGGCACGCCCGCTCCCACTGGGTTTGGTGTTTTGCACCTACTTCTGTGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                             <  1:1034943/149‑1 (MQ=255)
                                                                          gatttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                             <  2:6589503/147‑1 (MQ=255)
                                                                          gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                             <  2:1055329/150‑1 (MQ=255)
                                                                          gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCCCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                                                             <  2:1864213/150‑1 (MQ=255)
                                                                            aTTGTTTTTGTGGTCGTGAGTTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAgg                                                                           <  2:3585489/150‑1 (MQ=255)
                                                                             ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGc                                                                          <  2:5263186/150‑1 (MQ=255)
                                                                             ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGc                                                                          <  2:6237826/150‑1 (MQ=255)
                                                                             ttGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCACCCACAGGGGATTGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCCCCCGCGAAGGGCCCCCCAGCGGCCCCAGGGGCGTTAAAGGc                                                                          >  1:3403510/1‑150 (MQ=255)
                                                                               gTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCCGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCtt                                                                        <  2:2056046/150‑1 (MQ=255)
                                                                               gTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCtt                                                                        <  2:6823006/150‑1 (MQ=255)
                                                                                tttttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGGTGGTGGTTTGCCCCTACTTCTGTGGGAGGGGGGGCGCCCCCGATGGGCCGCCCAGCGGCCCCAGGGGCGTTAAAGGCGTg                                                                       >  1:1630046/1‑150 (MQ=38)
                                                                                tttttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGGGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                                                       <  2:2722845/150‑1 (MQ=255)
                                                                                tttttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCACCCACACGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCCGCGGCCCCCGGGGCGGTACACGCTTg                                                                       >  2:6332619/1‑150 (MQ=255)
                                                                                 ttttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCCCCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTgg                                                                      <  2:2639338/150‑1 (MQ=255)
                                                                                  ttggtgGTCGTGAGCTTCCGGGGGCGGCTGCGCAGCCCCTTCGCGGGCCCGCCCGCTCCCCCCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGc                                                                     <  2:5927414/147‑1 (MQ=255)
                                                                                   ttGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                                    <  2:4791575/150‑1 (MQ=255)
                                                                                   tggtgGTCGGGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCCGGGGTTGGTGTTTTGCCCCTACTTCGGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                                    <  1:3889394/148‑1 (MQ=255)
                                                                                   ggggggtggtGAGTTTCTGGGGGGGGGGGCGCCGCCCGTGCGGGGGCCCGCCCGCGCCCCCGGGGGTTGGTGTTTGCCCCCTCCTGGGGGGGGCGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCa                                                                    <  1:6428865/146‑1 (MQ=255)
                                                                                      tgGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCCCTCCCACAGGGGGTGGTGTTTTTCACCAAATTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCCCCCCGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTa                                                                 >  1:2689196/1‑150 (MQ=255)
                                                                                             gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgc                                                          >  2:2224181/1‑150 (MQ=255)
                                                                                               ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgccg                                                        >  2:3840081/8‑150 (MQ=39)
                                                                                               ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgccg                                                        >  2:6075968/8‑150 (MQ=39)
                                                                                               ccgcgaaggggCTGCGCAGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCGACGTCGGTGGGAGCGGGCGCGCCCGCGCAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCCGTACACccgccg                                                        >  1:5029865/8‑150 (MQ=255)
                                                                                                aTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCACCGGCCCCAGGGGCGGTACAGGCTTGGCAGTACACCCGCCgg                                                       >  2:1442999/1‑150 (MQ=255)
                                                                                                 ttCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCGCCCACAGGGGATGGTGTTTTGCACCAACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCCGCGGCCCCACGGGCGTTACAGGCTTGGCAGTACACCCGCCGgc                                                      >  1:4988314/1‑150 (MQ=255)
                                                                                                      ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATTGTGTTTTGCCCCTACGTGTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCCGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGgcagcac                                                 >  2:3494753/1‑149 (MQ=39)
                                                                                                      ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCCCGAAGGGCCGCCCAGCGGCCCCACGGGCGTTACAGGCTTGGCAGGACACCCCCCGGCAGCaa                                                 >  2:3067563/1‑150 (MQ=255)
                                                                                                         gctgctGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGGTGGTGGTTTGCACCGACTTCGGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCagccagcc                                              >  1:5574801/1‑150 (MQ=255)
                                                                                                            gctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCGGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGt                                           >  2:4514328/1‑150 (MQ=255)
                                                                                                              tgcGCAGCCCCTTCGCGGGCACGCCCGCGCCCACTGGGGATGGTGTTTTGCACCTACGTCGGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCCGTACACCCGCCGGCCGCAAGCCGGTcc                                         >  2:3127830/1‑149 (MQ=255)
                                                                                                                        cTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACGTGTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCCCCCGCAGCAAGCCGGTCAGACGAcccat                               >  1:2496961/1‑148 (MQ=255)
                                                                                                                         ttCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTg                              >  2:3218709/1‑150 (MQ=255)
                                                                                                                         tgcgcgGGCACGCCCGCTCCCACAGGGGTTGGTGTTTGGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTg                              <  2:655705/148‑1 (MQ=255)
                                                                                                                          tCGCGCGCACGCCCGCTCCCACAGGGGATGGTGTTTGGCACCTACTTCGGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTgg                             <  2:1478417/150‑1 (MQ=255)
                                                                                                                              gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTCCTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACCGGCTTGGCCGTCCACCCGCCGGCAGCACGCCGGTCAGACTACCCTGGCCacc                        >  2:5376554/1‑148 (MQ=39)
                                                                                                                              gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                                                                >  2:4078286/1‑47 (MQ=255)
                                                                                                                              gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAgg                                                                           <  2:3273032/100‑1 (MQ=255)
                                                                                                                              gggCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAgg                                                                           >  1:3273032/1‑100 (MQ=255)
                                                                                                                               ggCACGCCCGCTCCCACGGGGGATGGTGTTTGGCACCTACTTCTGTGGGCGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAAc                        <  2:1357238/150‑1 (MQ=255)
                                                                                                                                gCACGCCCGCTCCCACGGGGGTTGGTGTTTTGCACCTACTGCGGGGGGAGGGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACa                       <  1:960900/150‑1 (MQ=255)
                                                                                                                                gCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACa                       <  2:4772362/150‑1 (MQ=255)
                                                                                                                                  aCGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAgg                     <  2:5693788/150‑1 (MQ=255)
                                                                                                                                   cGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGt                    >  2:2465038/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCCACAGGGGATGGTGTTTTGCACCTACGTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCAGCCGGCAGCAAGCCGGTCAGACAACCCGTGGCCAACAGGTcc                  >  1:5989927/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCa                 >  1:5641356/1‑150 (MQ=255)
                                                                                                                                          tCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGc             >  1:4544726/1‑150 (MQ=255)
                                                                                                                                            ccACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCt           >  1:533922/1‑150 (MQ=255)
                                                                                                                                             cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt          >  1:1344526/1‑150 (MQ=255)
                                                                                                                                             cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt          >  2:3359263/1‑150 (MQ=255)
                                                                                                                                             cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt          >  1:1965750/1‑150 (MQ=255)
                                                                                                                                             cacaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtt          >  2:86661/1‑150 (MQ=255)
                                                                                                                                              acaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCATGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCCGACTACCCTTGGCCAACAGGTCCACGGCCttt         >  1:3142187/1‑150 (MQ=255)
                                                                                                                                              acaGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCttt         >  1:4417516/1‑150 (MQ=255)
                                                                                                                                               caGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCtttt        <  2:6709263/150‑1 (MQ=255)
                                                                                                                                                aGGGGTTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGACAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTc       <  2:5096422/150‑1 (MQ=255)
                                                                                                                                                aGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTc       >  2:5765952/1‑150 (MQ=255)
                                                                                                                                                aGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTc       >  2:3920024/1‑150 (MQ=255)
                                                                                                                                                aGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTc       >  1:4823171/1‑150 (MQ=255)
                                                                                                                                                     aTGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGatat  >  2:5491101/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGACTACCCTTGGCCAACAGGTCCACGGCCTTTTCGATAT  >  NC_002947/5738917‑5739215

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