Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,770,214 T→C 6.8% intergenic (‑80/+29) PP_2422 ← / ← lpxOA hypothetical protein/Fe(2+)/alpha‑ketoglutarate‑dependent dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,770,2140TC 6.8% 165.7 / 2.2 60intergenic (‑80/+29)PP_2422/lpxOAhypothetical protein/Fe(2+)/alpha‑ketoglutarate‑dependent dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base C (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ATCGCCGATGCCGGCCATTAACTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCTTGTAGATGCGGG  >  NC_002947/2770069‑2770361
                                                                                                                                                 |                                                                                                                                                    
atagCCGATGCCGGCCATTAACTGCTTACGGTAGGCGTCCAAGTCCATAATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGATATTCGAAGCGTGTCACGTGAGAGGCGGTACTTTGGTTGGTAAACGCAGGTCAc                                                                                                                                                   >  6:47950/4‑149 (MQ=255)
aTCGCCGATGCCGGCCATTAACTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGATTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCAc                                                                                                                                                   >  6:42749/1‑149 (MQ=255)
           cGGCCATTAACTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTgg                                                                                                                                       >  1:265651/1‑149 (MQ=255)
                aTTAACTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGc                                                                                                                                  <  5:320304/149‑1 (MQ=255)
                aTTAACTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGc                                                                                                                                  >  5:76957/1‑149 (MQ=255)
                    aCTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTg                                                                                                                              <  4:204281/149‑1 (MQ=255)
                     cTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGGAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGa                                                                                                                             <  1:237707/149‑1 (MQ=255)
                        cTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTTAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGcc                                                                                                                          <  5:27089/149‑1 (MQ=255)
                        cTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGcc                                                                                                                          >  6:27089/1‑149 (MQ=255)
                        cTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGc                                                                                                                           <  8:120473/148‑1 (MQ=255)
                         ttACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGt                                                                                                                                >  6:200436/1‑142 (MQ=255)
                         ttACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGt                                                                                                                                <  5:200436/142‑1 (MQ=255)
                               tAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGAt                                                                                                                   >  8:158118/1‑149 (MQ=255)
                                aGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATg                                                                                                                  >  1:196645/1‑149 (MQ=255)
                                 ggcggcCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGt                                                                                                                 <  3:242285/149‑1 (MQ=255)
                                 ggCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGCCACCGTACGAGTGACCTGCGt                                                                                                                                >  6:107651/1‑134 (MQ=255)
                                 ggCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGCCACCGTACGAGTGACCTGCGt                                                                                                                                <  5:107651/134‑1 (MQ=255)
                                  gCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTa                                                                                                                <  5:42749/149‑1 (MQ=255)
                                  gCGTCCCAGTACATCATGCTCATCGTCTTCACCTTGGGGGGTGCGAGTGGAGATTTCAGTCTAGTAGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCTTTGAGCCTTCAGATGTa                                                                                                                <  5:47950/149‑1 (MQ=255)
                                     tCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGAt                                                                                                             <  7:158118/149‑1 (MQ=255)
                                       ccAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATg                                                                                                            >  7:308683/1‑148 (MQ=255)
                                         aGTCTATCATGCTCATCGTCTTCATCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAac                                                                                                         <  3:333095/149‑1 (MQ=255)
                                           tCCTTCATGGTCATCGTCTTCACCTTGGGGGGTGCGAGTGGATCTTTCAGGCTATTCTAAGCGTGTCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAacac                                                                                                       <  6:17170/149‑1 (MQ=255)
                                           tCCATCTTGCTCCTCGTCTTCTCCTTGGGGGTTGCGCGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAacac                                                                                                       <  8:308683/149‑1 (MQ=255)
                                             catcatGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAacac                                                                                                       >  5:176092/1‑147 (MQ=255)
                                                 atGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCAc                                                                                                 >  7:209214/1‑149 (MQ=255)
                                                       aTCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAAcgc                                                                                           >  3:130732/1‑149 (MQ=255)
                                                        tCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAAcgcc                                                                                          >  8:71436/1‑149 (MQ=255)
                                                         cGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAAcgccg                                                                                         >  3:242703/1‑149 (MQ=255)
                                                                     gggggTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGcc                                                                              >  6:55028/1‑148 (MQ=255)
                                                                        ggTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCAtt                                                                          >  2:236807/1‑149 (MQ=255)
                                                                         gTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCAttt                                                                         >  3:183837/1‑149 (MQ=255)
                                                                            gCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCAttttt                                                                       <  7:284384/148‑1 (MQ=255)
                                                                              gAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCa                                                                    <  1:236807/149‑1 (MQ=255)
                                                                              gAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCa                                                                    <  7:71436/149‑1 (MQ=255)
                                                                                gTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTc                                                                                                                       <  8:133068/96‑1 (MQ=255)
                                                                                gTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTc                                                                                                                       >  7:133068/1‑96 (MQ=255)
                                                                                gTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAgg                                                                  <  2:196645/149‑1 (MQ=255)
                                                                                 tGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAgg                                                                  <  8:209214/148‑1 (MQ=255)
                                                                                    aGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGgtagt                                                              <  4:242703/149‑1 (MQ=255)
                                                                                     gCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGgtag                                                               >  5:103908/1‑147 (MQ=255)
                                                                                         tCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAgcg                                                          >  8:292805/1‑148 (MQ=255)
                                                                                               tAGTCGAAGCGTTGCACGTGCGAGGCGGTACTTTGGTGGGTAAAAGCAGGACACACCTACGGTGGCCTGCGTATAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCAATTTTCCAGGTAGAAGCGCATGCgg                                                   >  7:204038/1‑149 (MQ=255)
                                                                                               tAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGtagtag                                                            <  3:308186/140‑1 (MQ=255)
                                                                                               tAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGtagtag                                                            >  4:308186/1‑140 (MQ=255)
                                                                                                    gAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTgc                                               <  5:55028/148‑1 (MQ=255)
                                                                                                       gCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCAACCGCATGCGcta                                           <  6:76957/149‑2 (MQ=255)
                                                                                                       gCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTAGTACGGGGGCCTGCGTTGAGCCTTAAGATGTAGATGAACACCAGCACCAACGCCGGGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGctt                                           <  8:204038/149‑1 (MQ=255)
                                                                                                         gTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGcttct                                         <  4:183837/149‑1 (MQ=255)
                                                                                                          tGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGcttctt                                        <  4:130732/149‑1 (MQ=255)
                                                                                                            gCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCttcttc                                       <  6:103908/148‑1 (MQ=255)
                                                                                                                                tggtggGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAgc                                                           >  3:45934/1‑108 (MQ=255)
                                                                                                                                tggtggGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAgc                                                           <  4:45934/108‑1 (MQ=255)
                                                                                                                                    gggTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCt              >  6:53712/1‑149 (MQ=255)
                                                                                                                                       tcaaCGCAGGCCACCGTACGAGTGACCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTACCACGCTCGCGGATCtt             <  3:225578/145‑1 (MQ=255)
                                                                                                                                       tcaaCGCAGGCCACCGTACGAGTGACCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTACCACGCTCGCGGATCtt             >  4:225578/3‑147 (MQ=255)
                                                                                                                                       tAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCTTGt           <  2:265651/149‑1 (MQ=255)
                                                                                                                                         aaCGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCTTGTa          <  7:292805/148‑1 (MQ=255)
                                                                                                                                               ggTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCTTGTAGATGCgg   <  5:306433/149‑1 (MQ=255)
                                                                                                                                                gTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCTTGTAGATGCggg  >  7:190025/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATCGCCGATGCCGGCCATTAACTGCTTACGGTAGGCGTCCCAGTCCATCATGCTCATCGTCTTCACCTTGGGGGTTGCGAGTGGAGCTTTCAGGCTAGTCGAAGCGTGGCACGTGCGAGGCGGTACTTTGGTGGGTAAACGCAGGTCACTCGTACGGTGGCCTGCGTTGAGCCTTCAGATGTAGATGAACACCAGCACCAACGCCGCGACCACTGCCCATTTTTCCAGGTAGTAGCGCATGCGGTTGCGCTTCTTCAGCGCCTTGCCACGCTCGCGGATCTTGTAGATGCGGG  >  NC_002947/2770069‑2770361

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: