Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,816,909 A→G 10.9% intergenic (+39/+90) PP_5092 → / ← pilT NLP/P60 protein/twitching motility protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,816,9090AG10.9% 162.1 / 10.1 65intergenic (+39/+90)PP_5092/pilTNLP/P60 protein/twitching motility protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (3/4);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

AAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGT  >  NC_002947/5816763‑5817055
                                                                                                                                                  |                                                                                                                                                  
aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGc                                                                                                                                                 <  2:1580091/150‑1 (MQ=255)
aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGc                                                                                                                                                 >  2:1284950/1‑150 (MQ=255)
  gCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGccc                                                                                                                                               <  1:1616592/150‑1 (MQ=255)
      gTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTc                                                                                                                                           >  1:627656/1‑150 (MQ=255)
      gTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTc                                                                                                                                           >  2:731253/1‑150 (MQ=255)
        gAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCGTGCCCGcgaag                                                                                                                                         <  2:342071/150‑5 (MQ=255)
            aTCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                                     <  1:899534/150‑1 (MQ=255)
                    cTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGcgcg                                                                                                                             >  1:938978/1‑150 (MQ=255)
                       gAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCg                                                                                                                          <  2:627656/150‑1 (MQ=255)
                         gCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGct                                                                                                                        >  1:1279023/1‑150 (MQ=255)
                           ttcagtgatgcagggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGctct                                                                                                                      <  2:256153‑M1/40‑1 (MQ=255)
                           ttcagtgatgcagggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccccTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGctct                                                                                                                      <  2:84254‑M1/40‑1 (MQ=255)
                              ccTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCt                                                                                                                                            >  1:217797/1‑125 (MQ=255)
                              ccTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCt                                                                                                                                            <  2:217797/125‑1 (MQ=255)
                              ccTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTcgc                                                                                                                   <  2:553575/150‑1 (MQ=255)
                               cTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGctc                                                                                                                       >  2:870258/1‑145 (MQ=255)
                               cTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGctc                                                                                                                       <  1:870258/145‑1 (MQ=255)
                                    tGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAg                                                                                                             <  2:1279023/150‑1 (MQ=255)
                                      agggcggggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCt                                                                                                           <  2:634275‑M1/51‑1 (MQ=255)
                                        aGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcc                                                                                                         >  2:276792/1‑149 (MQ=255)
                                           aGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAt                                                                                                      <  1:1387583/150‑1 (MQ=255)
                                            gCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCCGGCACGCCCGCTCCCACACGTAGCGCGTCGCTCTCGCGCGCAGCCcacctc                                                                                                     >  2:617312/1‑147 (MQ=255)
                                             ggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATAc                                                                                                    >  2:1254116‑M1/93‑150 (MQ=255)
                                             cACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATAc                                                                                                    <  2:954123/150‑1 (MQ=255)
                                                  cggtggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccccTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg                                                                                               >  1:561681‑M1/88‑150 (MQ=255)
                                                  cGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg                                                                                               <  2:702973/150‑1 (MQ=255)
                                                   gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              <  1:731253/150‑1 (MQ=255)
                                                   gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              <  1:394217/150‑1 (MQ=255)
                                                   ggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              >  1:256153‑M1/87‑150 (MQ=255)
                                                   ggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              >  1:634275‑M1/87‑150 (MQ=255)
                                                      ggggatgctgacgatggagggtgcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccccTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAg                                                                                           >  1:84254‑M1/84‑150 (MQ=255)
                                                                      tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGGACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCcgcgca                                                                           >  1:531214/1‑148 (MQ=255)
                                                                      tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCACGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                           >  1:531707/1‑150 (MQ=255)
                                                                       gcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGaaa                                                                          >  1:1467240/1‑150 (MQ=255)
                                                                          gAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAggg                                                                       >  1:803615/1‑150 (MQ=255)
                                                                          gAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCAAGCCCGATCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAggg                                                                       >  2:545636/1‑150 (MQ=255)
                                                                                gcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc                                                                                    >  1:603727/1‑131 (MQ=255)
                                                                                gcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc                                                                                    <  2:603727/131‑1 (MQ=255)
                                                                                gcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTg                                                                 >  2:1415074/1‑150 (MQ=255)
                                                                                gcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCCCGGGCCCGCCCGCTCCCACCCCGAGCGCGCCGCCCTCACCCGCCGCTCCCACCCCCGTGGGCACGAGGACGCCCGCGAACggcactatg                                                                 >  2:1582939/1‑143 (MQ=255)
                                                                                   gcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAg                                                              >  1:1187754/1‑150 (MQ=255)
                                                                                   gcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAg                                                              <  1:1284950/150‑1 (MQ=255)
                                                                                    cagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAgg                                                             <  2:153657/150‑1 (MQ=255)
                                                                                     agcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGt                                                            <  2:1187754/150‑1 (MQ=255)
                                                                                          gagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              >  2:288097‑M1/48‑111 (MQ=255)
                                                                                          gagggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              <  1:288097‑M1/64‑1 (MQ=255)
                                                                                             ggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcac                                                                                                        >  2:92708‑M1/45‑98 (MQ=255)
                                                                                             ggagcgagggttgatcaaaggcttgtagtacctgcaaaggccctTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTcac                                                                                                        <  1:92708‑M1/54‑1 (MQ=255)
                                                                                                 ccAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCtt                                                <  2:360220/150‑1 (MQ=255)
                                                                                                 ccAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCtt                                                <  1:1415074/150‑1 (MQ=255)
                                                                                                  cAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCGCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCttt                                               <  1:276792/150‑1 (MQ=255)
                                                                                                   aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg                                                                        >  1:1173740/1‑124 (MQ=255)
                                                                                                   aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAgg                                                                        <  2:1173740/124‑1 (MQ=255)
                                                                                                          aaGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAcc                                       <  2:803615/150‑1 (MQ=255)
                                                                                                             cTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTc                                    <  2:1467240/150‑1 (MQ=255)
                                                                                                               gATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGc                                  <  1:545636/150‑1 (MQ=255)
                                                                                                                   ttgtagtacctgcaaaggccctTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg                                                                                               <  2:718672‑M1/63‑1 (MQ=255)
                                                                                                                   ttgtagtacctgcaaaggccctTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg                                                                                               >  1:718672‑M1/23‑85 (MQ=255)
                                                                                                                   ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              >  2:655154/1‑86 (MQ=255)
                                                                                                                   ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              <  1:655154/86‑1 (MQ=255)
                                                                                                                        tATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCTCGCGAAAGGGCCTTTGCAGGTACGACAAGCCTTGGATCAACCCTCGCTCCCTCAGc                         >  2:514866/1‑150 (MQ=255)
                                                                                                                                 ccGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGTAGCGCGTCGCTCCAGCGCGCAGCCCACACACCCGTGGGAGCGGGCGCGCCCACGAAAGGGCCTTTGCAGGTACCACAAGCCTTTGATCCACCCTCGCCCCCTCAGCCCCCGCAAc                >  2:137602/1‑150 (MQ=255)
                                                                                                                                      cctTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              >  2:671097‑M1/4‑67 (MQ=255)
                                                                                                                                      cctTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                              <  1:671097‑M1/64‑1 (MQ=255)
                                                                                                                                               GGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGt  <  1:917874‑M1/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGT  >  NC_002947/5816763‑5817055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: