Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,734,075 T→C 5.6% intergenic (‑45/+112) 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,0750TC 5.6% 431.6 / 9.5 141intergenic (‑45/+112)proY/hutHproline (histidine) APC transporter/histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (65/68);  new base C (3/5);  total (68/73)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 148 total aligned reads displayed.

GAGCCGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGG  >  NC_002947/5733933‑5734224
                                                                                                                                              |                                                                                                                                                     
   ccGTAGAAAAGGCTGGTGCCGATGGCGGAGCTGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGc                                                                                                                                             <  2:926703/150‑1 (MQ=255)
   ccGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGc                                                                                                                                             <  1:1394478/150‑1 (MQ=255)
             ggCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCTTTTTAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                   <  1:286984‑M1/150‑1 (MQ=255)
             ggCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                   <  1:2278626‑M1/150‑1 (MQ=255)
             ggCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGCGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                   <  2:164908‑M1/150‑1 (MQ=255)
                cGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCa                                                                                                                                >  1:317990/1‑150 (MQ=255)
                   tGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGc                                                                                                                             <  1:1008546/150‑1 (MQ=255)
                       gATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGc                                                                                                                         >  1:2306378/1‑150 (MQ=255)
                       gATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCCTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGc                                                                                                                         <  1:1692110/150‑1 (MQ=255)
                        aTGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCt                                                                                                                        >  1:1720054/1‑150 (MQ=255)
                          ggCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTAGCGGGCgtgcccgcgaa                                                                                                                      <  1:514485‑M1/150‑12 (MQ=255)
                                gCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                                <  2:275973/150‑1 (MQ=255)
                                 ccGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAggg                                                                                                               >  1:1921311/1‑150 (MQ=255)
                                    aGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGAt                                                                                                            <  2:1860555/150‑1 (MQ=255)
                                      ggCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATgg                                                                                                          >  2:41498/1‑150 (MQ=255)
                                        ccATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgtg                                                                                                        <  1:2352981/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)
                                           tGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCGTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCGCGCCCGCTCCCACCGGGGATGGTGttt                                                                                                     >  1:224593/1‑150 (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)
                                                  gATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCt                                                                                              >  1:366957/1‑150 (MQ=255)
                                                    tgtgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTAc                                                                                            >  2:2234792/1‑150 (MQ=255)
                                                      tgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatccc                                                                                          <  2:300147‑M1/150‑40 (MQ=255)
                                                      tgGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACtt                                                                                          >  1:1513047/1‑150 (MQ=255)
                                                          gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGTATGGTGTTTTGCACCTACTTCtgt                                                                                      <  1:573261/150‑1 (MQ=255)
                                                          gggCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCtgt                                                                                      >  2:1456824/1‑150 (MQ=255)
                                                            gCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacaggggatggtgcaaaacaccatcccctgtgg                                                                                    <  2:1732371‑M1/150‑46 (MQ=255)
                                                              aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtggga                                                                                  <  2:547905‑M1/150‑48 (MQ=255)
                                                              aCTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGa                                                                                  <  2:1513047/150‑1 (MQ=255)
                                                                ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGGATGGTGTTTTGCACCTACTTCTGTGGGatc                                                                                >  1:268217/1‑148 (MQ=255)
                                                                ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGc                                                                                >  1:699530/1‑150 (MQ=255)
                                                                ttAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCCGCTGGGCGGCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgggagcg                                                                               <  1:1821205‑M1/150‑51 (MQ=255)
                                                                      cgcgCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgcgc                                                                          >  2:550311/1‑150 (MQ=255)
                                                                           ttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                    <  2:212187‑M1/135‑46 (MQ=255)
                                                                           ttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                    >  1:212187‑M1/1‑90 (MQ=255)
                                                                           ttGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCgcccg                                                                     <  2:2196296‑M2/150‑6 (MQ=255)
                                                                              agaCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacagaagtaggtgcaaaacaccatcccctgtg                                                                                     >  2:397577‑M1/1‑87 (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)
                                                                                aCCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCAc                                                                                                <  2:158562/118‑1 (MQ=255)
                                                                                  cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCgcgcccgctcccacagaagtaggtgcaaaacaccatcccctgtgg                                                                                    <  1:1080238‑M1/128‑46 (MQ=255)
                                                                                  cTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGAAGCGGGAGCgcccgcgaagag                                                              <  2:1941117‑M2/150‑13 (MQ=255)
                                                                                   ttGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCgg                                                                              <  2:118588/133‑1 (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)
                                                                                                aTGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGccc                                                <  1:2196296/150‑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)
                                                                                                  gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGgc                                                                            >  2:1737954/1‑120 (MQ=255)
                                                                                                  gTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCa                                              >  1:837159/1‑150 (MQ=255)
                                                                                                   tCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAg                                             >  1:1033625/1‑150 (MQ=255)
                                                                                                       gTCCTTTAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGTCCCCAGGGGc                                         <  1:1299071/150‑1 (MQ=255)
                                                                                                       gTCCTTGAATTGTTTTTGTGGTCGGGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGc                                         >  2:2124008/1‑150 (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:2013335/1‑150 (MQ=255)
                                                                                                            tGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTAc                                    <  1:41498/150‑1 (MQ=255)
                                                                                                             gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                   <  1:2054011/150‑1 (MQ=255)
                                                                                                             gAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACa                                   <  1:2234792/150‑1 (MQ=255)
                                                                                                               aTTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGc                                                                                                  >  2:844458/1‑85 (MQ=255)
                                                                                                               aTTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGc                                                                                                  <  1:844458/85‑1 (MQ=255)
                                                                                                                  gTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCtt                              <  2:1928449/150‑1 (MQ=255)
                                                                                                                    ttttGTGGGCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTgg                            <  1:469771/150‑1 (MQ=255)
                                                                                                                        gtgGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGt                        <  2:2360274/150‑1 (MQ=255)
                                                                                                                           gTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTaca                     >  2:1098135/1‑150 (MQ=255)
                                                                                                                            tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCacgcccgctcccacaggagtaggtgcaaaaca                                                                                                 <  2:1956585‑M1/73‑33 (MQ=255)
                                                                                                                            tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGgt                                                                                                         <  2:295631/65‑1 (MQ=255)
                                                                                                                            tCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTacac                    >  1:1199364/1‑150 (MQ=255)
                                                                                                                             tgTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAcc                   <  2:2247955/149‑1 (MQ=255)
                                                                                                                              gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAccc                  >  1:2126838/1‑150 (MQ=255)
                                                                                                                              gTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACAccc                  >  2:1078424/1‑150 (MQ=255)
                                                                                                                               tGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGGGGTGGTGTTTTGCACCCACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccg                 >  1:139517/1‑150 (MQ=255)
                                                                                                                                gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgc                <  1:1187474/150‑1 (MQ=255)
                                                                                                                                gagaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGG                                                                                                                                   <  1:1182735‑M1/35‑1 (MQ=255)
                                                                                                                                 agaTTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACccgcc               >  1:1762241/1‑150 (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       >  1:1536041/1‑150 (MQ=255)
                                                                                                                                         ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa       >  1:785611/1‑150 (MQ=255)
                                                                                                                                         ggggCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCaa       >  2:788571/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)
                                                                                                                                              |                                                                                                                                                     
GAGCCGTAGAAAAGGCCGGTGCCGATGGCGGAGCCGAGGGCCATGAAACGGATGTGGCGGGCACTTAGCCCGCGCTTGAGACCTTGAGCTTGTTGCATGTCACGTCCTTGAATTGTTTTTGTGGTCGTGAGATTCAGGGGGCTGCTGCGCAGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGGATGGTGTTTTGCACCTACTTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCCCAGCGGCCCCAGGGGCGTTACAGGCTTGGCAGTACACCCGCCGGCAGCAAGCCGG  >  NC_002947/5733933‑5734224

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: