Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,039,659 T→A 8.4% P172P (CCT→CCA) ‡ PP_3560 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,039,6590TA 8.4% 183.1 / 5.8 66P172P (CCT→CCA) ‡PP_3560LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base A (2/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGC  >  NC_002947/4039514‑4039803
                                                                                                                                                 |                                                                                                                                                
ccTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTTAGCTCGATGCAGCGATCACCACCGGACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAg                                                                                                                                               >  1:98877/1‑149 (MQ=255)
     ggCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCCGAGCTGGc                                                                                                                                          >  2:87900/1‑149 (MQ=255)
      gCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGcc                                                                                                                                         >  5:133811/1‑149 (MQ=255)
        aTGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCgg                                                                                                                                       <  2:6150/149‑1 (MQ=255)
          gTCGGCGCAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCACTTCGAGGACCTTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCggg                                                                                                                                      >  4:53172/1‑148 (MQ=255)
           tCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCTCAGGGCTGGcaatggc                                                                                                                                    <  8:20749/149‑7 (MQ=255)
           tCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGACGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCTCAGGGCTGGcaatggc                                                                                                                                    <  8:265161/149‑7 (MQ=255)
             ggcggAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGtata                                                                                                                                  <  7:127384/149‑1 (MQ=255)
                    cTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGcc                                                                                                                           <  7:260100/149‑1 (MQ=255)
                        cgcggcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAg                                                                                                                       >  4:213441/1‑149 (MQ=255)
                         gcggcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGt                                                                                                                      >  4:39662/1‑149 (MQ=255)
                            gcgcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCCGAGCTGGCCGGGTATAGCGTGCCGCAGTTGc                                                                                                                   <  1:87900/149‑1 (MQ=255)
                              gcgTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTg                                                                                                                 >  1:293847/1‑149 (MQ=255)
                                   gAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccag                                                                                                            >  6:317272/1‑149 (MQ=255)
                                    aaGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTGCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccagc                                                                                                           <  5:85397/149‑1 (MQ=255)
                                     aGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGgccagcc                                                                                                          <  6:134457/149‑1 (MQ=255)
                                             aGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTc                                                                                                  <  1:308623/149‑1 (MQ=255)
                                                  gATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTg                                                                                               <  2:293847/147‑1 (MQ=255)
                                                   aTGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTgcgc                                                                                            <  4:201322/149‑1 (MQ=255)
                                                   aTGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGCACCGTTTCTGGGCCATTGCCAGCCCCGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTgcgc                                                                                            <  1:82441/149‑1 (MQ=255)
                                                              accaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGc                                                                                  <  6:37068/148‑1 (MQ=255)
                                                               ccaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCt                                                                                                                                             >  3:252180/1‑88 (MQ=255)
                                                               ccaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCt                                                                                                                                             <  4:252180/88‑1 (MQ=255)
                                                               ccaccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGt                                                                                >  5:172108/1‑149 (MQ=255)
                                                                caccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTg                                                                               >  1:281522/1‑149 (MQ=255)
                                                                caccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCAAGCCATTCCTGCGCTTCGACAAGCGTg                                                                               >  7:118792/1‑149 (MQ=255)
                                                                  ccGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGcc                                                                             <  3:39662/149‑1 (MQ=255)
                                                                    gAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTg                                                                           <  6:172108/149‑1 (MQ=255)
                                                                     aaCCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTgg                                                                          >  5:267990/1‑149 (MQ=255)
                                                                        cGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggcc                                                                       >  6:191354/1‑149 (MQ=255)
                                                                        cGGGCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTGCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggcc                                                                       <  4:95060/149‑1 (MQ=255)
                                                                          gTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccgg                                                                     >  5:243291/1‑149 (MQ=255)
                                                                             aGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGggccggcc                                                                   <  8:41007/148‑1 (MQ=255)
                                                                                 cTTACCCGCAGAACCTGGACTTCACACAGCTGTTCGAGGACCGTTTATGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGAGCCTGCGccggccgg                                                                 >  5:305650/1‑146 (MQ=255)
                                                                                    aCCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGa                                                             <  6:133811/147‑1 (MQ=255)
                                                                                     cccGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATc                                                           >  8:146807/1‑148 (MQ=255)
                                                                                     cccGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATc                                                           <  2:281522/148‑1 (MQ=255)
                                                                                       cGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATc                                                           >  5:250584/1‑146 (MQ=255)
                                                                                       cGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAg                                                        <  7:146807/149‑1 (MQ=255)
                                                                                       cGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAg                                                        >  2:114665/1‑149 (MQ=255)
                                                                                        gCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAg                                                        >  2:170282/1‑148 (MQ=255)
                                                                                         cAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCa                                                      <  6:243291/149‑1 (MQ=255)
                                                                                         cAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCa                                                      <  5:191354/149‑1 (MQ=255)
                                                                                             aCCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCGATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGaa                                                  <  2:98877/149‑1 (MQ=255)
                                                                                              ccTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAAc                                                 <  3:213441/149‑1 (MQ=255)
                                                                                               cTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAAc                                                 >  6:151164/1‑148 (MQ=255)
                                                                                               cTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAAc                                                 >  8:31365/1‑148 (MQ=255)
                                                                                                              cTGTTCGAGGACCGTTTCTGGGCCTTTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCAATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcaggt                                 <  3:53172/149‑1 (MQ=255)
                                                                                                              cTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCt                                                                                                >  4:71750/1‑86 (MQ=255)
                                                                                                              cTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCt                                                                                                <  3:71750/86‑1 (MQ=255)
                                                                                                                       gACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcaggt                                 >  1:16718/1‑140 (MQ=255)
                                                                                                                       gACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAgcaggt                                 <  2:16718/140‑1 (MQ=255)
                                                                                                                           gTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGaa                    <  1:114665/149‑1 (MQ=255)
                                                                                                                            tttCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAAc                   >  5:266777/1‑149 (MQ=255)
                                                                                                                             ttCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGaa                    <  6:250584/147‑1 (MQ=255)
                                                                                                                                 gggCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATg              >  8:159297/1‑149 (MQ=255)
                                                                                                                                  ggCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATgg             <  7:31365/149‑1 (MQ=255)
                                                                                                                                    ccATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGaa           >  6:168174/1‑149 (MQ=255)
                                                                                                                                      aTTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACt         <  8:118792/149‑1 (MQ=255)
                                                                                                                                       cggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAAc          >  6:35522/3‑147 (MQ=255)
                                                                                                                                       cggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAAc          <  5:35522/145‑1 (MQ=255)
                                                                                                                                        ggCCAGCTCAGGGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGCCAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGGCGCCAGCAGGGGCAGGTGTGCGAACTGCTGGAACTgg       >  7:240972/2‑149 (MQ=255)
                                                                                                                                         gCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTgg       <  7:180355/148‑1 (MQ=255)
                                                                                                                                          ccAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGAc     <  7:159297/149‑1 (MQ=255)
                                                                                                                                           cAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGa      >  4:267068/1‑147 (MQ=255)
                                                                                                                                               ccTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGc  >  3:20601/1‑147 (MQ=255)
                                                                                                                                                 tGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGaa           <  4:2164/136‑1 (MQ=255)
                                                                                                                                                 tGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGaa           >  3:2164/1‑136 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGC  >  NC_002947/4039514‑4039803

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