Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,729,020 G→A 6.2% A250T (GCC→ACC)  PP_2390 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,729,0200GA 6.2% 193.5 / 2.0 64A250T (GCC→ACC) PP_2390LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (2/2);  total (32/32)
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.75e-01

GGCGTTGGCTGGCGCGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAA  >  NC_002947/2728877‑2729168
                                                                                                                                                 |                                                                                                                                                    
ggcgTTGGCTGGCGCGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGgccgcc                                                                                                                                                 >  1:70244/1‑149 (MQ=255)
       gctggcGCGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGcc                                                                                                                                          >  7:41593/1‑149 (MQ=255)
           gcgcGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGt                                                                                                                                      >  7:10896/1‑149 (MQ=255)
                    gATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCttt                                                                                                                             <  8:44301/149‑1 (MQ=255)
                           tGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCAt                                                                                                                      >  6:125386/1‑149 (MQ=255)
                             gACATCAGCAGCGACCAGGTGCAGGTGTTCGATACCCTCGACCCTACCCCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCAt                                                                                                                      <  1:283740/149‑1 (MQ=255)
                                aTCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGt                                                                                                                 >  6:181599/1‑149 (MQ=255)
                                 tCAGCATCGACCAGGGGCAGGTGTTCGATAGCCTCGAAC‑‑ATGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTg                                                                                                                <  2:245639/149‑1 (MQ=255)
                                 tCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTg                                                                                                                <  6:124767/149‑1 (MQ=255)
                                 tCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTg                                                                                                                <  4:241292/149‑1 (MQ=255)
                                 tCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTg                                                                                                                >  1:199188/1‑149 (MQ=255)
                                 tCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTg                                                                                                                <  2:250205/149‑1 (MQ=255)
                                     cATCGACCAGGGGCAGGTGTTCGATACCCTCGCCC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCa                                                                                                            >  7:178826/1‑149 (MQ=255)
                                       tCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCAcc                                                                                                          >  2:141771/1‑149 (MQ=255)
                                         gACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCGAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCgg                                                                                                        <  5:10751/149‑1 (MQ=255)
                                          aCCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGc                                                                                                       <  1:141771/149‑1 (MQ=255)
                                             aGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGAc                                                                                                    <  6:43480/149‑1 (MQ=255)
                                             aGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGAc                                                                                                    <  1:300637/149‑1 (MQ=255)
                                                   aGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTtac                                                                                              <  5:181599/149‑1 (MQ=255)
                                                       gtTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGGTACTACa                                                                                          >  8:287740/1‑149 (MQ=255)
                                                        ttCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACa                                                                                          <  2:199188/148‑1 (MQ=255)
                                                        ttCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACAt                                                                                         <  8:41593/149‑1 (MQ=255)
                                                          cGATACCCTCGCCC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCAACCTGGCCCCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGa                                                                                       >  8:72886/1‑149 (MQ=255)
                                                               cccTAGACC‑‑AGTGCATATCGGCCGCACAGCAAGGCCTGGGCATTTCAGTTGTAGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTgg                                                                                  <  5:182308/149‑1 (MQ=255)
                                                                 cTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCt                                                                                >  6:87340/1‑149 (MQ=255)
                                                                  tCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGACCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCt                                                                                >  5:229705/1‑148 (MQ=255)
                                                                      cc‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAgg                                                                           >  3:37465/1‑149 (MQ=255)
                                                                           gggCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGc                                                                          >  3:203827/1‑147 (MQ=255)
                                                                            ggCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGc                                                                                                              >  1:166159/1‑110 (MQ=255)
                                                                            ggCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGc                                                                                                              <  2:166159/110‑1 (MQ=255)
                                                                              cATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGcc                                                                     <  6:340239/149‑1 (MQ=255)
                                                                               atatCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGcc                                                                     >  4:94179/1‑148 (MQ=255)
                                                                                       cgccCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTgg                                                                                                               <  1:95578/98‑1 (MQ=255)
                                                                                       cgccCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTgg                                                                                                               >  2:95578/1‑98 (MQ=255)
                                                                                       cgccCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAgcgc                                                            <  7:5547/149‑1 (MQ=255)
                                                                                        gccCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGcgcg                                                           <  5:125386/149‑1 (MQ=255)
                                                                                         cccAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGc                                                                                 >  3:294122/1‑126 (MQ=255)
                                                                                         cccAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGc                                                                                 <  4:294122/126‑1 (MQ=255)
                                                                                           cAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCa                                                        <  1:257262/149‑1 (MQ=255)
                                                                                            aGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCgg                                                          >  8:220172/1‑146 (MQ=255)
                                                                                                aGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATg                                                    >  3:123489/1‑148 (MQ=255)
                                                                                                    cTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGcac                                                 >  3:159806/1‑147 (MQ=255)
                                                                                                    cTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACaa                                               <  4:37465/149‑1 (MQ=255)
                                                                                                       ggCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCt                                            >  5:175894/1‑149 (MQ=255)
                                                                                                           tttCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGt                                        <  6:175894/149‑1 (MQ=255)
                                                                                                            ttCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGt                                        >  1:322784/1‑148 (MQ=255)
                                                                                                                 tGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAg                                                                       <  6:284475/112‑1 (MQ=255)
                                                                                                                 tGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTc                                  <  8:10896/149‑1 (MQ=255)
                                                                                                                   gTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTcc                                 <  2:70244/148‑1 (MQ=255)
                                                                                                                   gTCGATCTGGTACTGGCCAGCACCGACCTGACCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTcc                                 <  6:229705/148‑1 (MQ=255)
                                                                                                                        tCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCg                           >  8:47838/1‑149 (MQ=255)
                                                                                                                            gTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCa                       <  7:47838/149‑1 (MQ=255)
                                                                                                                                  gCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAcc                 >  3:174144/1‑149 (MQ=255)
                                                                                                                                  gCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGCCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAcc                 >  4:63975/1‑149 (MQ=255)
                                                                                                                                   ccAGCACCGACCTGACCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCa                                                                        <  2:24640/93‑1 (MQ=255)
                                                                                                                                   ccAGCACCGACCTGACCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCa                                                                        >  1:24640/1‑93 (MQ=255)
                                                                                                                                    cAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGcaca                                                >  8:8747/1‑116 (MQ=255)
                                                                                                                                    cAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAccgc               >  7:156851/1‑149 (MQ=255)
                                                                                                                                    cAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAccgc               >  4:188808/1‑149 (MQ=255)
                                                                                                                                    cAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAccgc               >  8:32147/1‑149 (MQ=255)
                                                                                                                                       cACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACcgccgc            <  3:188808/149‑1 (MQ=255)
                                                                                                                                       cACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACcgccgc            >  2:173708/1‑149 (MQ=255)
                                                                                                                                            aCCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGcc       <  4:159806/149‑1 (MQ=255)
                                                                                                                                               tGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTc                                  <  4:193150/119‑1 (MQ=255)
                                                                                                                                               tGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTc                                  >  3:193150/1‑119 (MQ=255)
                                                                                                                                                ggCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTAc    >  8:66359/1‑148 (MQ=255)
                                                                                                                                                 gccgccGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACaa  <  7:287740/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCGTTGGCTGGCGCGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACC‑‑AGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAA  >  NC_002947/2728877‑2729168

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