Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,967,633 A→T 6.0% L20L (CTT→CTA) ‡ amtB ← ammonium transporter AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,967,6330AT 6.0% 247.9 / 3.2 82L20L (CTT→CTA) ‡amtBammonium transporter AmtB
Reads supporting (aligned to +/- strand):  ref base A (36/41);  new base T (2/3);  total (38/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 = 9.61e-01

CCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGCGG  >  NC_002947/5967488‑5967778
                                                                                                                                                 |                                                                                                                                                 
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                                                                                                                                  >  4:216837/1‑149 (MQ=255)
              ttGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAg                                                                                                                                  >  7:201545/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                                                                                    >  7:157778/1‑149 (MQ=255)
                                                            tCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCAtgttg                                                                                    >  5:156017/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                                                            <  1:3607/149‑1 (MQ=255)
                                                                                    cGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTg                                                            >  3:64793/1‑149 (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)
                                                                                                                                              ggCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGCgg  <  2:214528/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCGACAGCACGTTCTTGGATCGCACCATGCCGCCATAGAACAGGGCCAAGCCCGGAATAGTCATGAACAGCACCAGCGCCGTTGCGGTCAGCATCCAGGCGGTGTCGCCGGAATTAAGCGCTGGGGCAGCTTCCTCTGCCAGGGCAAGCCCGGGCATTACGAGGGACAATAGGGCTCCTAGCCCTGCGATCTTACGCAGAGTCATGTTGTTTTCTCCTGGGGCGTTGGGTTTGGTGAGGCTTTTTTGCTGCTTAGATCGCGTCGGTATCGGTTTCGCCGGTACGGATGCGG  >  NC_002947/5967488‑5967778

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: