Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,967,628 A→C 6.9% 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 6.9% 257.8 / 8.6 87L22R (CTG→CGG) amtBammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base C (3/3);  total (43/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 = 5.23e-01

ATCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGAT  >  NC_002947/5967484‑5967774
                                                                                                                                                |                                                                                                                                                  
atcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGc                                                                                                                                                >  6:219749/1‑149 (MQ=255)
atcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAggg                                                                                                                                                 >  1:289166/1‑148 (MQ=255)
 tcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGc                                                                                                                                                >  5:99246/1‑148 (MQ=255)
 tcaCCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGc                                                                                                                                                >  6:276399/1‑148 (MQ=255)
    ccGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAgg                                                                                                                                                  >  3:273032/1‑143 (MQ=255)
    ccGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAgg                                                                                                                                                  <  4:273032/143‑1 (MQ=255)
    ccGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGc                                                                                                                                            <  4:37639/149‑1 (MQ=255)
      gACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGccc                                                                                                                                          <  2:289166/149‑1 (MQ=255)
         aGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                       <  8:175789/149‑1 (MQ=255)
          gCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCggg                                                                                                                                       >  7:246926/1‑148 (MQ=255)
           cACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTTGCCCTGGCa                                                                                                                                     <  5:219749/149‑1 (MQ=255)
            aCGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCAt                                                                                                                                    >  4:217771/1‑149 (MQ=255)
             cGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTTGCCCTGGCaga                                                                                                                                   <  5:121130/149‑3 (MQ=255)
                tCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTAc                                                                                                                                 >  7:72686/1‑148 (MQ=255)
                  ttGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                              >  7:201545/1‑149 (MQ=255)
                  ttGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                              >  4:216837/1‑149 (MQ=255)
                   tGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAgg                                                                                                                             >  8:6769/1‑149 (MQ=255)
                    ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGTATAGTCATGAACAGCACCAGCGCCGTTGCTGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGa                                                                                                                               <  3:226726/146‑1 (MQ=255)
                    ggATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                              <  5:257042/147‑1 (MQ=255)
                         gCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAAt                                                                                                                       <  3:216837/149‑1 (MQ=255)
                           aCCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAg                                                                                                                     <  2:154333/149‑1 (MQ=255)
                           aCCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAg                                                                                                                     >  2:3607/1‑149 (MQ=255)
                             cATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAggg                                                                                                                   <  3:217771/149‑1 (MQ=255)
                              aTGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGc                                                                                                                  <  4:164246/149‑1 (MQ=255)
                                gccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCCGGGCTTGCCCTGGCaga                                                                                                                                   >  8:65623/1‑128 (MQ=255)
                                gccgccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTc                                                                                                                >  1:294145/1‑149 (MQ=255)
                                    ccATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAg                                                                                                            <  5:32021/149‑1 (MQ=255)
                                         gAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTg                                                                                                       <  7:6769/149‑1 (MQ=255)
                                          aaCAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGc                                                                                                      >  8:167915/1‑149 (MQ=255)
                                           aCAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGc                                                                                                      >  6:49233/1‑148 (MQ=255)
                                           aCAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGc                                                                                                      >  5:186785/1‑148 (MQ=255)
                                             aGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTAGGAGCCCtattgt                                                                                                   <  6:99246/149‑6 (MQ=255)
                                              ggtccAAGCCCGGAATAGTCATGAACTGCACGAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGGGCTGGGGCAGCTTCCTCTGCCAGGGCAAGGCCGGGCATTACGAGGTACATTAGGGCTCCTAGCCCTGCGATc                                                                                                  <  2:294145/146‑1 (MQ=255)
                                                 ccAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTa                                                                                               <  7:167915/149‑1 (MQ=255)
                                                 ccAAGCCCGGAATAGTCATGAACAGCACAAGCGCCGTTGCCGTCAGCATACAGGCGGTGTCGCCGGAATTAAGAGCTGGGTCAGCTTCCTATGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCATAGCCCTGCGATCTTa                                                                                               >  1:159172/1‑149 (MQ=255)
                                                     gCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCa                                                                                           <  5:276399/149‑1 (MQ=255)
                                                             tAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtg                                                                                   <  5:49233/149‑1 (MQ=255)
                                                                tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                                >  5:156017/1‑149 (MQ=255)
                                                                tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                                >  7:157778/1‑149 (MQ=255)
                                                                tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                                >  2:8236/1‑149 (MQ=255)
                                                                   tGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCGTAGCCCTGCGATCTTACGCAGAGTCATgttgtt                                                                              <  5:75321/148‑1 (MQ=255)
                                                                   tGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATgttgtt                                                                              >  6:75321/1‑148 (MQ=255)
                                                                   tGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGttt                                                                             >  3:300487/1‑149 (MQ=255)
                                                                      aCAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtct                                                                          >  8:32918/1‑149 (MQ=255)
                                                                      aCAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTtct                                                                          >  8:37731/1‑149 (MQ=255)
                                                                        aGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTcc                                                                        >  8:54485/1‑149 (MQ=255)
                                                                            ccAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTggg                                                                    <  8:246926/149‑1 (MQ=255)
                                                                              aGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGc                                                                  <  7:37731/149‑1 (MQ=255)
                                                                              aGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGc                                                                  <  7:32918/149‑1 (MQ=255)
                                                                                cgcCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGt                                                                <  1:16252/149‑1 (MQ=255)
                                                                                  ccGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTg                                                              <  8:72686/149‑1 (MQ=255)
                                                                                   cGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGTGACAATAGTGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTgg                                                             <  2:159172/149‑1 (MQ=255)
                                                                                   cGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTgg                                                             >  3:61808/1‑149 (MQ=255)
                                                                                   cGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTgg                                                             >  7:180204/1‑149 (MQ=255)
                                                                                        cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTg                                                        >  3:64793/1‑149 (MQ=255)
                                                                                        cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTg                                                        <  1:3607/149‑1 (MQ=255)
                                                                                         ggTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTgg                                                       <  3:143936/149‑1 (MQ=255)
                                                                                          gTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGt                                                      <  8:157778/149‑1 (MQ=255)
                                                                                             aGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTCTTCTCCTGGGGCGTTGGGTTTGGTGa                                                    <  1:8236/148‑1 (MQ=255)
                                                                                               cATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATgttgtt                                                                              <  6:213064/120‑1 (MQ=255)
                                                                                               cATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATgttgtt                                                                              >  5:213064/1‑120 (MQ=255)
                                                                                               cATACAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCACGACCGCGCATTACGAGGGACAATAGGGCTCCTAGCCATGCGATCTTACGCAGAGTCAtgttggt                                                                              >  5:213067/1‑118 (MQ=255)
                                                                                                      gCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGtttt                                                                            <  4:155095/115‑1 (MQ=255)
                                                                                                            tCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGa                                                    <  5:226093/133‑1 (MQ=255)
                                                                                                            tCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGa                                                    >  6:226093/1‑133 (MQ=255)
                                                                                                              gCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTg                                                              >  6:165202/1‑121 (MQ=255)
                                                                                                               ccGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTc                                                                                      >  1:244856/1‑96 (MQ=255)
                                                                                                               ccGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTc                                                                                      <  2:244856/96‑1 (MQ=255)
                                                                                                                 ggAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATc                               <  8:180204/149‑1 (MQ=255)
                                                                                                                        aGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCg                         <  6:156017/148‑1 (MQ=255)
                                                                                                                           gcTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGtt                                                               >  3:81251/1‑107 (MQ=255)
                                                                                                                           gcTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGtt                                                               <  4:81251/107‑1 (MQ=255)
                                                                                                                            cTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATc                    <  4:300487/149‑1 (MQ=255)
                                                                                                                            cTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATc                    <  8:201545/149‑1 (MQ=255)
                                                                                                                                ggCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGtt                <  4:61808/149‑1 (MQ=255)
                                                                                                                                 gCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCtttttt                                           >  7:158132/1‑121 (MQ=255)
                                                                                                                                 gCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCtttttt                                           <  8:158132/121‑1 (MQ=255)
                                                                                                                                 gCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGttt               >  8:282829/1‑149 (MQ=255)
                                                                                                                                  cAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTc              <  7:54485/149‑1 (MQ=255)
                                                                                                                                      tcgtaatGCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATc                               >  3:117150/7‑128 (MQ=255)
                                                                                                                                      tcgtaatGCCCGGGCTTGCCCTGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATc                               <  4:117150/122‑1 (MQ=255)
                                                                                                                                         ctctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTa       >  6:65143/1‑149 (MQ=255)
                                                                                                                                          tctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTAc      >  1:180797/1‑149 (MQ=255)
                                                                                                                                          tctGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTAGGTa                      >  8:257468/1‑133 (MQ=255)
                                                                                                                                            tGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTAc      <  4:64793/147‑1 (MQ=255)
                                                                                                                                              ccAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGAt  >  1:214528/1‑149 (MQ=255)
                                                                                                                                              ccAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCATGCGATCTTACGCAGAGTAATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGATTAGATCGCGTCGGTATCGGTTTCGCAGGTACGGAt  >  7:49195/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ATCACCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGAT  >  NC_002947/5967484‑5967774

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