Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 886,082 T→C 6.3% W665R (TGG→CGG) ‡ PP_0767 → LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947886,0820TC 6.3% 186.1 / 2.4 63W665R (TGG→CGG) ‡PP_0767LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

GGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGCA  >  NC_002947/885943‑886225
                                                                                                                                           |                                                                                                                                               
ggggCTGGGCGAGGCGCGAGCGTGGCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGc                                                                                                                                        <  2:104386/149‑1 (MQ=255)
   gCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCCCAGATGTGcatc                                                                                                                                     <  7:143340/149‑4 (MQ=255)
     tGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCCCAGATGTGcatcaa                                                                                                                                   <  3:49552/149‑6 (MQ=255)
        tagAGGCTCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGTAGTTCACTTTGCCGAGGGCTTTGCTGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTc                                                                                                                                <  4:58933/147‑1 (MQ=255)
           aGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTctac                                                                                                                             >  8:95154/1‑149 (MQ=255)
           aGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTctac                                                                                                                             <  5:111708/149‑1 (MQ=255)
             gcgcGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTctacta                                                                                                                           >  8:51952/1‑149 (MQ=255)
                cgAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCt                                                                                                                        >  6:110432/1‑149 (MQ=255)
                      tGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCAt                                                                                                                  >  4:216129/1‑149 (MQ=255)
                        ccGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATg                                                                                                                 <  8:25977/148‑1 (MQ=255)
                          gggACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGAt                                                                                                               >  5:84947/1‑148 (MQ=255)
                           ggACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTa                                                                                                             >  7:219886/1‑149 (MQ=255)
                               aTCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCt                                                                                                         >  7:128629/1‑149 (MQ=255)
                               aTCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCt                                                                                                         >  7:125630/1‑149 (MQ=255)
                               aTCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCt                                                                                                         >  3:81157/1‑149 (MQ=255)
                                 cAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTgg                                                                                                       >  3:55128/1‑149 (MQ=255)
                                  aGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGt                                                                                                      >  6:41817/1‑149 (MQ=255)
                                  aGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGt                                                                                                      <  4:108441/149‑1 (MQ=255)
                                    tgtgCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCAGGCCATGATTACCCTGGTc                                                                                                     >  5:97388/1‑148 (MQ=255)
                                      tgCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTc                                                                                                     <  3:215501/146‑1 (MQ=255)
                                      tgCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTc                                                                                                     >  4:215501/1‑146 (MQ=255)
                                          cATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGc                                                                                               >  4:174896/1‑148 (MQ=255)
                                           aTCGGTCATTGCGTGATTGCCACTCTCGATGGCCGGGAAGGTCACTTTGCCGAGCCCTTTGCGGAGCTGGCAGAGGCCGAGCGCTTGAAGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCAATCATAACCATGGTCAAATCCg                                                                                              >  1:288588/1‑148 (MQ=255)
                                           aTCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCg                                                                                              >  6:353117/1‑148 (MQ=255)
                                           aTCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCg                                                                                              >  1:66221/1‑148 (MQ=255)
                                           aTCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGa                                                                                             <  5:110432/149‑1 (MQ=255)
                                            tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCa                                                                                                                   >  6:251004/1‑126 (MQ=255)
                                            tCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCa                                                                                                                   <  5:251004/126‑1 (MQ=255)
                                                  aTTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTAGCCTGGTCAAATGCGAACTGtgg                                                                                      <  7:51952/149‑1 (MQ=255)
                                                   ttGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGtggc                                                                                     <  3:174896/149‑1 (MQ=255)
                                                              ccACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAggg                                                                           <  5:353117/148‑1 (MQ=255)
                                                               cACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGc                                                                         >  4:128159/1‑149 (MQ=255)
                                                               cACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTCGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTGCCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGc                                                                         >  1:63883/1‑149 (MQ=255)
                                                                aCGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTg                                                                                                        <  2:16360/117‑1 (MQ=255)
                                                                aCGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTg                                                                                                        >  1:16360/1‑117 (MQ=255)
                                                                    tCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACg                                                                    >  1:321458/1‑149 (MQ=255)
                                                                      gATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGCTGCACATCTTGGATGTGCCCCCGGTCTACTACCTGGCCATGATTATCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGa                                                                  <  5:244298/149‑1 (MQ=255)
                                                                             gggAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTcgcg                                                           <  6:84947/149‑1 (MQ=255)
                                                                               gcaggggaCTTTGGCGCGGGCTTTGGGGTGCTGGGCGATGCCGATCGCTTGTTGCACATCTGGGATGTGGCCCCGGTCTACTTCCTTGCCATGTTTACCCTGGTCAAATTCGAACTTTGGCTGGCCCAGGGGCGTACGGACCTCGCGGa                                                         <  2:66221/142‑1 (MQ=255)
                                                                                   gTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTc                                                      <  8:128629/148‑1 (MQ=255)
                                                                                      aCTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTgg                                                  >  5:47953/1‑149 (MQ=255)
                                                                                        tttGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCt                                                >  6:131121/1‑149 (MQ=255)
                                                                                          tGCCGAGGCGTTTGGGGAGCTGGGCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCGCCGGTCTACTACCTGTCCATGATTTCCCTGGTCAAATGCGAAGTGTGGCTGGCCCAGGGGCGTACGGACGTCGCGGAGTCCTGGCTg                                               <  2:288588/148‑1 (MQ=255)
                                                                                          tGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCt                                                <  8:125630/147‑1 (MQ=255)
                                                                                             cGAGGCCTTTGCGGAGCTGGCCGAGGCCGCGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGTGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGc                                                 >  4:39472/1‑143 (MQ=255)
                                                                                             cGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGTGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGc                                                 <  3:39472/143‑1 (MQ=255)
                                                                                                        cGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAg                                <  3:216129/149‑1 (MQ=255)
                                                                                                         ggAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTg                                               <  5:193618/133‑1 (MQ=255)
                                                                                                         ggAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTg                                               >  6:193618/1‑133 (MQ=255)
                                                                                                             cTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTa                           <  8:219886/149‑1 (MQ=255)
                                                                                                                gccgaggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATg                         <  3:128159/148‑1 (MQ=255)
                                                                                                                   gaggccgagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATggcgg                     <  4:81157/149‑1 (MQ=255)
                                                                                                                         gagCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGagca               <  2:330419/149‑1 (MQ=255)
                                                                                                                           gcgcTCGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTGCCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGagcagc             <  2:63883/149‑1 (MQ=255)
                                                                                                                            cgcTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGcc            <  8:140241/149‑1 (MQ=255)
                                                                                                                               ttGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggc         >  3:109945/1‑149 (MQ=255)
                                                                                                                                 ggggCACATCCCAGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcgg       >  8:81966/4‑149 (MQ=255)
                                                                                                                                  aTGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcggc      >  3:56447/1‑149 (MQ=255)
                                                                                                                                   ggCACATCCCAGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTACGGCGGCGAGCAGCcggcggcg     >  2:73392/2‑149 (MQ=255)
                                                                                                                                    gCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTgg                                                  <  3:183982/103‑1 (MQ=255)
                                                                                                                                    gCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTgg                                                  >  4:183982/1‑103 (MQ=255)
                                                                                                                                    gCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCcggcggcg     <  2:321458/148‑1 (MQ=255)
                                                                                                                                      acaTCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGCa  >  4:269173/1‑149 (MQ=255)
                                                                                                                                           tGGGATGTGCCACCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTTGCTGGCCCAGGGGAGTAAGGACCTCGCCGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATgc                        >  6:156592/1‑121 (MQ=255)
                                                                                                                                           tGGGAGGTGCCCCCGGTCTACTACCTGCCCATGATTACCCTGGTCAAATGCGAAATGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCAGCAGCGCCTGGGCCAGACCTATgg                        <  5:156592/122‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
GGGGCTGGGCGAGGCGCGAGCGTGCCGGGACATCAGTGTGCTCATCGGTCATTGCGTGATTGCCACGCTCGATGGCCGGGAAGGTCACTTTGCCGAGGCCTTTGCGGAGCTGGCCGAGGCCGAGCGCTTGATGCACATCTGGGATGTGCCCCCGGTCTACTACCTGGCCATGATTACCCTGGTCAAATGCGAACTGTGGCTGGCCCAGGGGCGTACGGACCTCGCGGAGTCCTGGCTGCAGCGCCTGGGCCAGACCTATGGCGGCGAGCAGCCGGCGGCGGCA  >  NC_002947/885943‑886225

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: ‑