Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,080,200 G→A 6.3% H193Y (CAC→TAC)  amaC ← D‑lysine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,080,2000GA 6.3% 186.2 / 2.1 63H193Y (CAC→TAC) amaCD‑lysine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

ACCAGGCACGGCACGCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAGGTGTGCACCT  >  NC_002947/4080055‑4080346
                                                                                                                                                 |                                                                                                                                                  
aCCAGGCACGGCACGCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGtgctg                                                                                                                                                 >  2:321835/1‑149 (MQ=255)
  cAGGCACGGCACGCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGAGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGtgccg                                                                                                                                                 >  3:321448/1‑145 (MQ=255)
    ggcacggcacGCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGtgctgctg                                                                                                                                              >  3:38891/1‑148 (MQ=255)
     gcacggcacgCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGtgctgctg                                                                                                                                              >  2:4514/1‑147 (MQ=255)
         ggcactCCGGCACGGGCAATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGTAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATc                                                                                                                                        <  3:68629/149‑1 (MQ=255)
             cgCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCa                                                                                                                                    <  4:38891/149‑1 (MQ=255)
                cGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCg                                                                                                                                 <  1:321835/149‑1 (MQ=255)
                cGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAAAGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCg                                                                                                                                 >  6:161683/1‑149 (MQ=255)
                                gcgGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGgcagc                                                                                                                 <  7:11615/149‑1 (MQ=255)
                                gcgGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGgcag                                                                                                                  >  2:176560/1‑148 (MQ=255)
                                  ggATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCAc                                                                                                               >  3:136839/1‑149 (MQ=255)
                                             gCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGcagcag                                                                                                    >  4:290433/1‑149 (MQ=255)
                                                  ttCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTACTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTgca                                                                                                        <  7:131817/140‑1 (MQ=255)
                                                  ttCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTACTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTgca                                                                                                        >  8:131817/1‑140 (MQ=255)
                                                     caccaGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGgca                                                                                                                   >  2:19870/1‑126 (MQ=255)
                                                     caccaGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGgca                                                                                                                   <  1:19870/126‑1 (MQ=255)
                                                       ccaGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTACTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTg                                                                                          >  2:183643/1‑149 (MQ=255)
                                                         aGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGtt                                                                                        <  1:176560/149‑1 (MQ=255)
                                                           gCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTgg                                                                                      >  3:169301/1‑149 (MQ=255)
                                                            ccTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTACTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGc                                                                                     <  1:183643/149‑1 (MQ=255)
                                                             cTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGcc                                                                                    >  2:121089/1‑149 (MQ=255)
                                                             cTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGcc                                                                                    >  3:346941/1‑149 (MQ=255)
                                                              ttCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCg                                                                                   >  1:78876/1‑149 (MQ=255)
                                                                   cGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGtgctgc                                                                                                                                                <  7:12389/83‑1 (MQ=255)
                                                                   cGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGtgctgc                                                                                                                                                >  8:12389/1‑83 (MQ=255)
                                                                    gAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGc                                                                             >  1:160201/1‑149 (MQ=255)
                                                                                         gAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGCCCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGc                                                        >  7:195615/1‑149 (MQ=255)
                                                                                         gAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGc                                                        >  8:44724/1‑149 (MQ=255)
                                                                                         gAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGATGCCACTGCTGCTGATGCAGATCGGCACCGGTAGGGTTCTGGCAGCACGGGTGAAGCAGTACCACGCTGTTGGCCGGCAGGGTATGCAAGGTGGCCAGCATGc                                                        >  6:124252/1‑149 (MQ=255)
                                                                                           ggAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGTCAGGCTCTGCAAGGTGGCCAGCATGCCg                                                      <  3:290433/149‑1 (MQ=255)
                                                                                             aaCGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTc                                                    <  1:5903/149‑1 (MQ=255)
                                                                                             aaCGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTc                                                    <  7:44724/149‑1 (MQ=255)
                                                                                                gggATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGaa                                                 <  4:346941/149‑1 (MQ=255)
                                                                                                  gATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGt                                               >  8:55727/1‑149 (MQ=255)
                                                                                                          ggCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTccacgcca                                       >  2:200054/1‑149 (MQ=255)
                                                                                                          ggCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTccacgcca                                       >  4:227452/1‑149 (MQ=255)
                                                                                                             gTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACgg                                    <  4:268073/149‑1 (MQ=255)
                                                                                                             gTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACgg                                    >  3:268073/1‑149 (MQ=255)
                                                                                                                ccTTGACCTCTTGGACCACTTCTTGCCACTGGTGGTGCTGCAGATCGGCTCCGGTAGTGTTGTGGCAGCACGTGTGGAGGAGTTCCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCGAGCATGCCGTCGAAGTCCACGCCACGGGt                                  <  5:124252/148‑1 (MQ=255)
                                                                                                                 cTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTgg                                >  3:259494/1‑149 (MQ=255)
                                                                                                                     aCCACTTCGCCCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTg                            <  8:195615/149‑1 (MQ=255)
                                                                                                                     aCCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGGAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTTGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCt                             <  5:161683/148‑1 (MQ=255)
                                                                                                                      ccACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTgg                           <  1:4514/149‑1 (MQ=255)
                                                                                                                       cACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAg                                                                                                                             >  8:208010/1‑50 (MQ=255)
                                                                                                                       cACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAg                                                                                                                             <  7:208010/50‑1 (MQ=255)
                                                                                                                       cACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGTAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGt                          >  8:236583/1‑149 (MQ=255)
                                                                                                                       cACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGt                          >  6:111646/1‑149 (MQ=255)
                                                                                                                       cACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACCcgtggcgtgga                          <  3:227452/149‑2 (MQ=255)
                                                                                                                       cACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAACTCCACGCCACCcgtggcgtgga                          <  7:55727/149‑2 (MQ=255)
                                                                                                                          ttCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTg                                                                  >  2:171442/1‑106 (MQ=255)
                                                                                                                          ttCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTg                                                                  <  1:171442/106‑1 (MQ=255)
                                                                                                                            cGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGTAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCaaa                      <  7:236583/148‑1 (MQ=255)
                                                                                                                             gACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCaaa                      <  4:136839/147‑1 (MQ=255)
                                                                                                                             gACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAg                     <  2:345514/148‑1 (MQ=255)
                                                                                                                             gACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAg                     >  5:223153/1‑148 (MQ=255)
                                                                                                                             gACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAg                     <  6:18620/148‑1 (MQ=255)
                                                                                                                              aCCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTa                   >  5:315125/1‑149 (MQ=255)
                                                                                                                               ccACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTGCACGCCACGGGTGGCCTGGTCAAAGTAc                  <  1:200054/149‑1 (MQ=255)
                                                                                                                                  cTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGa               <  4:259494/149‑1 (MQ=255)
                                                                                                                                      ttgCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAg            <  2:160201/148‑1 (MQ=255)
                                                                                                                                      ttgCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAg            <  6:315125/148‑1 (MQ=255)
                                                                                                                                        gCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAGgt          <  1:121089/148‑1 (MQ=255)
                                                                                                                                         ccACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAGgtgt        >  5:356016/1‑149 (MQ=255)
                                                                                                                                         ccACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAGgtg         <  4:169301/148‑1 (MQ=255)
                                                                                                                                          cACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGTAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAg                                                >  8:171408/1‑108 (MQ=255)
                                                                                                                                          cACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAg                                                <  7:171408/108‑1 (MQ=255)
                                                                                                                                           aCTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCaaa                      >  8:250826/1‑133 (MQ=255)
                                                                                                                                           aCTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCaaa                      <  7:250826/133‑1 (MQ=255)
                                                                                                                                               gTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAGGTGTGCACCt  >  8:227402/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACCAGGCACGGCACGCCGGCACGGGCCATTTCGCGGATGGCGTAGGCATCTTCCACCAGGCCTTCGGCGAAGCCTTGATAAGCGATGTCGAGGAACGGGATCAGCTGGCGTGCCTTGACCACTTCGACCACTTGTTGCCACTGGTGCTGCTGCAGATCGGCACCGGTAGGGTTGTGGCAGCACGGGTGCAGCAGTACCACGCTGTTGGCCGGCAGGCTCTGCAAGGTGGCCAGCATGCCGTCGAAGTCCACGCCACGGGTGGCCTGGTCAAAGTACGGATAGGTGTGCACCT  >  NC_002947/4080055‑4080346

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: