Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,039,661 A→G 7.1% E173G (GAG→GGG)  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,6610AG 7.1% 201.8 / 4.4 70E173G (GAG→GGG) PP_3560LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base G (2/3);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAGGAA  >  NC_002947/4039514‑4039809
                                                                                                                                                   |                                                                                                                                                    
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                                                            <  5:191354/149‑1 (MQ=255)
                                                                                         cAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCa                                                            <  6:243291/149‑1 (MQ=255)
                                                                                             aCCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCGATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGaa                                                        <  2:98877/149‑1 (MQ=255)
                                                                                              ccTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAAc                                                       <  3:213441/149‑1 (MQ=255)
                                                                                               cTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAAc                                                       >  8:31365/1‑148 (MQ=255)
                                                                                               cTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAAc                                                       >  6:151164/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                 >  3:2164/1‑136 (MQ=255)
                                                                                                                                                 tGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGaa                 <  4:2164/136‑1 (MQ=255)
                                                                                                                                                  gAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGaa                                                        >  3:200181/1‑96 (MQ=255)
                                                                                                                                                  gAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGaa                                                        <  4:200181/96‑1 (MQ=255)
                                                                                                                                                   aGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGGCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAGGaa  <  7:217586/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTCGGGCATGTCGGCGGAGCTGGCGCGGCGCGTGGAAGACGGTGAGCTCGATGCAGCGATCACCACCGAACCGGTCAGCCCTTACCCGCAGAACCTGGACTTCACACCGCTGTTCGAGGACCGTTTCTGGGCCATTGCCAGCCCTGAGCTGGCCGGGTATAGCGTGCCGCAGTTGCTGGCCAGCCAGCCATTCCTGCGCTTCGACAAGCGTGCCTGGGCCGGCCGGCAGATCGAGCAGGAACTACGCCGCCAGCAGGTGCAGGTGAGCGAACTGATGGAACTGGACAGCCAGGAA  >  NC_002947/4039514‑4039809

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