Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 889,448 G→T 9.8% intergenic (+185/‑9) 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,4480GT 9.8% 168.4 / 7.3 58intergenic (+185/‑9)PP_0769/PP_0770sensor histidine kinase/PemI‑like protein
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base T (2/4);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGTCAG  >  NC_002947/889300‑889596
                                                                                                                                                    |                                                                                                                                                    
ggAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTggg                                                                                                                                                      >  7:148209/1‑149 (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                                                                                                                                               >  3:380658/1‑129 (MQ=255)
                          gggCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAca                                                                                                                                               <  4:380658/130‑2 (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                                                                      <  4:190797/149‑1 (MQ=255)
                                                                                tgtAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATg                                                                      >  1:132917/1‑149 (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                                                                                              >  6:100611/1‑116 (MQ=255)
                                                                                         gacgTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTg                                                                                              <  5:100611/116‑1 (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                                                                <  5:80480/129‑1 (MQ=255)
                                                                                                          cTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTc                                                                >  6:80480/1‑129 (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)
                                                                                                                                                  gggAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGt     <  3:249584/148‑1 (MQ=255)
                                                                                                                                                   ggAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGt     >  7:232919/1‑147 (MQ=255)
                                                                                                                                                   ggAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGTc    <  3:131164/148‑1 (MQ=255)
                                                                                                                                                   ggAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGTc    <  6:230406/148‑1 (MQ=255)
                                                                                                                                                    gagaTACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGTCAg  <  5:402827/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAGAGTGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGGTGCTCAAACCCGTCAG  >  NC_002947/889300‑889596

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: