Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,553,216 T→A 5.5% L277Q (CTG→CAG)  PP_3137 → group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,553,2160TA 5.5% 270.4 / 4.7 91L277Q (CTG→CAG) PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base A (3/2);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

AGCCACCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTACAAGTGA  >  NC_002947/3553069‑3553364
                                                                                                                                                   |                                                                                                                                                    
agccaCCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTgg                                                                                                                                                    <  1:141120/150‑1 (MQ=255)
 gccaCCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGgc                                                                                                                                                   <  2:1380332/150‑1 (MQ=255)
  ccaCCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGgcg                                                                                                                                                  >  2:240607/1‑150 (MQ=255)
  ccaCCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGACTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTAGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGgcg                                                                                                                                                  >  1:1001552/1‑150 (MQ=255)
   caCCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGgcgc                                                                                                                                                 <  2:82944/150‑1 (MQ=255)
      ccGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGcgcc                                                                                                                                              <  1:1759580/150‑1 (MQ=255)
           cGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCgcggc                                                                                                                                         <  2:169412/150‑1 (MQ=255)
             gCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCGGGGCATAGCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGGTCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGcggcgg                                                                                                                                       >  2:118298/1‑150 (MQ=255)
                 tGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAg                                                                                                                                   >  1:1414341/1‑150 (MQ=255)
                  ggCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGt                                                                                                                                  <  2:2341027/150‑1 (MQ=255)
                   gCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGgcgcgc                                                                                                                                               <  2:557803/136‑1 (MQ=255)
                   gCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGgcgcgc                                                                                                                                               >  1:557803/1‑136 (MQ=255)
                   gCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTg                                                                                                                                 >  2:1381225/1‑150 (MQ=255)
                      gCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGTCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGAc                                                                                                                              <  1:663326/150‑1 (MQ=255)
                      gCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGAc                                                                                                                              >  2:1002898/1‑150 (MQ=255)
                        gTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGAcgc                                                                                                                            >  1:291862/1‑150 (MQ=255)
                           aGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGAcgcgcg                                                                                                                         <  2:2193911/150‑1 (MQ=255)
                            gTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCga                                                                                                                        <  1:1000781/150‑1 (MQ=255)
                             tCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGAGTTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCgag                                                                                                                       >  1:632905/1‑150 (MQ=255)
                                  tCGGGCGTTGTTGTGGGGGCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTg                                                                                                                  <  2:516433/150‑1 (MQ=255)
                                   cGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTgg                                                                                                                 <  1:1203346/150‑1 (MQ=255)
                                    gggCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGc                                                                                                                <  1:424053/150‑1 (MQ=255)
                                    gggCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGc                                                                                                                <  1:2077354/150‑1 (MQ=255)
                                     ggCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCt                                                                                                               <  2:1847492/150‑1 (MQ=255)
                                      gCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTg                                                                                                              >  1:1991876/1‑150 (MQ=255)
                                        cTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCa                                                                                                            >  2:1066548/1‑150 (MQ=255)
                                            ttGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGcc                                                                                                        >  1:583329/1‑150 (MQ=255)
                                              gtgGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCgt                                                                                                      <  2:429433/150‑1 (MQ=255)
                                                     cTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGCGCGCCAGgcgcacc                                                                                                                                              >  1:312123/1‑100 (MQ=255)
                                                     cTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGCGCGCCAGGCGcacc                                                                                                                                              <  2:312123/103‑1 (MQ=255)
                                                        ggCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTg                                                                                                              >  1:893802/1‑132 (MQ=255)
                                                        ggCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTg                                                                                                              <  2:893802/132‑1 (MQ=255)
                                                           aTACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTg                                                                                         >  2:412989/1‑150 (MQ=255)
                                                           aTACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTAGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTg                                                                                         <  1:2299395/150‑1 (MQ=255)
                                                           aTACCTTTGGTGATCGCGCTGGCCTGCCTGCTGAAGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTg                                                                                         >  1:466261/1‑150 (MQ=255)
                                                             aCCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTggg                                                                                       <  2:583329/150‑1 (MQ=255)
                                                               cTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGcc                                                                                                        >  2:526407/1‑131 (MQ=255)
                                                               cTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGcc                                                                                                        <  1:526407/131‑1 (MQ=255)
                                                                tttGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCaa                                                                                    <  2:211223/150‑1 (MQ=255)
                                                                tttGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCaa                                                                                    >  1:1968623/1‑150 (MQ=255)
                                                                      gATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTgg                                                                                                                 >  1:1078287/1‑115 (MQ=255)
                                                                      gATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTgg                                                                                                                 <  2:1078287/115‑1 (MQ=255)
                                                                      gATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTccc                                                                              >  1:862510/1‑150 (MQ=255)
                                                                       aTCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCGGCTGTTGCTGTTGCGCTACGCACTGGAGGCGGTGCGCCTGGGGCGCCGCGGGGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCg                                                                             >  2:1128002/1‑150 (MQ=255)
                                                                           cgcTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGaaa                                                                         <  2:1968623/150‑1 (MQ=255)
                                                                            gcTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAAt                                                                        >  1:1177611/1‑150 (MQ=255)
                                                                                 cctgcctgCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCg                                                                   <  2:1414341/150‑1 (MQ=255)
                                                                                  ctgcctgcTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCgg                                                                  >  2:1710457/1‑150 (MQ=255)
                                                                                       tgctgctgGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGgt                                                                                             <  2:2372747/118‑1 (MQ=255)
                                                                                       tgctgctgGGCGGCTGGGGCCTGCTGTTGCTGTCGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGgt                                                                                             >  1:2372747/1‑118 (MQ=255)
                                                                                        gctgctgGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACt                                                            <  1:1002898/150‑1 (MQ=255)
                                                                                             tgGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGt                                                       <  2:1177611/150‑1 (MQ=255)
                                                                                                  ggcggttgCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTg                                                  >  2:1801171/8‑150 (MQ=255)
                                                                                                       gggCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCg                                                                                                       <  1:2224409/92‑1 (MQ=255)
                                                                                                       gggCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCg                                                                                                       >  2:2224409/1‑92 (MQ=255)
                                                                                                         gCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGCCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCAcc                                           >  2:1735971/1‑150 (MQ=255)
                                                                                                              ctgttgctgttgATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCtt                                      <  1:1710457/150‑1 (MQ=255)
                                                                                                                gttgctgttgATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGc                                    >  2:1435961/1‑150 (MQ=255)
                                                                                                                gttgctgttgATCTACCCACTGCAGGCGGGGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTACAGGCCGTGTTCCTGGTGTTGGGCACGTTCCCCGAAATGCTCGGCGGACTGAATTGCCTGCGGCACCGCCGTGc                                    <  2:1001552/150‑1 (MQ=255)
                                                                                                                 ttgctgttgATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCg                                   >  2:2238177/1‑150 (MQ=255)
                                                                                                                  ggctgttgATCTACCCACGGCAGGCGGTGCGGCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCATCTGAAGTTCCTGCGCCACCGCCTTGCgg                                  <  2:721386/149‑1 (MQ=255)
                                                                                                                    ctgttgATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCt                                >  1:354139/1‑150 (MQ=255)
                                                                                                                     tgttgATCTACCCACTGCGGGCGGTGCGCCGGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGTTg                               <  2:501935/150‑1 (MQ=255)
                                                                                                                     tgttgATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTg                               <  1:412989/150‑1 (MQ=255)
                                                                                                                     tgttgATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTg                               <  1:1735971/150‑1 (MQ=255)
                                                                                                                          aTCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCaaa                          >  2:740456/1‑150 (MQ=255)
                                                                                                                          aTCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCaaa                          >  1:1909601/1‑150 (MQ=255)
                                                                                                                           tCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAAt                         >  1:1020339/1‑150 (MQ=255)
                                                                                                                            tctcccACTGCAGGCGGTGCGCCTGGCGCGCCTCGGTGGCAAGTGGTCGCGCGAGAACTGGCTGCTGGCCGTGTTCCTGGTGTTGGGCAAGTCCCCGCACATGCTCGGCCAAAAGAAGTTGCTGCCCCACCGCCTTGCGGCGGCaaaatc                       <  2:87776/147‑4 (MQ=255)
                                                                                                                              aCCCACTGCAGGCGGTGAGCCTGGCGCGACGCGGCTGCAAGTGGACGCGCTAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCgg                      <  2:632905/150‑1 (MQ=255)
                                                                                                                               cccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTccc                                                                              <  2:2337023/93‑1 (MQ=255)
                                                                                                                               cccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTccc                                                                              >  1:2337023/1‑93 (MQ=255)
                                                                                                                               cccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGc                     >  1:1199924/1‑150 (MQ=255)
                                                                                                                               cccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGc                     <  1:240607/150‑1 (MQ=255)
                                                                                                                               cccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGc                     >  2:108425/1‑150 (MQ=255)
                                                                                                                                ccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCgaga                    >  1:87776/1‑147 (MQ=255)
                                                                                                                                ccACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCt                    >  1:1860037/1‑150 (MQ=255)
                                                                                                                                  aCTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCttt                  >  1:503424/1‑150 (MQ=255)
                                                                                                                                  aCTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCttt                  >  2:66727/1‑150 (MQ=255)
                                                                                                                                      gccgcggcgCGCCAGGCGCACCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCGACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATc              >  1:2108300/4‑150 (MQ=255)
                                                                                                                                      cAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATc              <  2:354139/150‑1 (MQ=255)
                                                                                                                                        ggCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGa            <  2:503424/150‑1 (MQ=255)
                                                                                                                                        ggCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGa            >  1:2361246/1‑150 (MQ=255)
                                                                                                                                        ggCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGa            <  1:1381225/150‑1 (MQ=255)
                                                                                                                                        ggCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGa            <  2:1538487/150‑1 (MQ=255)
                                                                                                                                         gcggcgCGCCAGGCGCACCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTg                               <  1:1994556/130‑1 (MQ=255)
                                                                                                                                         gcggcgCGCCAGGCGCACCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTg                               >  2:1994556/1‑130 (MQ=255)
                                                                                                                                                gcCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTACAAGt    >  1:1385165/1‑150 (MQ=255)
                                                                                                                                                gcCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTACAAGt    <  1:2238177/150‑1 (MQ=255)
                                                                                                                                                 ccTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTACAAGTg   <  2:291862/150‑1 (MQ=255)
                                                                                                                                                  cTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGCGTACAAGTga  <  2:1199924/150‑1 (MQ=255)
                                                                                                                                                   tGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTa         >  1:195407/1‑142 (MQ=255)
                                                                                                                                                   tGGCGCGCCGCGGCGGCAAGGGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTa         <  2:195407/142‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCACCCGAGCGGCATTGGCTGCGTGAGTCACGTCGGGCCTGGTTGTGGGGCCTGGGCATACCTTTGGTGATCGCGCTGGCCTGCCTGCTGCTGGGCGGCTGGGGCCTGCTGTTGCTGTTGATCTACCCACTGCAGGCGGTGCGCCTGGCGCGCCGCGGCGGCAAGTGGACGCGCGAGAACTGGCTGCAGGCCGTGTTCCTGGTGTTGGGCAAGTTCCCGGAAATGCTCGGCCAACTGAAGTTCCTGCGCCACCGCCTTGCGGCTGGCAAATCGGCTTTGATCGAGTACAAGTGA  >  NC_002947/3553069‑3553364

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: