Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,966,333 A→T 10.8% intergenic (‑402/+28) PP_5232 ← / ← amtB hypothetical protein/ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,966,3330AT10.8% 230.3 / 15.6 83intergenic (‑402/+28)PP_5232/amtBhypothetical protein/ammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (35/39);  new base T (4/5);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

CGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGA  >  NC_002947/5966191‑5966476
                                                                                                                                               |                                                                                                                                               
cgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCggg                                                                                                                                           <  5:28688/149‑1 (MQ=255)
cgcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCggg                                                                                                                                           <  3:158807/149‑1 (MQ=255)
 gcATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCggg                                                                                                                                           >  5:115906/1‑148 (MQ=255)
 gcATTCTAGCGCGCCCGTCGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAATACTGGCAAAAAAATGCCCGGCAAGCCggg                                                                                                                                           >  6:98212/1‑148 (MQ=255)
  cATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCa                                                                                                                                         <  5:14874/149‑1 (MQ=255)
  cATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCa                                                                                                                                         >  1:335688/1‑149 (MQ=255)
         gcgcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCTTGCCGGGCATtttttt                                                                                                                                  <  2:335688/149‑1 (MQ=255)
          cgcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTa                                                                                                                                 <  5:30211/149‑1 (MQ=255)
           gcgcCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGCCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTAc                                                                                                                                <  3:242034/149‑1 (MQ=255)
                      aaaaGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGt                                                                                                                     >  2:104908/1‑149 (MQ=255)
                      aaaaGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGt                                                                                                                     >  6:238231/1‑149 (MQ=255)
                       aaagcaaagcACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGt                                                                                                                     >  6:85905/1‑148 (MQ=255)
                         agcaaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTAc                                                                                                                  <  5:85905/149‑1 (MQ=255)
                          gcaaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAATACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGGCGCAGTTACa                                                                                                                 <  5:98212/149‑1 (MQ=255)
                           caaagcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAg                                                                                                                >  5:56849/1‑149 (MQ=255)
                              agcaCCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGtt                                                                                                             <  7:39236/149‑1 (MQ=255)
                                     tAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGCGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGcc                                                                                                       <  1:86030/149‑1 (MQ=255)
                                     tAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCcg                                                                                                      <  3:306410/149‑1 (MQ=255)
                                             cTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGtt                                                                                              <  8:42233/149‑1 (MQ=255)
                                             cTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGtt                                                                                              <  5:349310/149‑1 (MQ=255)
                                                cGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtg                                                                                             <  6:214179/147‑1 (MQ=255)
                                                 ggCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCg                                                                                                <  2:63129/143‑1 (MQ=255)
                                                 ggCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCg                                                                                                >  1:63129/1‑143 (MQ=255)
                                                  gCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtgtg                                                                                           >  2:216563/1‑147 (MQ=255)
                                                  gCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTtgtg                                                                                           <  1:216563/147‑1 (MQ=255)
                                                     tGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAg                                                                                      >  5:124953/1‑149 (MQ=255)
                                                     tGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAg                                                                                      >  1:109884/1‑149 (MQ=255)
                                                      gTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGc                                                                                                                        <  2:135708/114‑1 (MQ=255)
                                                      gTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGa                                                                                       <  5:121830/147‑1 (MQ=255)
                                                       tAAGTTTTTTACCTACGTCACTAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                                     <  7:57143/148‑1 (MQ=255)
                                                       tAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                                     >  7:267903/1‑148 (MQ=255)
                                                        aaGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGtt                                                                                              >  5:277101/1‑138 (MQ=255)
                                                        aaGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGtt                                                                                              <  6:277101/138‑1 (MQ=255)
                                                        aaGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGc                                                                                     >  8:57143/1‑147 (MQ=255)
                                                          gTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAg                                                                                      <  6:168257/144‑1 (MQ=255)
                                                          gTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAg                                                                                      >  5:168257/1‑144 (MQ=255)
                                                           ttttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGt                                                                                >  3:179866/1‑149 (MQ=255)
                                                           ttttttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAAATGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGt                                                                                <  6:298941/149‑1 (MQ=255)
                                                               ttACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGa                                                                             <  2:213295/148‑1 (MQ=255)
                                                                 acatccGTCACGAACGTGTCGAATAAACCGTTTGCCGGGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAgg                                                                           <  6:5481/143‑1 (MQ=255)
                                                                  ccTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGcc                                                                         >  8:30580/1‑149 (MQ=255)
                                                                    tACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccga                                                                       >  4:33447/1‑149 (MQ=255)
                                                                    tACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccga                                                                       >  5:72605/1‑149 (MQ=255)
                                                                        tCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGccgaccga                                                                   <  5:238231/149‑1 (MQ=255)
                                                                         cACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGAGCTCGTTGTGCTGAGCCAGGTCGAGGCcgaccgac                                                                  >  5:337795/1‑149 (MQ=255)
                                                                          aCGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGAct                                                                 >  7:334432/1‑149 (MQ=255)
                                                                                 tgtCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACtcttcttc                                                          >  3:40321/1‑149 (MQ=255)
                                                                                   tCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTcttcttctt                                                        <  7:30580/149‑1 (MQ=255)
                                                                                     gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAgg                                                                                 >  1:292567/1‑122 (MQ=255)
                                                                                     gAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAgg                                                                                 <  2:292567/122‑1 (MQ=255)
                                                                                          aaCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGAc                                                 >  6:94205/1‑149 (MQ=255)
                                                                                          aaCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGAc                                                 <  1:104908/149‑1 (MQ=255)
                                                                                          aaCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGAc                                                 >  7:232745/1‑149 (MQ=255)
                                                                                           aCCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGAcc                                                >  3:302212/1‑149 (MQ=255)
                                                                                               tttGACGTGCGGGAAATTCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTATTTACCGAGTCGCAGTTACAGGTTGTAGCCGCGCTCTTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACACg                                              >  6:287349/1‑147 (MQ=255)
                                                                                                ttGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAg                                           >  7:35474/1‑149 (MQ=255)
                                                                                                    cGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCTGCAGGcc                                        >  5:333250/1‑148 (MQ=255)
                                                                                                    cGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGccc                                       <  6:337795/149‑1 (MQ=255)
                                                                                                     gTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCa                                      >  4:161909/1‑149 (MQ=255)
                                                                                                      tGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGccc                                       >  8:354452/1‑147 (MQ=255)
                                                                                                      tGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGccc                                       <  7:354452/147‑1 (MQ=255)
                                                                                                        cGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATc                                    <  6:115906/148‑1 (MQ=255)
                                                                                                         gggAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcac                                  <  5:306514/149‑1 (MQ=255)
                                                                                                            aaaTGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATcaccac                               <  2:289888/149‑1 (MQ=255)
                                                                                                                 cTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATcc                           <  2:109884/148‑1 (MQ=255)
                                                                                                                  tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAg                         >  4:82305/1‑149 (MQ=255)
                                                                                                                  tGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAg                         >  1:17457/1‑149 (MQ=255)
                                                                                                                      aaaCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCAcc                     <  6:56849/149‑1 (MQ=255)
                                                                                                                        aCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCtt                   <  4:40321/149‑1 (MQ=255)
                                                                                                                        aCTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCtt                   <  4:179866/149‑1 (MQ=255)
                                                                                                                          tGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTga                 <  3:33447/149‑1 (MQ=255)
                                                                                                                          tGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGAACGACTCGTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCtt                   >  2:265375/1‑147 (MQ=255)
                                                                                                                            gCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTga                 <  8:334432/147‑1 (MQ=255)
                                                                                                                                aaaaaTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCa           >  4:20799/1‑149 (MQ=255)
                                                                                                                                   cttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGt        >  7:279008/3‑149 (MQ=255)
                                                                                                                                    ttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTgaga               >  7:274594/2‑141 (MQ=255)
                                                                                                                                    ttGCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTgaga               <  8:274594/140‑1 (MQ=255)
                                                                                                                                    aTGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACg         >  6:86369/1‑147 (MQ=255)
                                                                                                                                     tGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg      <  4:302212/149‑1 (MQ=255)
                                                                                                                                     tGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAg      <  3:82305/149‑1 (MQ=255)
                                                                                                                                      gCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTc                                                                                                 <  3:162629/57‑1 (MQ=255)
                                                                                                                                      gCCCGGCTTGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTc                                                                                                 >  4:162629/1‑57 (MQ=255)
                                                                                                                                         cGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGa  <  6:72605/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CGCATTCTAGCGCGCCCGTGGAAAAAGCAAAGCACCATAAACGCCCTGCGGCTTGTAAGTTTTTTACCTACGTCACGAACGTGTCGAATAAACCGTTTGCCGTGCGGGAAATGCTGAAAAACTGGCAAAAAAATGCCCGGCAAGCCGGGCAAGTTTTTACCGCGTCGCAGTTACAGGTTGTAGCCGCGCTCGTTGTGCTGAGCCAGGTCGAGGCCGACCGACTCTTCTTCTTCGTTGACCCGCAGGCCCATCACCACATCCAGCACCTTGAGAATCACGTAGGTGA  >  NC_002947/5966191‑5966476

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