Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,510,778 C→G 7.6% P138R (CCG→CGG) ‡ urtA → urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,510,7780CG 7.6% 548.8 / 27.6 170P138R (CCG→CGG) ‡urtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (95/61);  new base G (7/6);  total (102/68)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 173 total aligned reads displayed.

AAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCTGCGC  >  NC_002947/5510629‑5510926
                                                                                                                                                     |                                                                                                                                                    
aaaGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGcc                                                                                                                                                      <  1:3451639/150‑1 (MQ=255)
           aaCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCt                                                                                                                                           >  1:3586969/1‑150 (MQ=255)
            aCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTa                                                                                                                                          >  1:4316618/1‑150 (MQ=255)
            aCTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTa                                                                                                                                          >  1:3159329/1‑150 (MQ=255)
              tgctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTaca                                                                                                                                        >  2:2090652/1‑150 (MQ=255)
                ctgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTAcaca                                                                                                                                      <  2:4777762/150‑1 (MQ=255)
                 tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                                                                                                                                     >  2:2560809/1‑150 (MQ=255)
                 tgACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAg                                                                                                                                     >  1:820805/1‑150 (MQ=255)
                     ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                                                                                                                                 >  1:2967894/1‑150 (MQ=255)
                     ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                                                                                                                                 <  2:3716216/150‑1 (MQ=255)
                     ccAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGgcgc                                                                                                                                 >  2:1802853/1‑150 (MQ=255)
                          aCAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGccc                                                                                                                            <  2:4386524/150‑1 (MQ=255)
                            aaaGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCg                                                                                                                          >  2:1088833/1‑150 (MQ=255)
                              aGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGaa                                                                                                                        >  1:2240546/1‑150 (MQ=255)
                               gTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAc                                                                                                                       <  2:3607468/150‑1 (MQ=255)
                                tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                                                                                                                      >  1:1300329/1‑150 (MQ=255)
                                tGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTAGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAAcc                                                                                                                      >  2:5121353/1‑150 (MQ=255)
                                 ggCCGTGGTTTTCGGCTGCTGGACATCCGTGTCGAGCAAGTCTGTTCTGAAGGTGTTAGAAGAACTAAACTGAATGCTGTTCTACAAGGTGCAGTATGAGGTCGAATAGATGTCGCCTAACGTGTTCTACAAAGGAGACGCCCCTAAACa                                                                                                                     >  2:2630234/1‑150 (MQ=255)
                                 ggCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCa                                                                                                                     <  2:1025797/150‑1 (MQ=255)
                                    agccgtgtTCGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCa                                                                                                                  <  1:1709933/146‑1 (MQ=255)
                                     gtggtgTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCaa                                                                                                                 <  2:3568303/150‑1 (MQ=255)
                                         tgtTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCa                                                                                                             >  2:4687177/1‑150 (MQ=255)
                                            tCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATcc                                                                                                          <  1:591744/150‑1 (MQ=255)
                                             cGGTTGCTGGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATccc                                                                                                         <  2:2283227/150‑1 (MQ=255)
                                               gttgctgGACGTCGGTGTCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATccctg                                                                                                       >  2:272999/3‑148 (MQ=255)
                                                   ctgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGt                                                                                                   <  1:4771909/150‑1 (MQ=255)
                                                    tgGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTg                                                                                                  >  2:1921545/1‑150 (MQ=255)
                                                      gACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                                                                                                >  1:5442699/1‑150 (MQ=255)
                                                      gACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGa                                                                                                >  1:2383576/1‑150 (MQ=255)
                                                         aTCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGTCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGACGCCCCGAACCAGCAAGCCATCCCGGCCTTGGAATa                                                                                             >  1:786003/1‑150 (MQ=255)
                                                           cggtgtCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGACTatc                                                                                           <  1:793379/148‑3 (MQ=255)
                                                           cggtgtCGCGCAAATCGGTTCTCCCGGTCTTCAAGGAACTCAACAATATTCTGTTCTACCCGGTGCAGTATGAAGGTGAGGAATCCGAGCGCAACGTTTTCTATACGGGCGCGGCACCGAACCAGCAGGCAATCCCTGCTGTTGACTatc                                                                                           >  1:1093673/3‑148 (MQ=255)
                                                           ccGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATAcc                                                                                           >  1:702789/1‑150 (MQ=255)
                                                            cGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCt                                                                                          <  1:1810689/150‑1 (MQ=255)
                                                                tCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgatg                                                                                      >  1:3567681/1‑150 (MQ=255)
                                                                tCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCtgatg                                                                                      >  1:1081113/1‑150 (MQ=255)
                                                                  gcgcAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGACTACCTGATGAg                                                                                    >  2:3932390/1‑150 (MQ=255)
                                                                   cgcAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGc                                                                                   <  1:4807774/150‑1 (MQ=255)
                                                                     cAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGa                                                                                 <  1:2208667/150‑1 (MQ=255)
                                                                     cAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGACTGCTGTTCTACCCGGTGCAGTATGAGTTCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGa                                                                                 <  1:4933492/150‑1 (MQ=255)
                                                                      aaGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAg                                                                                >  1:4122523/1‑150 (MQ=255)
                                                                           tGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAcgg                                                                           <  1:1165634/150‑1 (MQ=255)
                                                                                tGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGAcggatgtc                                                                      <  2:2383576/150‑6 (MQ=255)
                                                                                tGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCa                                                                      >  1:36115/1‑150 (MQ=255)
                                                                                  ccGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAgc                                                                    >  1:2497912/1‑150 (MQ=255)
                                                                                  ccGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAgc                                                                    <  2:3406583/150‑1 (MQ=255)
                                                                                        ttCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAg                                                              >  1:4994043/1‑150 (MQ=255)
                                                                                          cGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCa                                                                <  1:1581292/146‑1 (MQ=255)
                                                                                           gaagaaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAgcgc                                                           >  1:2804754/1‑150 (MQ=255)
                                                                                           gaagaaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAgcgc                                                           >  2:1001732/1‑150 (MQ=255)
                                                                                           gaagaaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAgcgc                                                           >  1:4314188/1‑150 (MQ=255)
                                                                                            aagaaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGct                                                          <  1:1677545/150‑1 (MQ=255)
                                                                                               aaCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGAGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCTTGGAATACCAGATGAGCGAGGACGGCGGCAGCGCCAAGCGcttct                                                       >  1:5303415/1‑150 (MQ=255)
                                                                                                aCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGcttctt                                                      >  2:4087941/1‑150 (MQ=255)
                                                                                                aCTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGcttctt                                                      >  2:787686/1‑150 (MQ=255)
                                                                                                 cTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCttcttc                                                     <  1:3557877/150‑1 (MQ=255)
                                                                                                 cTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCttcttc                                                     <  1:2313135/150‑1 (MQ=255)
                                                                                                 cTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCttcttc                                                     <  1:2106457/150‑1 (MQ=255)
                                                                                                  tCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCa                                                                      >  2:4805207/1‑132 (MQ=255)
                                                                                                  tCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTcc                                                    >  1:3821592/1‑150 (MQ=255)
                                                                                                    aaCGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAgcgc                                                           >  1:101661/1‑141 (MQ=255)
                                                                                                     aCGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCaa                                                                                                                 >  1:2916192/1‑86 (MQ=255)
                                                                                                        gCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGt                                                                                                   >  2:5380444/1‑97 (MQ=255)
                                                                                                        gCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGt                                                                                                   <  1:5380444/97‑1 (MQ=255)
                                                                                                            gctgTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCAc                                          <  2:4942008/150‑1 (MQ=255)
                                                                                                               gTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGa                                       <  2:79715/150‑1 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGTGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   >  2:2280039/1‑150 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   >  2:3626385/1‑150 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   >  1:134509/1‑150 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   >  1:2014385/1‑150 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   >  1:3430674/1‑150 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   <  2:3430674/150‑1 (MQ=255)
                                                                                                                   tACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTAc                                   >  2:2507104/1‑150 (MQ=255)
                                                                                                                     cccGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCctgctg                                               <  1:3091124/136‑1 (MQ=255)
                                                                                                                     cccGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGt                                 <  2:4195718/150‑1 (MQ=255)
                                                                                                                      ccGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTc                                <  2:1668118/150‑1 (MQ=255)
                                                                                                                      ccGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTc                                >  2:1293007/1‑150 (MQ=255)
                                                                                                                        ggTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTa                              >  2:3654870/1‑150 (MQ=255)
                                                                                                                         gTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTAc                             >  2:635637/1‑150 (MQ=255)
                                                                                                                         gTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTAc                             <  2:2930767/150‑1 (MQ=255)
                                                                                                                             aGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCAc                         >  1:2470191/1‑150 (MQ=255)
                                                                                                                              gTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACg                        >  1:1559010/1‑150 (MQ=255)
                                                                                                                               tATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGc                       <  1:3436776/150‑1 (MQ=255)
                                                                                                                                  gAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGcacc                    >  2:3282007/1‑150 (MQ=255)
                                                                                                                                  gAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGcacc                    >  2:3281939/1‑150 (MQ=255)
                                                                                                                                    gggCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGcaccac                  <  1:1234492/150‑1 (MQ=255)
                                                                                                                                        gAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCCTCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCCGCGCCCAGCGCTTCTTCCTGCTGGGCCCCGACTACGTCTACCCACGCCCCACcaac              >  2:3815136/1‑150 (MQ=255)
                                                                                                                                        gAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACcaac              >  1:2925331/1‑150 (MQ=255)
                                                                                                                                          agagATGTTGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACcaacaa            <  2:3119883/150‑1 (MQ=255)
                                                                                                                                           gagaTGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAg           <  1:3313176/150‑1 (MQ=255)
                                                                                                                                           gagaTGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCAAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAg           >  1:5130429/1‑150 (MQ=255)
                                                                                                                                                gTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCt      <  2:820805/150‑1 (MQ=255)
                                                                                                                                                gTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCt      <  2:3159329/150‑1 (MQ=255)
                                                                                                                                                gTCGCAGATCGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATctt      <  2:4316618/150‑3 (MQ=255)
                                                                                                                                                 tCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGcc                                                                 <  1:2671284/90‑1 (MQ=255)
                                                                                                                                                 tCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACTCACCACCAACAAGATCCTg     >  2:4314312/1‑150 (MQ=255)
                                                                                                                                                 tCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCAAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCTg     <  2:5130429/150‑1 (MQ=255)
                                                                                                                                                   gCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCTgcg   >  1:774195/1‑150 (MQ=255)
                                                                                                                                                   gCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCTgcg   >  1:5222874/1‑150 (MQ=255)
                                                                                                                                                    ccGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCTgcgc  >  2:765372/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAGAGCCGTCAACTGCTGACCCAGGACAAAGTGGCCGTGGTGTTCGGCTGCTGGACATCCGTGTCGCGCAAGTCTGTTCTGCCGGTGTTCGAAGAACTCAACGGCCTGCTGTTCTACCCGGTGCAGTATGAGGGCGAAGAGATGTCGCCGAACGTGTTCTACACAGGCGCCGCCCCGAACCAGCAAGCCATCCCGGCCGTGGAATACCTGATGAGCGAGGACGGCGGCAGCGCCAAGCGCTTCTTCCTGCTGGGCACCGACTACGTCTACCCACGCACCACCAACAAGATCCTGCGC  >  NC_002947/5510629‑5510926

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