Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,149,783 G→C 5.8% V79L (GTG→CTG)  PP_3654 → leucine‑responsive regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,149,7830GC 5.8% 218.8 / 3.5 69V79L (GTG→CTG) PP_3654leucine‑responsive regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base C (2/2);  total (35/34)
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.14e-01

CCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTC  >  NC_002947/4149645‑4149930
                                                                                                                                          |                                                                                                                                                   
cccTGCCCTCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACCCACGGCTGGCACCGAGCGCCCTTTGGCTGAACTTTGCCGCGCTGGTATTCCTGACGCTACGCGAGGTAACCCGCCACCCGGTGGCGGACtt                                                                                                                                           <  2:318346/149‑1 (MQ=255)
cccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCTCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACtt                                                                                                                                           <  8:262591/149‑1 (MQ=255)
 ccTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACttt                                                                                                                                          >  4:241615/1‑149 (MQ=255)
  cTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTg                                                                                                                                         <  6:135876/149‑1 (MQ=255)
   tGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTg                                                                                                                                         >  3:115811/1‑148 (MQ=255)
   tGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTg                                                                                                                                         >  5:156952/1‑148 (MQ=255)
      cATCGGCGCTTGAAGGCACTGGAAGGGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAAc                                                                                                                                     >  1:314368/1‑149 (MQ=255)
        tCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACgg                                                                                                                                   <  5:81649/149‑1 (MQ=255)
        tCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACgg                                                                                                                                   >  6:53287/1‑149 (MQ=255)
          tggtcTTGAAGGCAATGGTAGTGTTCGGAGCCATTCTTGGGTATTACGGACGGCAGGCACCGACCGGCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGAAGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGc                                                                                                                                  <  4:175069/144‑1 (MQ=255)
               ttGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCa                                                                                                                                 >  7:247332/1‑144 (MQ=255)
               ttGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCa                                                                                                                                 <  8:247332/144‑1 (MQ=255)
                    ggCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAg                                                                                                                        >  4:175995/1‑148 (MQ=255)
                      cACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATc                                                                                                                     <  5:53287/149‑1 (MQ=255)
                       aCTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATcc                                                                                                                    >  1:15691/1‑149 (MQ=255)
                         tGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCg                                                                                                                  >  5:217400/1‑149 (MQ=255)
                          ggAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCg                                                                                                                  >  5:6300/1‑148 (MQ=255)
                           gAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCg                                                                                                                  >  2:114279/1‑147 (MQ=255)
                              gagTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACAGGAGCAGATCCCGCAGAt                                                                                                             >  2:220988/1‑149 (MQ=255)
                                 tCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGt                                                                                                          >  6:79887/1‑149 (MQ=255)
                                 tCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGt                                                                                                          <  3:175995/149‑1 (MQ=255)
                                   cGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTc                                                                                                         <  3:241615/148‑1 (MQ=255)
                                     gAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGAGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCg                                                                                                        >  7:103688/1‑147 (MQ=255)
                                              tGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTc                                                                                                         <  3:237769/137‑1 (MQ=255)
                                              tGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTc                                                                                                         >  4:237769/1‑137 (MQ=255)
                                               gggCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGc                                                                                            >  2:264064/1‑149 (MQ=255)
                                                ggCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCgg                                                                                             <  6:156952/147‑1 (MQ=255)
                                                       acgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTTTTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGc                                                                                    <  1:124280/149‑1 (MQ=255)
                                                        cgcacgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGc                                                                                    >  1:50445/1‑148 (MQ=255)
                                                           acgGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAcc                                                                                >  2:289613/1‑149 (MQ=255)
                                                              gCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCgg                                                                             <  2:314368/149‑1 (MQ=255)
                                                                   cACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTAc                                                                        <  6:6300/149‑1 (MQ=255)
                                                                      cGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACctg                                                                     >  6:172704/1‑149 (MQ=255)
                                                                       gAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACctgc                                                                    <  2:15691/149‑1 (MQ=255)
                                                                           gcCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAc                                                                                 <  3:98999/132‑1 (MQ=255)
                                                                           gcCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAc                                                                                 >  4:98999/1‑132 (MQ=255)
                                                                             ccTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGc                                                                                    <  6:141472/127‑1 (MQ=255)
                                                                             ccTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGc                                                                                    >  5:141472/1‑127 (MQ=255)
                                                                                   gCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTa                                                         <  5:79887/148‑1 (MQ=255)
                                                                                   gCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAg                                                        >  6:171075/1‑149 (MQ=255)
                                                                                     tGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAccc                                                                               <  6:177885/124‑1 (MQ=255)
                                                                                     tGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAccc                                                                               >  5:177885/1‑124 (MQ=255)
                                                                                       aaCTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAc                                                                                 <  1:281846/120‑1 (MQ=255)
                                                                                       aaCTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGAc                                                                                 >  2:281846/1‑120 (MQ=255)
                                                                                            tGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGa                                                 <  1:114279/147‑1 (MQ=255)
                                                                                             gCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCt                                              >  4:110046/1‑149 (MQ=255)
                                                                                                  gcTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCgg                                         <  1:127793/149‑1 (MQ=255)
                                                                                                      gTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGtt                                     >  4:136271/1‑149 (MQ=255)
                                                                                                       tATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGttt                                    <  2:50445/149‑1 (MQ=255)
                                                                                                         ttCGTGACGCTACGCGAGGTAACCCGCCAGCCGTTGGCGGACTTTGAACCGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGGAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTc                                   <  5:250322/148‑1 (MQ=255)
                                                                                                            gTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGa                                <  5:172704/148‑1 (MQ=255)
                                                                                                             tGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAg                              >  3:97261/1‑149 (MQ=255)
                                                                                                                     cgcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGttt                                    >  2:63634/1‑135 (MQ=255)
                                                                                                                     cgcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGttt                                    <  1:63634/135‑1 (MQ=255)
                                                                                                                     cgcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGt                          >  8:297867/1‑145 (MQ=255)
                                                                                                                     cgcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGt                          <  7:297867/145‑1 (MQ=255)
                                                                                                                     cgcgAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATg                       <  1:264064/148‑1 (MQ=255)
                                                                                                                           aaagtccGCCACCGGCTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAg                  <  2:37921/142‑1 (MQ=255)
                                                                                                                           aaagtccGCCACCGGCTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAg                  >  1:37921/6‑147 (MQ=255)
                                                                                                                            aagtccGCCACCGGCTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAg                                 >  3:104658/5‑131 (MQ=255)
                                                                                                                            aagtccGCCACCGGCTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAg                                 <  4:104658/127‑1 (MQ=255)
                                                                                                                               cccGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGt                          >  7:125692/1‑135 (MQ=255)
                                                                                                                               cccGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGt                          <  8:125692/135‑1 (MQ=255)
                                                                                                                                ccGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGa           >  7:144754/1‑149 (MQ=255)
                                                                                                                                 cGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGAc          >  3:192989/1‑149 (MQ=255)
                                                                                                                                   ccAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCa        <  4:192989/149‑1 (MQ=255)
                                                                                                                                     aGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGt      <  5:171075/149‑1 (MQ=255)
                                                                                                                                        cGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTAtt   <  3:110046/149‑1 (MQ=255)
                                                                                                                                        cGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTAtt   >  3:52255/1‑149 (MQ=255)
                                                                                                                                         ggTGGCGGACTTTGAATCGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCGGGATTAGCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTc  <  1:220988/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CCCTGCCATCGGCGCTTGAAGGCACTGGAAGAGTCCGGAGCCATTCTGGGCTATCACGCACGGCTGGCACCGAGCGCCCTTGGGCTGAACTTTGCCGCGCTGGTATTCGTGACGCTACGCGAGGTAACCCGCCAGCCGGTGGCGGACTTTGAAACGGCACTGGCGCAGATCCCGCAGATTGTCGAAGCGCAGCGGCTGTTCGGCGACCCGGATTACCTGCTGCATGTGGTAGCCAAGGACCTGCCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTC  >  NC_002947/4149645‑4149930

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: