Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,119,196 G→C 8.5% E324Q (GAG→CAG)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,119,1960GC 8.5% 178.8 / 7.4 59E324Q (GAG→CAG) PP_4505sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base C (2/3);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

GCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCA  >  NC_002947/5119050‑5119343
                                                                                                                                                  |                                                                                                                                                   
gCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAg                                                                                                                                                   <  6:278900/149‑1 (MQ=255)
         aCCTGATCCGCGAAATCCTCGCCGTGGCCCGGGTGGATGCCTAACAAGCCCATGCCGACCCGGTCTACCTCAATCCCCTACTCGGCAGCGGGCGAAAGGACGCCCTGCTGAGCGCGCCGGATCAGTAGGTGCGCCTCCAGGCGCacgtc                                                                                                                                          <  6:38256/149‑5 (MQ=255)
           cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAgg                                                                                                                                        >  1:311066/1‑149 (MQ=255)
                   cGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCt                                                                                                                                >  2:76935/1‑149 (MQ=255)
                   cGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCACGTTTGAGCt                                                                                                                                >  5:323386/1‑149 (MQ=255)
                     aaaTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTg                                                                                                                              >  6:190361/1‑149 (MQ=255)
                        tCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAg                                                                                                                           <  3:304740/149‑1 (MQ=255)
                          cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAgg                                                                                                                          <  6:149073/148‑1 (MQ=255)
                          cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGCACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAggg                                                                                                                         >  7:128744/1‑149 (MQ=255)
                             gCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                       <  8:238486/148‑1 (MQ=255)
                                   gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGtt                                                                                                                                      <  7:77576/127‑1 (MQ=255)
                                   gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGtt                                                                                                                                      >  8:77576/1‑127 (MQ=255)
                                    cccGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                               >  4:176854/1‑149 (MQ=255)
                                            gATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                   >  6:246777/1‑121 (MQ=255)
                                                    aCAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCg                                                                                               >  7:18640/1‑149 (MQ=255)
                                                      aaGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGtt                                                                                             >  8:88638/1‑149 (MQ=255)
                                                          ccATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGAc                                                                                          >  2:264124/1‑148 (MQ=255)
                                                            aTGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAAc                                                                                       >  4:286337/1‑149 (MQ=255)
                                                                cGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgc                                                                                   >  6:46339/1‑149 (MQ=255)
                                                                 gAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgct                                                                                  >  3:133538/1‑149 (MQ=255)
                                                                       gTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCaa                                                                            >  8:65490/1‑149 (MQ=255)
                                                                       gTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCaa                                                                            >  4:266719/1‑149 (MQ=255)
                                                                             cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                      <  2:311066/149‑1 (MQ=255)
                                                                                  tGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCt                                                                 <  7:88638/149‑1 (MQ=255)
                                                                                   gCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGATCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCtt                                                                >  2:70931/1‑149 (MQ=255)
                                                                                       tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAc                                                            >  5:326225/1‑149 (MQ=255)
                                                                                       tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAc                                                            >  3:85170/1‑149 (MQ=255)
                                                                                          tCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcg                                                         >  7:146026/1‑149 (MQ=255)
                                                                                          tCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGA‑AACCTGCTGCGCAATGCCCTGCGCTTCAACCcgg                                                        >  7:244892/1‑149 (MQ=255)
                                                                                                gCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGc                                                   <  8:146026/149‑1 (MQ=255)
                                                                                                      gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                  >  8:152024/1‑144 (MQ=255)
                                                                                                      gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                  <  7:152024/144‑1 (MQ=255)
                                                                                                      gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCg                                             >  8:23721/1‑149 (MQ=255)
                                                                                                           gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                        <  7:65490/149‑1 (MQ=255)
                                                                                                            aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                       <  8:128744/149‑1 (MQ=255)
                                                                                                            aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                       <  8:18640/149‑1 (MQ=255)
                                                                                                            aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                       <  3:286337/149‑1 (MQ=255)
                                                                                                            aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                       <  1:76935/149‑1 (MQ=255)
                                                                                                             cGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAgg                                      <  4:21210/149‑1 (MQ=255)
                                                                                                              gCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGt                                     >  1:185043/1‑149 (MQ=255)
                                                                                                                ccTGCTGAGCGCGCCGGATCAGGACGTGCGCATCGAGGAGCAGCGACGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCGGTTGACAACATGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGGCAGCCGGTCGAGGt                                     >  7:153000/1‑147 (MQ=255)
                                                                                                                  tgctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAgc                                 >  6:74902/1‑149 (MQ=255)
                                                                                                                   gctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGc                                               <  6:300749/134‑1 (MQ=255)
                                                                                                                   gctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGc                                               >  5:300749/1‑134 (MQ=255)
                                                                                                                       aGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGcgcc                              <  6:7505/147‑1 (MQ=255)
                                                                                                                        gcgcgcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGA‑AACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCcgcgc                          <  8:244892/149‑1 (MQ=255)
                                                                                                                            gcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCcgcgcgcg                       <  5:46339/149‑1 (MQ=255)
                                                                                                                             ccGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGa                      <  5:290285/149‑1 (MQ=255)
                                                                                                                              cGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAg                     >  1:196618/1‑149 (MQ=255)
                                                                                                                              cGGATCAGGACGTGCGCATGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAg                     >  3:141374/1‑149 (MQ=255)
                                                                                                                                gATCAGGACGTGCGCATGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCa                   <  4:141374/149‑1 (MQ=255)
                                                                                                                                 aTCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAg                  <  7:23721/149‑1 (MQ=255)
                                                                                                                                 aTCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAg                  <  2:196618/149‑1 (MQ=255)
                                                                                                                                     ctggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                  <  6:214931/108‑1 (MQ=255)
                                                                                                                                     ctggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGcc                                                  >  5:214931/6‑113 (MQ=255)
                                                                                                                                            cgcCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCg                                         <  7:35772/115‑1 (MQ=255)
                                                                                                                                            cgcCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCg                                         >  8:35772/1‑115 (MQ=255)
                                                                                                                                             gcCTGGAGGCTCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTg      <  1:70931/149‑1 (MQ=255)
                                                                                                                                               cTGGAGGCGCAGCCAGGTTCGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGAt    >  2:295954/1‑149 (MQ=255)
                                                                                                                                                  gAGGCGCAGCCAGGTTTGAGCTTGCAGCGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCa  <  5:190361/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCA  >  NC_002947/5119050‑5119343

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