Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,250,260 A→G 7.5% E56G (GAA→GGA)  ureE → urease accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,250,2600AG 7.5% 146.4 / 2.8 54E56G (GAA→GGA) ureEurease accessory protein
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base G (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTATATCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTGCAGATCGAACCC  >  NC_002947/3250113‑3250408
                                                                                                                                                     |                                                                                                                                                      
gTATATCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CTGAGGTGaa                                                                                                                                                       <  6:209386/149‑1 (MQ=255)
  atatCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAAct                                                                                                                                                     >  1:180932/1‑149 (MQ=255)
    atCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGCCGATGGCCGCGAAGCCGGGCTGATGCTCGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAActg‑‑c                                                                                                                                                 >  6:107920/1‑149 (MQ=255)
     tCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAActg‑‑c                                                                                                                                                 <  4:271317/148‑1 (MQ=255)
     tCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAActg‑‑c                                                                                                                                                 >  3:271317/1‑148 (MQ=255)
     tCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAACCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAActg‑‑ct                                                                                                                                                >  6:330348/1‑149 (MQ=255)
        gAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAACCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGc                                                                                                                                             <  5:330348/149‑1 (MQ=255)
             cGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATg                                                                                                                                        >  1:183324/1‑149 (MQ=255)
                                  cACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATcc                                                                                                                   >  3:104341/1‑149 (MQ=255)
                                      gTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGt                                                                                                               >  1:289278/1‑149 (MQ=255)
                                      gTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGAAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAg                                                                                                                          >  8:257666/1‑138 (MQ=255)
                                      gTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGAAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAg                                                                                                                          <  7:257666/138‑1 (MQ=255)
                                         aCCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgct                                                                                                            >  7:217354/1‑149 (MQ=255)
                                         aCCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgct                                                                                                            >  6:233443/1‑149 (MQ=255)
                                          cccTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgctg                                                                                                           <  1:289152/149‑1 (MQ=255)
                                           ccTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGGACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgctgg                                                                                                          >  8:151127/1‑149 (MQ=255)
                                           ccTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGGACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgctgg                                                                                                          >  2:236228/1‑149 (MQ=255)
                                           ccTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgctgg                                                                                                          <  2:43972/149‑1 (MQ=255)
                                               gACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTgctggctg                                                                                                       <  6:193834/148‑1 (MQ=255)
                                                aCGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTgc                                                                                                      <  4:44374/148‑1 (MQ=255)
                                                 cGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTgcgc                                                                                                    <  5:233443/149‑1 (MQ=255)
                                                     gACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGccc                                                                                                  >  3:100742/1‑147 (MQ=255)
                                                      aCAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAg                                                                                               >  6:189150/1‑149 (MQ=255)
                                                        aGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGc                                                                                             >  2:21024/1‑149 (MQ=255)
                                                         gCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCg                                                                                            <  2:180932/149‑1 (MQ=255)
                                                           cGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGGACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGgt                                                                                          <  7:151127/149‑1 (MQ=255)
                                                           cGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGGACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGgt                                                                                          <  1:236228/149‑1 (MQ=255)
                                                                  agagCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGTGGAGGGAGGTGAACTGAAGTGGCCCATGCCGAGGGCCCCCCCCTGATCCCCGTGCTGGCTGCGCCCGAGGCGGTTTc                                                                                       >  1:267672/1‑149 (MQ=255)
                                                                      ccgccTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTgc                                                                               >  1:179645/1‑149 (MQ=255)
                                                                        gccTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTgcgc                                                                             <  5:189150/149‑1 (MQ=255)
                                                                             cGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCtgtg                                                                         <  2:289278/148‑1 (MQ=255)
                                                                                gTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGCACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCg                                                                      >  7:291017/1‑148 (MQ=255)
                                                                                     ccTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCcac                                                                <  8:217354/149‑1 (MQ=255)
                                                                                     ccTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCcac                                                                >  5:58088/1‑149 (MQ=255)
                                                                                         gACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCt                                                            <  4:104341/149‑1 (MQ=255)
                                                                                           cGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGt                                                          <  6:271804/149‑1 (MQ=255)
                                                                                             aTGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTg                                                        <  3:37765/149‑1 (MQ=255)
                                                                                             aTGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTg                                                        <  3:157645/149‑1 (MQ=255)
                                                                                              tGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTgg                                                       >  7:117852/1‑149 (MQ=255)
                                                                                                   gcgAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCAcgc                                                  <  4:100742/149‑1 (MQ=255)
                                                                                                      aaGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCTGTGTCGACGGTGCGCTGTGCTGATCCACACCTGTTGGCACGcgcc                                               <  2:179645/149‑1 (MQ=255)
                                                                                                           gggCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTa                                          >  6:80766/1‑149 (MQ=255)
                                                                                                              cTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTacca                                       >  5:303194/1‑149 (MQ=255)
                                                                                                                 aTGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCt                                    >  1:327722/1‑149 (MQ=255)
                                                                                                                 aTGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCt                                    <  2:267672/149‑1 (MQ=255)
                                                                                                                    cTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTggg                                 <  5:80766/149‑1 (MQ=255)
                                                                                                                    cTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTgg                                  <  4:159884/148‑1 (MQ=255)
                                                                                                                       gAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGt                                <  6:58088/147‑1 (MQ=255)
                                                                                                                               ccACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTa                                          <  7:56782/129‑1 (MQ=255)
                                                                                                                               ccACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTa                                          >  8:56782/1‑129 (MQ=255)
                                                                                                                                      ctgcG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGc                <  1:80611/148‑1 (MQ=255)
                                                                                                                                          g‑‑cgGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTg              >  4:103661/1‑146 (MQ=255)
                                                                                                                                               gAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTGCAGATc        <  1:21024/149‑1 (MQ=255)
                                                                                                                                                 ggTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTGCAGATCg       >  1:233360/1‑148 (MQ=255)
                                                                                                                                                   tGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTGCAGATCGAAc    <  8:291017/149‑1 (MQ=255)
                                                                                                                                                     aaCTG‑‑CTGGCCGATGCCGAGGGCACCAAGATGATCCGCGTGATGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCCCGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTGCAGATCGAAccc  >  3:193641/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTATATCCGAACCCGGTACGCGGGTTGTCAGTGGCACTGTCACCCTTGACGTGGACAGCCGTATCAAGAGCCGCCTGCGAGTAACCCTGGACGATGGCCGCGAAGCCGGGCTGATGCTGGAGCGCGGCCACCTGCTGCG‑‑CGGAGGTGAACTG‑‑CTGGCCGATGCCGAGGGCACCCAGCTGATCCGCGTGCTGGCTGCGCCCGAGGCGGTGTCGACGGTGCGCTGTGCCGATCCACACCTGTTGGCACGCGCCGCCTACCACCTGGGTAACCGCCATGTGCCGCTGCAGATCGAACCC  >  NC_002947/3250113‑3250408

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