Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,997 T→A 10.4% K246M (AAG→ATG) ‡ PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,9970TA10.4% 214.7 / 13.0 77K246M (AAG→ATG) ‡PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (4/4);  total (37/40)
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.63e-01

GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCG  >  NC_002947/5182851‑5183142
                                                                                                                                                  |                                                                                                                                                 
gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg                                                                                                                                                 >  1:298174/1‑149 (MQ=255)
  gcgcagcagcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACGGCTTGCt                                                                                                                                               <  3:43285/145‑1 (MQ=255)
   tcgcgcagcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTg                                                                                                                                              <  6:416803/145‑1 (MQ=255)
      agcagcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGAt                                                                                                                                            >  7:140167/1‑148 (MQ=255)
         agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                        <  6:281606/149‑1 (MQ=255)
         agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                        <  8:266350/149‑1 (MQ=255)
           cTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACg                                                                                                                                       >  7:161399/1‑148 (MQ=255)
              gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                    <  3:82308/148‑4 (MQ=255)
              gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATCACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCa                                                                                                                                     <  8:264172/147‑1 (MQ=255)
                cTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                                 >  3:276089/1‑149 (MQ=255)
                 ttCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTTAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                                 >  2:401749/1‑148 (MQ=255)
                 ttCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                                 >  7:273473/1‑148 (MQ=255)
                 ttCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                                 >  5:232586/1‑148 (MQ=255)
                    cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCAct                                                                                                                             <  6:819/149‑1 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCaa                                                                                                                                         <  7:222423/136‑1 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCaa                                                                                                                                         >  8:222423/1‑136 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                            <  5:161198/149‑1 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                            <  6:13123/149‑1 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                            <  6:19108/149‑1 (MQ=255)
                       aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                          <  4:39567/149‑1 (MQ=255)
                         cTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGt                                                                                                                        <  6:244245/149‑1 (MQ=255)
                          ttCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                       >  8:135163/1‑149 (MQ=255)
                           tCGTACAGGCGTAAGGCCTTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                       <  8:140167/148‑1 (MQ=255)
                                 aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcc                                                                                                                >  8:258272/1‑149 (MQ=255)
                                   gCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                       >  5:25984/1‑140 (MQ=255)
                                   gCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                       <  6:25984/140‑1 (MQ=255)
                                    cGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAggg                                                                                                                                  <  5:294529/128‑1 (MQ=255)
                                    cGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAggg                                                                                                                                  >  6:294529/1‑128 (MQ=255)
                                         ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAacac                                                                                                        >  7:271340/1‑149 (MQ=255)
                                         ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAacac                                                                                                        <  4:276089/149‑1 (MQ=255)
                                                aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGcc                                                                                                  <  7:135163/148‑1 (MQ=255)
                                                aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCt                                                                                                 >  6:181316/1‑149 (MQ=255)
                                                  gCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCCCACCCCCAACACCCCGCCTTc                                                                                               >  5:329327/1‑149 (MQ=255)
                                                       ctcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACACCGCCTTCGGCtt                                                                                          >  2:237877/1‑149 (MQ=255)
                                                        tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                         >  5:35673/1‑149 (MQ=255)
                                                          gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAcccc                                                                                                     >  5:210118/1‑135 (MQ=255)
                                                          gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAcccc                                                                                                     <  6:210118/135‑1 (MQ=255)
                                                            cctgATCGTCATCCATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAgg                                                                                     >  8:92478/4‑149 (MQ=255)
                                                              cGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgc                                                                                    <  3:120335/148‑1 (MQ=255)
                                                               gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgc                                                                                  >  4:340507/1‑149 (MQ=255)
                                                               gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgc                                                                                  <  3:340507/149‑1 (MQ=255)
                                                                 tCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgcgc                                                                                <  3:48896/149‑1 (MQ=255)
                                                                              cGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGctct                                                                   <  8:273473/149‑1 (MQ=255)
                                                                               gATCAGGACCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCtt                                                                  >  8:352026/1‑149 (MQ=255)
                                                                                   aGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAg                                                              <  6:232586/149‑1 (MQ=255)
                                                                                    ggCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                             <  4:25611/149‑1 (MQ=255)
                                                                                    ggCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                             <  6:385145/149‑1 (MQ=255)
                                                                                       ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                          >  6:14615/1‑149 (MQ=255)
                                                                                           cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                       <  8:161399/148‑1 (MQ=255)
                                                                                           cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                       <  5:14615/148‑1 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                                     <  5:181316/147‑1 (MQ=255)
                                                                                               aacaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCa                                                   <  2:298174/148‑1 (MQ=255)
                                                                                               aacaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCa                                                   <  7:258272/148‑1 (MQ=255)
                                                                                               aacaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAg                                                  >  6:88878/1‑149 (MQ=255)
                                                                                                    cccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGCTCACTCTTGTAGCACACCCCCAACCCCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCATGCAgg                                                             >  7:376280/1‑133 (MQ=255)
                                                                                                          aaTGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTa                                       <  6:35673/149‑1 (MQ=255)
                                                                                                                cttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                             <  3:103244/121‑1 (MQ=255)
                                                                                                                 ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTg                                >  8:232187/1‑149 (MQ=255)
                                                                                                                      gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGc                           >  1:188698/1‑149 (MQ=255)
                                                                                                                      gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGc                           >  5:387299/1‑149 (MQ=255)
                                                                                                                        aGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCt                          <  4:84452/148‑1 (MQ=255)
                                                                                                                               cTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaaa                  <  2:125721/149‑1 (MQ=255)
                                                                                                                                  gCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACACGCCACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCt               >  1:103877/1‑149 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCTGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg             >  1:229656/1‑149 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaa                   >  4:382052/1‑143 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaa                   <  3:382052/143‑1 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg             >  4:124767/1‑149 (MQ=255)
                                                                                                                                    ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACCCCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCCGGCTCAGTGAAATCTTg             >  6:81120/4‑149 (MQ=255)
                                                                                                                                    ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg             >  8:334817/4‑149 (MQ=255)
                                                                                                                                     tgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGc            >  7:133106/3‑149 (MQ=255)
                                                                                                                                      aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc           >  8:125674/1‑149 (MQ=255)
                                                                                                                                      aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc           >  4:321790/1‑149 (MQ=255)
                                                                                                                                         ttGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTc        <  8:271340/149‑1 (MQ=255)
                                                                                                                                         ttGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTc        <  6:329327/149‑1 (MQ=255)
                                                                                                                                           gAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCt      >  4:389811/1‑149 (MQ=255)
                                                                                                                                              cAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGgc   <  1:401749/149‑1 (MQ=255)
                                                                                                                                               aGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGgcg  <  2:188698/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCG  >  NC_002947/5182851‑5183142

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