Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,729,030 G→A 6.8% R253H (CGC→CAC)  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,0300GA 6.8% 174.4 / 2.0 59R253H (CGC→CAC) PP_2390LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base A (2/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

CGCGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAAGGACATC  >  NC_002947/2728889‑2729175
                                                                                                                                             |                                                                                                                                                 
cgcgAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCAGGCGCCTGGTc                                                                                                                                            <  6:58898/149‑1 (MQ=255)
 gcgAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTgg                                                                                                                                              >  4:317923/1‑146 (MQ=255)
 gcgAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTgg                                                                                                                                              <  3:317923/146‑1 (MQ=255)
             gcTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAg                                                                                                                                >  5:181656/1‑148 (MQ=255)
                     tCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTg                                                                                                                       >  5:81965/1‑149 (MQ=255)
                        gCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCaa                                                                                                                                 >  4:74594/1‑136 (MQ=255)
                        gCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCaa                                                                                                                                 >  4:187732/1‑136 (MQ=255)
                        gCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCaa                                                                                                                                 <  3:74594/136‑1 (MQ=255)
                        gCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCaa                                                                                                                                 <  3:187732/136‑1 (MQ=255)
                          aTCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCAc                                                                                                                  >  7:163325/1‑149 (MQ=255)
                              aCCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGc                                                                                                              >  6:281956/1‑149 (MQ=255)
                                cAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGa                                                                                                            >  7:77890/1‑149 (MQ=255)
                                     gCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTTGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGtt                                                                                                       <  7:259856/149‑1 (MQ=255)
                                       aGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGACCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTtac                                                                                                     >  6:21749/1‑149 (MQ=255)
                                          tgtTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTtactac                                                                                                  <  6:81965/149‑1 (MQ=255)
                                               gATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGAc                                                                                             <  7:269320/149‑1 (MQ=255)
                                                aTACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGAcc                                                                                            >  7:282483/1‑149 (MQ=255)
                                                        gACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTc                                                                                      <  8:163325/147‑1 (MQ=255)
                                                            aGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGc                                                                                 <  7:194303/148‑1 (MQ=255)
                                                            aGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGc                                                                                 >  8:194303/1‑148 (MQ=255)
                                                               gCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGc                                                                             >  6:321730/1‑149 (MQ=255)
                                                                 atatCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGcc                                                                            >  5:296310/1‑148 (MQ=255)
                                                                    tCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCaa                                                                         >  2:83097/1‑148 (MQ=255)
                                                                      ggCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTAGATGACCTGGCTCAAGGCCAGCACaaaa                                                                       <  8:178017/148‑1 (MQ=255)
                                                                      ggCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCaaaa                                                                       <  8:182818/148‑1 (MQ=255)
                                                                      ggCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCaa                                                                         >  7:178017/1‑146 (MQ=255)
                                                                      ggCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACAGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCaa                                                                         >  7:182818/1‑146 (MQ=255)
                                                                         cgccCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAgcgc                                                                   >  4:273311/1‑149 (MQ=255)
                                                                          gccCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCTGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGcgcg                                                                  >  2:20295/1‑149 (MQ=255)
                                                                            ccAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTtac                                                                                                     >  3:270038/1‑112 (MQ=255)
                                                                            ccAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTtac                                                                                                     <  4:270038/112‑1 (MQ=255)
                                                                             cAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCGAGCCCAAAACCGCgg                                                                 >  4:135528/1‑147 (MQ=255)
                                                                              aGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCgg                                                                 >  8:202663/1‑146 (MQ=255)
                                                                              aGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAg                                                              <  4:315853/149‑1 (MQ=255)
                                                                                 aaGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGACCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATg                                                           <  5:21749/149‑1 (MQ=255)
                                                                                      cTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACaa                                                      <  8:77890/149‑1 (MQ=255)
                                                                                      cTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACaa                                                      <  6:181656/149‑1 (MQ=255)
                                                                                        gggCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGc                                                    >  1:231235/1‑149 (MQ=255)
                                                                                          gCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTg                                                  >  2:28958/1‑149 (MQ=255)
                                                                                                 aGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACt                                           <  6:296310/149‑1 (MQ=255)
                                                                                                 aGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACt                                           >  1:223769/1‑149 (MQ=255)
                                                                                                   tGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCTGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTc                                         <  1:20295/149‑1 (MQ=255)
                                                                                                   tGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATTACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACaa                                                      >  2:204978/1‑136 (MQ=255)
                                                                                                   tGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATTACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACaa                                                      <  1:204978/136‑1 (MQ=255)
                                                                                                      tCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTg                                      >  3:217852/1‑149 (MQ=255)
                                                                                                              gTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCg                                  <  3:274035/145‑1 (MQ=255)
                                                                                                              gTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGc                                <  2:231235/147‑1 (MQ=255)
                                                                                                               tACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAg                             <  8:282483/149‑1 (MQ=255)
                                                                                                                      cAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGc                                                                                                              <  4:4049/61‑1 (MQ=255)
                                                                                                                      cAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGc                                                                                                              >  3:4049/1‑61 (MQ=255)
                                                                                                                      cAGCACCGACCTGGCCGCCGGGCACCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAccgc                      >  5:5259/1‑149 (MQ=255)
                                                                                                                      cAGCACCGACCTGGCCGCCGGGCACCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTAccgc                      >  5:563/1‑149 (MQ=255)
                                                                                                                      cAGCACCGACCATGCTGCCGGGCGCCAGGTCACGCCTATCAAGCATCCCCTGGCCAAAGGCGACGGTTACTACATGACCTGGCTCAAGACCAGACCAAAAGCGCGGCAGATGAAAAAGCTGCTTGACTTCCTGCTCGGCCAGGTAcagc                      >  2:47790/1‑149 (MQ=255)
                                                                                                                          aCCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCt                  >  6:181646/1‑149 (MQ=255)
                                                                                                                              aCCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCACGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGcc              >  2:315722/1‑149 (MQ=255)
                                                                                                                                   gccgccGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTAc           >  3:235629/1‑147 (MQ=255)
                                                                                                                                   gccgccGGGCACCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTAc           <  6:563/147‑1 (MQ=255)
                                                                                                                                   gccgccGGGCACCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTAc           <  6:5259/147‑1 (MQ=255)
                                                                                                                                    ccgccgGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAAg        <  2:223769/149‑1 (MQ=255)
                                                                                                                                           gcgcCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAAGGACATc  <  5:289312/148‑1 (MQ=255)
                                                                                                                                           gcgcCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAAGGACATc  >  6:289312/1‑148 (MQ=255)
                                                                                                                                           gcgcCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAAGGACATc  >  6:190394/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CGCGAATGGATGCGCTGGACATCAGCATCGACCAGGGGCAGGTGTTCGATACCCTCGACCAGGGCATATCGGCCGCCCAGCAAGGCCTGGGCATTTCAGTGGTCGATCTGGTACTGGCCAGCACCGACCTGGCCGCCGGGCGCCTGGTCACGCCTTTCAAGCATGCCGTGGCCACCGGCGACGGTTACTACATGACCTGGCTCAAGGCCAGCCCAAAAGCGCGGCAGATGCACAAGCTGCGTGACTTCCTGCTCGGCCAGGTACCGCCGCTGGCCTACAAGGACATC  >  NC_002947/2728889‑2729175

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