Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,734,083 Δ1 bp 5.4% intergenic (‑53/+104) 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,734,0830A. 5.4% 482.2 / 19.2 149intergenic (‑53/+104)proY/hutHproline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (69/71);  new base . (3/5);  total (73/76)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 149 total aligned reads displayed.

CCGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGAC  >  NC_002947/5733936‑5734230
                                                                                                                                                   |                                                                                                                                                   
ccGTAGAAAAGGCTGGTGCCGATGGCGGAGCTGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGc                                                                                                                                                   <  2:926703/150‑1 (MQ=255)
ccGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGc                                                                                                                                                   <  1:646282/150‑1 (MQ=255)
ccGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGc                                                                                                                                                   <  1:1394478/150‑1 (MQ=255)
          ggCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                         <  1:2278626‑M1/150‑1 (MQ=255)
             cGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCa                                                                                                                                      >  2:1946017/1‑150 (MQ=255)
                tGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGc                                                                                                                                   <  1:1008546/150‑1 (MQ=255)
                    gATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCCTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGc                                                                                                                               <  1:1692110/150‑1 (MQ=255)
                     aTGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCt                                                                                                                              >  1:1720054/1‑150 (MQ=255)
                       ggCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTAGCGGGCgtgcccgcgaa                                                                                                                            <  1:514485‑M1/150‑12 (MQ=255)
                         cGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTcccc                                                                                                                          >  1:275973/1‑149 (MQ=255)
                             gCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                                      <  2:275973/150‑1 (MQ=255)
                              ccGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAggg                                                                                                                     >  1:1921311/1‑150 (MQ=255)
                                   ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATgg                                                                                                                >  2:41498/1‑150 (MQ=255)
                                   ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATgg                                                                                                                <  2:1994723/150‑1 (MQ=255)
                                      cATGAAACGGATGTGGCGGGCACTTCGCCCGCGCGTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgtgt                                                                                                             >  1:666090/1‑150 (MQ=255)
                                      cATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgc                                                                                                             >  1:300147‑M1/1‑124 (MQ=255)
                                      cATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgtgt                                                                                                             >  1:216844/1‑150 (MQ=255)
                                       aTGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGtt                                                                                                            <  2:2090610/150‑1 (MQ=255)
                                       aTGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaa                                                                                                           >  2:730079‑M1/1‑122 (MQ=255)
                                        tGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGttt                                                                                                           <  1:283711/150‑1 (MQ=255)
                                         gAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaa                                                                                                          >  2:1907641‑M1/1‑120 (MQ=255)
                                          aaaCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTg                                                                                                         >  1:147317/1‑150 (MQ=255)
                                           aaCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgtgcaaaac                                                                                                        >  1:1732371/1‑144 (MQ=255)
                                            aCGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTATGCGTGAGATTCAGGCGGCTGTTGTCGGGGCTCTTCGCTGTCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                       >  1:2341978/1‑150 (MQ=255)
                                               gATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCGCGCCCGCTCCCACAGGGGATGGTGTTTTGCACCt                                                                                                    <  2:332074/150‑1 (MQ=255)
                                                 tgtgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTAc                                                                                                  <  1:2217035/150‑1 (MQ=255)
                                                 tgtgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTAc                                                                                                  >  2:2234792/1‑150 (MQ=255)
                                                   tgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatccc                                                                                                <  2:300147‑M1/150‑40 (MQ=255)
                                                   tgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACtt                                                                                                >  1:1513047/1‑150 (MQ=255)
                                                       gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                            <  1:2316933/150‑1 (MQ=255)
                                                       gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                            <  2:1921311/150‑1 (MQ=255)
                                                        ggCACTTAGCCCGCCCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatcccctgtg                                                                                           <  1:500897‑M1/150‑45 (MQ=255)
                                                           aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtggga                                                                                        <  2:547905‑M1/150‑48 (MQ=255)
                                                           aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                        <  2:1513047/150‑1 (MQ=255)
                                                           aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                        >  1:2247955/1‑150 (MQ=255)
                                                             ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgggagcg                                                                                     <  1:1821205‑M1/150‑51 (MQ=255)
                                                                   cgcgCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgcgc                                                                                >  2:550311/1‑150 (MQ=255)
                                                                        ttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                          >  1:212187‑M1/1‑90 (MQ=255)
                                                                        ttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                          <  2:212187‑M1/135‑46 (MQ=255)
                                                                           agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatcccctgtg                                                                                           >  2:397577‑M1/1‑87 (MQ=255)
                                                                           agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatcccctgtg                                                                                           <  1:397577‑M1/131‑45 (MQ=255)
                                                                           agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacac                                                                                                      <  1:2394037‑M1/119‑34 (MQ=255)
                                                                           agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacac                                                                                                      >  2:2394037‑M1/1‑86 (MQ=255)
                                                                            gaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCgcccgcgaa                                                                       >  2:216844‑M2/1‑141 (MQ=255)
                                                                             aCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCAc                                                                                                      >  1:158562/1‑118 (MQ=255)
                                                                               cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                          >  2:1080238‑M1/1‑83 (MQ=255)
                                                                               cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                          <  1:1080238‑M1/128‑46 (MQ=255)
                                                                               cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGc                                                                                                        >  2:214630/1‑114 (MQ=255)
                                                                               cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGc                                                                                                        <  1:214630/114‑1 (MQ=255)
                                                                               cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGAAGCGGGAGCgcccgcgaagag                                                                    <  2:1941117‑M2/150‑13 (MQ=255)
                                                                                ttGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGgc                                                                   >  2:1696489/1‑150 (MQ=255)
                                                                                        gttgCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGc                                                           <  1:1946017/150‑1 (MQ=255)
                                                                                         ttgCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                          >  1:737228/1‑86 (MQ=255)
                                                                                         ttgCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                          <  2:737228/86‑1 (MQ=255)
                                                                                         ttgCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCa                                                                                                       <  2:1711560/105‑1 (MQ=255)
                                                                                              tGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgggagcgg                                                                                    >  2:1974279‑M1/1‑68 (MQ=255)
                                                                                              tGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgggagcgg                                                                                    <  1:1974279‑M1/119‑52 (MQ=255)
                                                                                               gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgt                                                                                                                                     <  1:1641767‑M1/69‑3 (MQ=255)
                                                                                               gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgt                                                                                                                                     >  2:1641767‑M1/1‑67 (MQ=255)
                                                                                               gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                                  >  2:1737954/1‑120 (MQ=255)
                                                                                               gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCa                                                    >  1:837159/1‑150 (MQ=255)
                                                                                                tCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAg                                                   >  1:1033625/1‑150 (MQ=255)
                                                                                                tCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgggagcgggcgtgcccgcgaaggggctgcgcagcagccccct                                                  <  1:730079‑M1/150‑86 (MQ=255)
                                                                                                    gTCCTTTAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGTCCCCAGGGGc                                               <  1:1299071/150‑1 (MQ=255)
                                                                                                    gTCATTGATTTGTTGTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCCGGGGTTGTGTTTTGCACATACTTCTGTGGGAGCGGTCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGc                                               <  1:1375546/150‑1 (MQ=255)
                                                                                                        ttGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCTGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                            <  1:332263/101‑1 (MQ=255)
                                                                                                        ttGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                            >  2:332263/1‑101 (MQ=255)
                                                                                                         tGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTAc                                          <  1:41498/150‑1 (MQ=255)
                                                                                                         tGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTAc                                          >  1:2013335/1‑150 (MQ=255)
                                                                                                          gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                         <  1:2054011/150‑1 (MQ=255)
                                                                                                          gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                         <  1:2234792/150‑1 (MQ=255)
                                                                                                            aTTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGc                                                                                                        >  2:844458/1‑85 (MQ=255)
                                                                                                               gTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCtt                                    >  1:323493/1‑150 (MQ=255)
                                                                                                                     gtgGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGt                              <  2:2360274/150‑1 (MQ=255)
                                                                                                                         tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgt                                                                                                               <  2:295631/65‑1 (MQ=255)
                                                                                                                         tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgt                                                                                                               >  1:295631/1‑65 (MQ=255)
                                                                                                                         tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTacac                          >  1:1199364/1‑150 (MQ=255)
                                                                                                                          tgTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAcc                         <  2:2247955/149‑1 (MQ=255)
                                                                                                                           gggagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAccc                        <  2:1546797/148‑1 (MQ=255)
                                                                                                                           gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAccc                        >  2:1078424/1‑150 (MQ=255)
                                                                                                                            tGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGGTGGTGTTTTGCACCCACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccg                       >  1:139517/1‑150 (MQ=255)
                                                                                                                            tGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGGGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccg                       >  1:79457/1‑150 (MQ=255)
                                                                                                                             gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAg                                                                                       <  1:871835/85‑1 (MQ=255)
                                                                                                                             gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgc                      <  1:1187474/150‑1 (MQ=255)
                                                                                                                             gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                         <  1:1182735‑M1/35‑1 (MQ=255)
                                                                                                                             gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                         >  2:1182735‑M1/1‑35 (MQ=255)
                                                                                                                              agaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgcc                     >  1:1762241/1‑150 (MQ=255)
                                                                                                                               gaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAAccgccg                    <  2:1201689/150‑1 (MQ=255)
                                                                                                                                  tCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCATCCGGCTGCTTGCCAGCGGCCCCAGGGGCGTTACAGGCTTTGCAGTACACCCGCCGgca                 <  2:1191949/150‑1 (MQ=38)
                                                                                                                                    aGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCgcccgcgaagggccgcccagcggccccctgaatctcacgaccacaaaaacaattcaaggacgtga               <  1:730079‑M2/150‑66 (MQ=255)
                                                                                                                                     gggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCgcccgcgaagggccgcccagcggccccctgaatctcacgaccacaaaaacaattcaaggacgtgac              <  2:2288786‑M2/150‑67 (MQ=255)
                                                                                                                                      ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                   <  1:894775/128‑1 (MQ=255)
                                                                                                                                      ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                   >  2:894775/1‑128 (MQ=255)
                                                                                                                                      ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa             >  2:788571/1‑150 (MQ=255)
                                                                                                                                      ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa             >  1:1536041/1‑150 (MQ=255)
                                                                                                                                          ctgctgCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCg         >  1:243915/1‑150 (MQ=255)
                                                                                                                                           tgctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                   <  1:1145361/123‑1 (MQ=255)
                                                                                                                                           tgctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTg                                   >  2:1145361/1‑123 (MQ=255)
                                                                                                                                           tgctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACGCCCGCCGGCAGCAAGCCgg        >  2:1590739/1‑150 (MQ=255)
                                                                                                                                           tgctgcGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAGGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCgg        >  2:1748313/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGGTCAGAC  >  NC_002947/5733936‑5734230

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: