Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 889,444 T→A 9.8% intergenic (+181/‑13) 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,4440TA 9.8% 163.7 / 8.2 61intergenic (+181/‑13)PP_0769/PP_0770sensor histidine kinase/PemI‑like protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base A (2/4);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCG  >  NC_002947/889300‑889592
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                        >  8:178370/1‑148 (MQ=255)
                           ggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCaa                                                                                                                        >  5:135995/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)
                                                                                                                                              ccTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAAcc    >  1:114462/1‑149 (MQ=255)
                                                                                                                                                ttGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCg  >  1:398205/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCG  >  NC_002947/889300‑889592

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: