Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 889,441 A→C 10.0% intergenic (+178/‑16) PP_0769 → / → PP_0770 sensor histidine kinase/PemI‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947889,4410AC10.0% 159.4 / 9.3 60intergenic (+178/‑16)PP_0769/PP_0770sensor histidine kinase/PemI‑like protein
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base C (2/4);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GCAGTTGGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAA  >  NC_002947/889294‑889588
                                                                                                                                                   |                                                                                                                                                   
gCAGTTGGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTAc                                                                                                                                                    >  6:306002/1‑149 (MQ=255)
   gTTGGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTAcc                                                                                                                                                   >  5:123559/1‑147 (MQ=255)
      ggAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTggg                                                                                                                                              >  7:148209/1‑149 (MQ=255)
      ggAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTgg                                                                                                                                               >  3:210852/1‑148 (MQ=255)
       gagaGGGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGga                                                                                                                                             <  1:91532/149‑1 (MQ=255)
        tgagTTGCTTCGCATGGTGATGTAGGTCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTGAGGGCAGATGACAGTGTATATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGgag                                                                                                                                            <  6:123559/148‑1 (MQ=255)
          agTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGAt                                                                                                                                          <  8:137791/149‑1 (MQ=255)
              gCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTc                                                                                                                                      <  4:193405/149‑1 (MQ=255)
              gCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTCGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTc                                                                                                                                      <  1:60310/149‑1 (MQ=255)
               cTTCCCATGGTGATGTCTGGCATTTCCTTGTTGTGGTTGCAGGGCTTTTGTTCTGTCAGGGCTGTTTAAATTGTTTATACCATGTTTCTACATTTTGCTTTCAATCATTGGTTGGGTTTCTACTAAGTATCTCCCAAGGTAGATAcata                                                                                                                                     <  7:255685/149‑4 (MQ=255)
               cTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAcata                                                                                                                                     <  3:284814/149‑4 (MQ=255)
                  cccAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATg                                                                                                                                   <  1:102390/148‑1 (MQ=255)
                   ccAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAcatagtag                                                                                                                                 <  6:30360/149‑8 (MQ=255)
                      aGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGAt                                                                                                                              <  5:306002/149‑1 (MQ=255)
                       ggTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCATGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTAGCTTGGGAGATACTCATGCAGATc                                                                                                                             <  8:34295/149‑1 (MQ=255)
                            tGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGa                                                                                                                         >  4:186640/1‑148 (MQ=255)
                               aGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTc                                                                                                                      >  8:146995/1‑148 (MQ=255)
                                gggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAca                                                                                                                                       <  4:380658/130‑2 (MQ=255)
                                gggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAca                                                                                                                                       >  3:380658/1‑129 (MQ=255)
                                gggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCaa                                                                                                                    >  3:183462/1‑149 (MQ=255)
                                 ggCATTTCCTTCTTGAGGTTGCATGTCTTTTGTTCGTTCAGGGCATATGACAGTTTAGATACGACGTTTCTACATTTCGCTTGGAATCATTTGGAGGGTTTCTACTATGTATCTACCTTTGGAGATACTCATGCAGATAAAGATTaaa                                                                                                                    <  4:210852/148‑1 (MQ=255)
                                 ggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCaa                                                                                                                    >  5:135995/1‑148 (MQ=255)
                                 ggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCaa                                                                                                                    >  8:178370/1‑148 (MQ=255)
                                    aTTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGt                                                                                                                >  5:414762/1‑149 (MQ=255)
                                     tttCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTg                                                                                                               <  5:292875/149‑1 (MQ=255)
                                        ccTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTggg                                                                                                             <  6:414762/148‑1 (MQ=255)
                                             ttGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCaa                                                                                                         >  4:25418/1‑147 (MQ=255)
                                                       aGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGAAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGAt                                                                                             >  8:25531/1‑149 (MQ=255)
                                                               tGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTg                                                                                      <  3:186640/148‑1 (MQ=255)
                                                                 ttCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGccc                                                                                   >  3:133131/1‑149 (MQ=255)
                                                                          ggCAGATGACAGTGTAGATAAGACTTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTa                                                                          >  8:45303/1‑149 (MQ=255)
                                                                            cAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTATCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTAc                                                                         >  2:198329/1‑148 (MQ=255)
                                                                              gATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCa                                                                      <  4:133131/149‑1 (MQ=255)
                                                                                tGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCaa                                                                     <  7:146995/148‑1 (MQ=255)
                                                                                      tgtAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATg                                                              >  1:132917/1‑149 (MQ=255)
                                                                                      tgtAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATg                                                              <  4:190797/149‑1 (MQ=255)
                                                                                        tAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATgc                                                             <  6:117438/148‑1 (MQ=255)
                                                                                           aTACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCt                                                         >  6:349664/1‑149 (MQ=255)
                                                                                            tACGACGTTTCTACATTTCGCTTGCAATTATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTc                                                        >  4:255542/1‑149 (MQ=255)
                                                                                              cgaagTTTCTACATTTCGCTTGCAATTATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCgg                                                      <  3:255542/149‑1 (MQ=255)
                                                                                               gacgTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTg                                                                                      <  5:100611/116‑1 (MQ=255)
                                                                                               gacgTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTg                                                                                      >  6:100611/1‑116 (MQ=255)
                                                                                                 cgTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGgtgt                                                   <  8:148209/149‑1 (MQ=255)
                                                                                                  gTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTc                                                  <  7:25531/149‑1 (MQ=255)
                                                                                                   tttCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCg                                                 >  6:402827/1‑149 (MQ=255)
                                                                                                        aCATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCACGCAGATGCGCCTCGGTGTCGGCTcc                                            >  7:380669/1‑149 (MQ=255)
                                                                                                          aTTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCAc                                          >  7:74215/1‑149 (MQ=255)
                                                                                                              cGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCt                                       >  5:394827/1‑148 (MQ=255)
                                                                                                                cTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTc                                                        >  6:80480/1‑129 (MQ=255)
                                                                                                                cTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTc                                                        <  5:80480/129‑1 (MQ=255)
                                                                                                                    cAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCtt                                <  2:132917/149‑1 (MQ=255)
                                                                                                                      aTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGa                              >  4:131164/1‑149 (MQ=255)
                                                                                                                         aTTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGacac                           >  8:160877/1‑149 (MQ=255)
                                                                                                                         aTTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGaca                            >  7:101881/1‑148 (MQ=255)
                                                                                                                                  ttttCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTga                  <  7:375792/148‑1 (MQ=255)
                                                                                                                                  gTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTga                  <  2:175536/149‑1 (MQ=255)
                                                                                                                                  gTTTCTACTATGTATCTACCTTGGCAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCCCCGCAGTACTCAAGCAGATGCGACTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACAGGTg                   >  3:312933/1‑148 (MQ=255)
                                                                                                                                       tactaTGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGa             <  1:198329/149‑1 (MQ=255)
                                                                                                                                         atgagTATCTCCCAAGGTAGATACTCATGCAGATCAAGATTCAACAGTGGTGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCACCACCCTGAGCCTTGACACAACGGGTGAGACGATg           >  8:161818/5‑149 (MQ=255)
                                                                                                                                                  tACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCaaa  >  4:249584/1‑149 (MQ=255)
                                                                                                                                                   aCCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACg                       <  4:373794/127‑1 (MQ=255)
                                                                                                                                                   aCCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACg                       >  3:373794/1‑127 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAGTTGGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAA  >  NC_002947/889294‑889588

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