Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,967,628 A→C 5.1% L22R (CTG→CGG)  amtB ← ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,967,6280AC 5.1% 312.3 / 6.6 98L22R (CTG→CGG) amtBammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (42/51);  new base C (3/2);  total (45/53)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

ATCATCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGC  >  NC_002947/5967481‑5967776
                                                                                                                                                   |                                                                                                                                                     
atcatcACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAg                                                                                                                                                      <  6:302109/149‑1 (MQ=255)
atcatcACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAg                                                                                                                                                      <  5:179081/149‑1 (MQ=255)
 tcatcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGACCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAgg                                                                                                                                                     >  1:326474/1‑149 (MQ=255)
    tcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGc                                                                                                                                                   >  2:65218/1‑148 (MQ=255)
     caCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGc                                                                                                                                                   >  2:156418/1‑147 (MQ=255)
       ccGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGc                                                                                                                                               >  7:166905/1‑149 (MQ=255)
           cAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCgg                                                                                                                                          <  7:13174/149‑1 (MQ=255)
            aGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                         >  8:211465/1‑149 (MQ=255)
            aGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                         <  1:394101/149‑1 (MQ=255)
            aGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                         <  4:426967/149‑1 (MQ=255)
             gCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTTGCCCTGGc                                                                                                                                        <  7:223180/149‑1 (MQ=255)
                  ttCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTAc                                                                                                                                   >  5:302717/1‑149 (MQ=255)
                       ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGa                                                                                                                                 <  8:13590/146‑1 (MQ=255)
                       ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                                <  7:216332/147‑1 (MQ=255)
                          tCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGAc                                                                                                                            >  1:168288/1‑148 (MQ=255)
                             cACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCAGCAGCGGCGTTGCGGTCAGCATGCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGAGAATa                                                                                                                        <  2:233897/149‑1 (MQ=255)
                             cACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACaa                                                                                                                          >  5:417805/1‑147 (MQ=255)
                             cACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACaa                                                                                                                          >  6:233032/1‑147 (MQ=255)
                             cACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATa                                                                                                                        <  2:168288/149‑1 (MQ=255)
                              aCCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAg                                                                                                                       >  4:234171/1‑149 (MQ=255)
                                    ccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTcc                                                                                                                 >  6:23275/1‑149 (MQ=255)
                                    ccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTc                                                                                                                  >  6:17130/1‑148 (MQ=255)
                                    ccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTc                                                                                                                  <  3:399691/148‑1 (MQ=255)
                                            gAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTg                                                                                                         <  1:424291/149‑1 (MQ=255)
                                             aaCAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGc                                                                                                        <  3:423667/149‑1 (MQ=255)
                                             aaCAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTAGGAGCCCtat                                                                                                        >  4:87400/1‑147 (MQ=255)
                                               cAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGa                                                                                                      <  8:96724/149‑1 (MQ=255)
                                               cAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGa                                                                                                      <  1:65218/149‑1 (MQ=255)
                                                aGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGAt                                                                                                     <  3:53335/149‑1 (MQ=255)
                                                 gggCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTAGGAGCCCtattgtc                                                                                                    <  3:87400/149‑7 (MQ=255)
                                                  ggCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCt                                                                                                   >  4:295462/1‑149 (MQ=255)
                                                  ggCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCt                                                                                                   <  5:23275/149‑1 (MQ=255)
                                                  ggCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCt                                                                                                   <  5:17130/149‑1 (MQ=255)
                                                    ccAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTa                                                                                                 >  8:338098/1‑149 (MQ=255)
                                                     cAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTAc                                                                                                <  3:234171/149‑1 (MQ=255)
                                                     cAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTAc                                                                                                <  1:99815/149‑1 (MQ=255)
                                                      aaGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACg                                                                                               <  7:211465/149‑1 (MQ=255)
                                                           cGGAATATTCATTTAAAGCACCACGGCCGTTGCGGTTATGATCCAGGCGGGGTCGCCGGAATTAATCGCTGGGGGATCTTCCTCTTGCATGGCAAGCCCGGGCATTAAGATGGACAATACGGCTCCTAGCCCTGCGATCTTACGCagag                                                                                          <  8:379327/149‑1 (MQ=255)
                                                             gAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTc                                                                                        <  3:295462/149‑1 (MQ=255)
                                                               aTAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAt                                                                                      >  7:403457/1‑149 (MQ=255)
                                                               aTAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAt                                                                                      <  4:19190/149‑1 (MQ=255)
                                                                tAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtg                                                                                     <  7:338098/149‑1 (MQ=255)
                                                                tAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtg                                                                                     >  4:309456/1‑149 (MQ=255)
                                                                 aGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgt                                                                                    >  6:182484/1‑149 (MQ=255)
                                                                   tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                                  >  3:378456/1‑149 (MQ=255)
                                                                   tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                                  >  1:177002/1‑149 (MQ=255)
                                                                     aTGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCagag                                                                                          >  4:124185/1‑139 (MQ=255)
                                                                     aTGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCagag                                                                                          <  3:124185/139‑1 (MQ=255)
                                                                       gAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGtttt                                                                              >  7:39671/1‑149 (MQ=255)
                                                                        aaCAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtc                                                                             >  5:352321/1‑149 (MQ=255)
                                                                          cAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtctc                                                                           <  2:352600/149‑1 (MQ=255)
                                                                            gCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCATGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCt                                                                         <  5:298759/149‑1 (MQ=255)
                                                                            gCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCt                                                                         >  4:104144/1‑149 (MQ=255)
                                                                             cACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTg                                                                        >  1:69784/1‑149 (MQ=255)
                                                                             cACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTg                                                                        >  7:124478/1‑149 (MQ=255)
                                                                               ccAGTGCCGTTGCGGTCAGCATCCAGGCGGTGTCGGCGGTATTAAGCGCTGGGGGAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                      <  8:371130/149‑1 (MQ=255)
                                                                                 aGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGc                                                                    <  3:309456/149‑1 (MQ=255)
                                                                                      cGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTgg                                                               <  8:3339/149‑1 (MQ=255)
                                                                                      cGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTgg                                                               <  8:166905/149‑1 (MQ=255)
                                                                                          gCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGtt                                                            <  6:302717/148‑1 (MQ=255)
                                                                                           cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGttt                                                           <  4:194339/148‑1 (MQ=255)
                                                                                           cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTg                                                          <  4:378456/149‑1 (MQ=255)
                                                                                           cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTg                                                          <  2:326474/149‑1 (MQ=255)
                                                                                           cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTg                                                          <  3:104144/149‑1 (MQ=255)
                                                                                             gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTttggtttggt                                                        >  2:169194/1‑149 (MQ=255)
                                                                                             gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                        >  6:266907/1‑149 (MQ=255)
                                                                                             gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                        <  8:403457/149‑1 (MQ=255)
                                                                                             gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                        <  6:352321/149‑1 (MQ=255)
                                                                                                aGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGa                                                      <  7:91491/148‑1 (MQ=255)
                                                                                                aGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGa                                                      <  4:277145/148‑1 (MQ=255)
                                                                                                aGCATCCAGGCGGTGTCACCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAg                                                     <  7:209732/149‑1 (MQ=255)
                                                                                                  cATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGc                                                   >  7:109496/1‑149 (MQ=255)
                                                                                                   aTCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtctc                                                                           <  4:365988/124‑1 (MQ=255)
                                                                                                   aTCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtctc                                                                           >  3:365988/1‑124 (MQ=255)
                                                                                                    tCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCtt                                                 >  1:285418/1‑149 (MQ=255)
                                                                                                        ggCGGGGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCtttttt                                             <  5:43015/149‑1 (MQ=255)
                                                                                                             tgtCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTtgctgc                                        >  7:288667/1‑149 (MQ=255)
                                                                                                                 gCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTCCGCAGAGTCATGTTGTTTTCTCGTGGGGCGTTGGGTTTGGTGAGGCTTTTTTgctgct                                       >  4:320964/1‑146 (MQ=255)
                                                                                                                 gCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                        >  2:308015/1‑129 (MQ=255)
                                                                                                                 gCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGAACTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                        <  1:308015/129‑1 (MQ=255)
                                                                                                                 cccccaaTTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGTCAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTa                                     <  1:156418/143‑1 (MQ=255)
                                                                                                                      aaTTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATcgc                               >  3:96381/1‑149 (MQ=255)
                                                                                                                       aTTAATCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGATGGAAAATAGGGATCCTAGCCCTGCGATCTTACGAAGAGTCAGGTTGGTTTCTCCTCGCGCGTTGGGATTGGCGAGGATTTTTTTCTGCTTAGATcgc                               >  2:263173/1‑148 (MQ=255)
                                                                                                                        ttAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGt                             <  6:43998/149‑1 (MQ=255)
                                                                                                                        ttAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGt                             <  8:288667/149‑1 (MQ=255)
                                                                                                                          aaGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGATGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCg                           >  2:377761/1‑149 (MQ=255)
                                                                                                                           aGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCg                           <  8:124478/148‑1 (MQ=255)
                                                                                                                                tGGGGCAGCTTCCTCTGCCAGGGCAAGCTCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATc                      >  1:248034/1‑149 (MQ=255)
                                                                                                                                      aGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCg               >  1:391386/1‑149 (MQ=255)
                                                                                                                                        cTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGcc             <  5:25008/149‑1 (MQ=255)
                                                                                                                                           ccTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCgg           >  6:390101/1‑148 (MQ=255)
                                                                                                                                             tctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGt                             >  7:335206/1‑128 (MQ=255)
                                                                                                                                             tctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGt                             <  8:335206/128‑1 (MQ=255)
                                                                                                                                              atGCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACg       >  4:309617/2‑149 (MQ=255)
                                                                                                                                               tGCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCttttt                                              <  2:339403/109‑1 (MQ=255)
                                                                                                                                               tGCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCttttt                                              >  1:339403/1‑109 (MQ=255)
                                                                                                                                                gCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGa     >  3:188330/1‑149 (MQ=255)
                                                                                                                                                  cAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATg   >  6:171666/1‑149 (MQ=255)
                                                                                                                                                   aGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGc  <  6:417805/149‑1 (MQ=255)
                                                                                                                                                   aGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGc  <  5:266907/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATCATCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGC  >  NC_002947/5967481‑5967776

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: