Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,141,427 T→A 8.9% L374M (TTG→ATG) ‡ PP_1897 → DNA internalization‑related competence protein ComEC/Rec2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,141,4270TA 8.9% 239.7 / 12.7 81L374M (TTG→ATG) ‡PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base A (4/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

GCCTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTG  >  NC_002947/2141279‑2141576
                                                                                                                                                     |                                                                                                                                                     
gcCTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCtt                                                                                                                                                      >  2:1437632/1‑150 (MQ=255)
gcCTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCtt                                                                                                                                                      >  1:817933/1‑150 (MQ=255)
    gggATCTGGCGGTCGTTCTCCGCCGTGGCCACGCTCATCTATTGCTTGTGGGGGCTGTTGGGGGGTTGGCGGCCCTGGCAGGCCTGGACAAGAGCGCAATGGTGGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCTGCTTGCCg                                                                                                                                                  >  2:342682/1‑150 (MQ=255)
       ttCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCAAGCCGGTg                                                                                                                                               <  2:1275400/150‑1 (MQ=255)
          tGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACATGAGCGCAATGGATGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGc                                                                                                                                            <  2:1733724/150‑1 (MQ=255)
          tGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGc                                                                                                                                            <  1:1286868/150‑1 (MQ=255)
          tGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTGGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCTCACCGGCAAGCCGGTggc                                                                                                                                             <  1:310602/150‑3 (MQ=255)
                  gTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTgc                                                                                                                                    >  2:378608/1‑150 (MQ=255)
                   ttCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTgcg                                                                                                                                   >  1:1606473/1‑150 (MQ=255)
                       ccgccgTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGcgccgc                                                                                                                               >  1:927809/1‑150 (MQ=255)
                         gccgTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTg                                                                                                                             <  2:110637/150‑1 (MQ=255)
                          ccgTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTgg                                                                                                                            <  1:133465/150‑1 (MQ=255)
                                 cACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCt                                                                                                                     >  1:211294/1‑150 (MQ=255)
                                 cACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCAGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCt                                                                                                                     >  1:1180215/1‑150 (MQ=255)
                                   cGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCtt                                                                                                                                                      <  2:614807/115‑1 (MQ=255)
                                   cGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCtt                                                                                                                                                      >  1:614807/1‑115 (MQ=255)
                                   cGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTtg                                                                                                                   <  2:66162/150‑1 (MQ=255)
                                    gCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTtgt                                                                                                                  <  2:817933/150‑1 (MQ=255)
                                      tCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTgg                                                                                                                <  2:211294/150‑1 (MQ=255)
                                         tCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCgg                                                                                                             >  2:1225734/1‑150 (MQ=255)
                                            attGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGACAGGCATGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCTCCGCTGGCCAACCTTGTGGCGGTGc                                                                                                          <  1:57150/148‑1 (MQ=255)
                                            aCTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGc                                                                                                          >  1:123767/1‑150 (MQ=255)
                                              tGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCtt                                                                                                                    >  2:2032224/1‑138 (MQ=255)
                                              tGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCtt                                                                                                                    <  1:2032224/138‑1 (MQ=255)
                                              tGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACATTGTGACCGTGCCg                                                                                                        >  2:1377827/1‑150 (MQ=255)
                                                cTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTg                                                                                                      >  1:362989/1‑150 (MQ=255)
                                                cTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTg                                                                                                      >  2:255362/1‑150 (MQ=255)
                                                 ttCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTgg                                                                                                     <  1:255362/150‑1 (MQ=255)
                                                 ttCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTgg                                                                                                     <  2:2223015/150‑1 (MQ=255)
                                                  tCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTggg                                                                                                    <  1:1248116/150‑1 (MQ=255)
                                                  tCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTggg                                                                                                    <  2:1811963/150‑1 (MQ=255)
                                                   cAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGt                                                                                                   >  1:336184/1‑150 (MQ=255)
                                                      tGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCAAAGGCCTCCCCTGGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAg                                                                                                <  2:593133/150‑1 (MQ=255)
                                                       gCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGc                                                                                               <  1:1437632/150‑1 (MQ=255)
                                                       gCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCCCCTCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGc                                                                                               <  1:1998780/150‑1 (MQ=255)
                                                           ggTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTTGCGGTTCCGTGGGTCAGCCTgg                                                                                           <  2:1835970/150‑1 (MQ=255)
                                                                ggtgcatGGCGCCCCTGGAAGGCCTGGACACGAGCGCAAGGGGTGGTCGCCGTCGGCTTGCTAGTGGTGTTGCTGGCTACAAGGTTGCCGGTGAGCCTCAGTGCGGCGCTGGCGAACCTTGTGGCGGTGCGGTGGGGCAGCCTGGCGATc                                                                                      >  2:500963/7‑150 (MQ=255)
                                                                    gtTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGc                                                                                  <  2:1606473/150‑1 (MQ=255)
                                                                     ttGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGt                                                                                                            >  1:2259696/1‑123 (MQ=255)
                                                                     ttGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGt                                                                                                            <  2:2259696/123‑1 (MQ=255)
                                                                        gcggcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGtt                                                                              <  1:1178122/150‑1 (MQ=255)
                                                                          ggcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTgg                                                                            <  2:659288/150‑1 (MQ=255)
                                                                          ggcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCAACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTgg                                                                            >  2:2396879/1‑150 (MQ=255)
                                                                          ggcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTAGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTgg                                                                            >  1:2048358/1‑150 (MQ=255)
                                                                           gcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGc                                                                           >  1:933969/1‑150 (MQ=255)
                                                                           gcCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGc                                                                           >  2:2044276/1‑150 (MQ=255)
                                                                               tGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGcc                                                                                                                                                   <  2:653480/74‑1 (MQ=255)
                                                                               tGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGcc                                                                                                                                                   >  1:653480/1‑74 (MQ=255)
                                                                                  cAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTg                                                                    <  1:2169115/150‑1 (MQ=255)
                                                                                   aGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGTCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTgg                                                                   <  1:562233/150‑1 (MQ=255)
                                                                                   aGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGGCTGGCAGGCATGCCGTTGGCGTTGCttg                                                                   >  2:1157595/1‑148 (MQ=255)
                                                                                   aGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTgg                                                                   <  2:480691/150‑1 (MQ=255)
                                                                                      ccTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCa                                                                >  1:1165296/1‑150 (MQ=255)
                                                                                         ggACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGACGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGc                                                             <  2:2047801/150‑1 (MQ=255)
                                                                                          gACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCt                                                            >  2:2134531/1‑150 (MQ=255)
                                                                                                 gcgcAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGc                                                                                          <  1:1582660/113‑1 (MQ=255)
                                                                                                 gcgcAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGc                                                                                          >  2:1582660/1‑113 (MQ=255)
                                                                                                 gcgcAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTg                                                     <  2:336184/150‑1 (MQ=255)
                                                                                                     aaTGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGc                                                 >  1:747248/1‑150 (MQ=255)
                                                                                                     aaTGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGc                                                 <  1:342682/150‑1 (MQ=255)
                                                                                                      aTGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGTCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCt                                                <  2:55043/150‑1 (MQ=255)
                                                                                                          gTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGt                                                                                                                                      <  1:2341277/60‑1 (MQ=255)
                                                                                                          gTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGt                                                                                                                                      >  2:2341277/1‑60 (MQ=255)
                                                                                                           tGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCgg                                           <  1:378608/150‑1 (MQ=255)
                                                                                                                gCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTGGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGtt                                      >  2:2333749/1‑150 (MQ=255)
                                                                                                                      ggCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAg                                <  1:2134531/150‑1 (MQ=255)
                                                                                                                       gCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAgg                               >  2:211281/1‑150 (MQ=255)
                                                                                                                        cTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGc                              <  2:123767/150‑1 (MQ=255)
                                                                                                                         ttGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGAGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGAGAGGCg                             >  2:916688/1‑150 (MQ=255)
                                                                                                                          tGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGt                            >  1:2108252/1‑150 (MQ=255)
                                                                                                                              aCCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGct                        >  2:1689585/1‑150 (MQ=255)
                                                                                                                               ccGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGctc                       <  2:747248/150‑1 (MQ=255)
                                                                                                                                       gCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCtggctggc               >  1:365646/1‑150 (MQ=255)
                                                                                                                                        cTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCa              <  1:1229366/150‑1 (MQ=255)
                                                                                                                                           gctcACCGGCAAGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAgg            >  2:1804524/4‑150 (MQ=255)
                                                                                                                                           gCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGt           >  2:1377249/1‑150 (MQ=255)
                                                                                                                                             tcACCGGCAAGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTg                                                                             >  2:687575/2‑83 (MQ=255)
                                                                                                                                             tcACCGGCAAGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTg                                                                             <  1:687575/82‑1 (MQ=255)
                                                                                                                                             tcACCGGCAAGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTg          >  2:1143618/2‑150 (MQ=255)
                                                                                                                                             tcACCGGCAAGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTg          >  1:1614873/2‑150 (MQ=255)
                                                                                                                                              cACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGACGTTGGCGGTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTgg         >  1:2332872/1‑150 (MQ=255)
                                                                                                                                                 cGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCct      >  2:1050696/1‑150 (MQ=255)
                                                                                                                                                   gCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGt           >  2:560002/1‑143 (MQ=255)
                                                                                                                                                   gCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGt           <  1:560002/143‑1 (MQ=255)
                                                                                                                                                   gCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTAAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCctgc    <  1:500963/150‑1 (MQ=255)
                                                                                                                                                     ttGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCctgctg  >  1:1917306/1‑150 (MQ=255)
                                                                                                                                                     ttGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCctgctg  <  2:927809/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTGGGTTCTGGCTGTCGTTCGCCGCCGTGGCCACGCTCATCTACTGCTTCAGTGCCCGGTTGGGTGGTTGGCGGCCCTGGCAGGCCTGGACACGAGCGCAATGGGTGATCGCCGTCGGCTTGCTACCGGTGTTGCTGGCCACCGGCTTGCCGGTGAGCCTCAGTGCGCCGCTGGCCAACCTTGTGGCGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTG  >  NC_002947/2141279‑2141576

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: