Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,029 G→T 6.8% Y87* (TAC→TAA) ‡ PP_3717 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,0290GT 6.8% 213.0 / 3.3 68Y87* (TAC→TAA) ‡PP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (28/35);  new base T (2/3);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTTGTCGGGGGCGATTTCCTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACT  >  NC_002947/4239882‑4240177
                                                                                                                                                   |                                                                                                                                                    
gaagTTGTCGGGGGCGATTTCCTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTgg                                                                                                                                                      >  4:63107/1‑148 (MQ=255)
gaagTTGTCGGGGGCGATTTCCGGCAGGTCGTAGAAGCCTGCGGCCAGGCCGCTTTGCACGGCAAGGCAGAACGGGGCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGt                                                                                                                                                     <  6:73294/149‑1 (MQ=255)
  agTTGTCGGGGGCGATTTCCTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTAc                                                                                                                                                   <  7:263980/149‑1 (MQ=255)
           gggCGATTTCCTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTgc                                                                                                                                          >  5:246527/1‑149 (MQ=255)
                    ccTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCAGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCAGGTACAGCAACTGCGGcgcccgc                                                                                                                                 >  2:105670/1‑149 (MQ=255)
                     cTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGc                                                                                                                                                <  8:265477/133‑1 (MQ=255)
                     cTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGc                                                                                                                                                >  7:265477/1‑133 (MQ=255)
                     cTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGcgcccgc                                                                                                                                 <  3:385/148‑1 (MQ=255)
                       gCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGTACCAGCaaggcat                                                                                                                                         >  4:48290/1‑133 (MQ=255)
                       gCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGTACCAGCaaggcat                                                                                                                                         <  3:48290/138‑6 (MQ=255)
                          ggtggtAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTgg                                                                                                                           <  4:109483/149‑1 (MQ=255)
                           gtggtAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTgg                                                                                                                           <  8:165450/148‑1 (MQ=255)
                           gtggtAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTgg                                                                                                                           >  7:165450/1‑148 (MQ=255)
                           gtgctAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGc                                                                                                                          <  4:108504/149‑1 (MQ=255)
                             ggtAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCAc                                                                                                                        >  2:352173/1‑149 (MQ=255)
                                aGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTc                                                                                                                      >  7:141121/1‑148 (MQ=255)
                                 gAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGTAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCaa                                                                                                                    <  5:341015/149‑1 (MQ=255)
                                  aaGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAg                                                                                                                   <  2:283872/149‑1 (MQ=255)
                                          ggCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGatc                                                                                                            >  6:99394/1‑148 (MQ=255)
                                           gCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGTACCAGCaag                                                                                                                                             <  7:200390/114‑2 (MQ=255)
                                           gCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGTACCAGCaag                                                                                                                                             >  8:200390/1‑113 (MQ=255)
                                             cAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTg                                                                                                               >  6:201194/1‑142 (MQ=255)
                                             cAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTg                                                                                                               <  5:201194/142‑1 (MQ=255)
                                              aGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGatcatcat                                                                                                       >  3:221830/1‑149 (MQ=255)
                                               ggCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGatcatcatc                                                                                                      <  6:51318/149‑1 (MQ=255)
                                                  cGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGAc                                                                                                   >  5:158402/1‑149 (MQ=255)
                                                   gCTTTGGACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGatca                                                                                                           <  7:251295/140‑1 (MQ=255)
                                                   gCTTTGGACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCCCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGatca                                                                                                           >  8:251295/1‑140 (MQ=255)
                                                    cTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTc                                                                                                 <  3:63107/149‑1 (MQ=255)
                                                       tGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTcc                                                                                                <  8:141121/147‑1 (MQ=255)
                                                            ggCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCa                                                                                         >  5:256404/1‑149 (MQ=255)
                                                             gCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCAc                                                                                        <  6:256404/149‑1 (MQ=255)
                                                              cAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCAc                                                                                        >  7:265560/1‑148 (MQ=255)
                                                                  gCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACt                                                                                   >  5:171891/1‑149 (MQ=255)
                                                                    aGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCa                                                                                                                         >  1:36499/1‑109 (MQ=255)
                                                                    aGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCa                                                                                                                         <  2:36499/109‑1 (MQ=255)
                                                                     gAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGc                                                                                <  6:158402/149‑1 (MQ=255)
                                                                      aaaGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGTCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGcc                                                                               >  5:2507/1‑149 (MQ=255)
                                                                      aaaGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGcc                                                                               <  5:201133/149‑1 (MQ=255)
                                                                      aaaGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGcc                                                                               >  6:201133/1‑149 (MQ=255)
                                                                       aaGGGTCCAGCAGGTAGCCGGCGCGGAACTAGCGCTCCCGCACCGCGGCCCGATGCAGTTCGGGGATGCCTTGCTGTTACAGCAACTGCGGCGCCCGCCCCTGGCACTCACGGCTGAGCATCATCGACTCCACCTGCACCAACTGGCCg                                                                              <  8:99890/149‑1 (MQ=255)
                                                                          ggTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATg                                                                           <  1:352173/149‑1 (MQ=255)
                                                                           ggccAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATGATCATCGACTCCACCTGCACCAACTGGCCGATgg                                                                          <  6:85120/147‑1 (MQ=255)
                                                                                  agGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGc                                                                   <  8:280841/149‑1 (MQ=255)
                                                                                  agGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCGGAGCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGc                                                                   <  4:297770/149‑1 (MQ=255)
                                                                                       gCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACAt                                                              <  4:221830/149‑1 (MQ=255)
                                                                                          ggCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACt                                                           <  2:183386/149‑1 (MQ=255)
                                                                                             gcTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAg                                                                                                                   >  7:162656/1‑90 (MQ=255)
                                                                                             gcTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAg                                                                                                                   <  8:162656/90‑1 (MQ=255)
                                                                                             gcTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTc                                                          >  1:330052/1‑147 (MQ=255)
                                                                                                 aaaTAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAgcg                                                     <  6:171891/148‑1 (MQ=255)
                                                                                                 aaaTAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAgc                                                      <  6:246527/147‑1 (MQ=255)
                                                                                                  aaaaGCGCTCCAGCACCGCGGCCCGAAGCCGTTAGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCAAATACTCCAgcgc                                                    >  8:298192/4‑148 (MQ=255)
                                                                                                       cgcTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCggg                                              <  8:265560/149‑1 (MQ=255)
                                                                                                         cTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCCACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCggggg                                            >  7:233373/1‑149 (MQ=255)
                                                                                                         cTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATgg                                                                          <  4:296527/119‑1 (MQ=255)
                                                                                                         cTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATgg                                                                          >  3:296527/1‑119 (MQ=255)
                                                                                                             aGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATgg                                                                          <  1:231781/115‑1 (MQ=255)
                                                                                                             aGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATgg                                                                          >  2:231781/1‑115 (MQ=255)
                                                                                                             aGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTc                                        <  5:99394/149‑1 (MQ=255)
                                                                                                                        ccGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCa                             <  8:233373/149‑1 (MQ=255)
                                                                                                                                    ggggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGCTGCCATCGGCcagttgg                                                            <  3:349024/106‑6 (MQ=255)
                                                                                                                                    ggggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGCTGCCATCGGCcagttgg                                                            >  4:349024/1‑101 (MQ=255)
                                                                                                                                    ggggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGc                  >  3:179651/1‑148 (MQ=255)
                                                                                                                                       gATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACg              >  1:66172/1‑149 (MQ=255)
                                                                                                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGcc            <  7:89917/146‑1 (MQ=255)
                                                                                                                                            cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGcc            >  8:89917/1‑146 (MQ=255)
                                                                                                                                             ttGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCa                                    <  3:341585/121‑1 (MQ=255)
                                                                                                                                             ttGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCa                                    >  4:341585/1‑121 (MQ=255)
                                                                                                                                                   gTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGCCTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGGTTCCCCATGGGCCAGGGCGCGCTGCAACGCCCCGTGCCAct  >  2:315266/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAAGTTGTCGGGGGCGATTTCCTGCAGGTGGTAGAAGCCTTCGGCCAGGCCGCTTTGCACGGCAAGGCAGAAAGGGTCCAGCAGGTAGCCGGCGCTGAAATAGCGCTCCAGCACCGCGGCCCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACT  >  NC_002947/4239882‑4240177

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: