Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,250,655 C→G 6.9% L14V (CTG→GTG)  ureJ → urease accessory protein UreJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,250,6550CG 6.9% 210.0 / 4.4 69L14V (CTG→GTG) ureJurease accessory protein UreJ
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTACCAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTCGCAACCTGACATTGCCCCTG  >  NC_002947/3250511‑3250801
                                                                                                                                                |                                                                                                                                                   
gcTTACCAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCAGTGc                                                                                                                                                  >  8:279064/1‑148 (MQ=255)
  ttACCAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAg                                                                                                                                               >  3:330280/1‑149 (MQ=255)
  ttACCAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAg                                                                                                                                               >  3:103403/1‑149 (MQ=255)
     ccAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCAGTGCCAGca                                                                                                                                             <  7:279064/148‑2 (MQ=255)
      cAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCAGTGCCAGcatc                                                                                                                                           <  2:37208/149‑4 (MQ=255)
      cAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCtt                                                                                                                                           >  4:348643/1‑149 (MQ=255)
      cAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCtt                                                                                                                                           >  3:319449/1‑149 (MQ=255)
      cAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCtt                                                                                                                                           >  8:156878/1‑149 (MQ=255)
        aaGTGCACCGCATGGCCAGAGCCACAGTGATGCCCATGGCCACGATGACCAGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAATGACTTTGGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTg                                                                                                                                         <  3:361680/149‑1 (MQ=255)
        aaGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCAGTGCCAGcatcaa                                                                                                                                         <  5:67728/149‑6 (MQ=255)
        aaGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTg                                                                                                                                         <  8:64422/149‑1 (MQ=255)
         aGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGc                                                                                                                                        <  4:330280/149‑1 (MQ=255)
          gTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGcc                                                                                                                                       <  5:250237/149‑1 (MQ=255)
           tGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGcc                                                                                                                                       >  1:342605/1‑148 (MQ=255)
             ccccGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTCCGACTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATCAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGACCAt                                                                                                                                    >  8:319125/3‑149 (MQ=255)
                         acagccacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTTAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGc                                                                                                                                        <  8:364606/133‑1 (MQ=255)
                         acagccacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGAc                                                                                                                       <  6:256442/149‑1 (MQ=255)
                         acagccacagTCATGCCCATGGACACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTga                                                                                                                                        >  5:271328/1‑132 (MQ=255)
                         acacccacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTAGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGc                                                                                                                                        >  7:364606/1‑133 (MQ=255)
                         acacccacagTCATGCCAATGGCCACGATCACCCGATCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGc                                                                                                                                        <  6:271328/133‑1 (MQ=255)
                          cagccacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACa                                                                                                                      <  4:103403/149‑1 (MQ=255)
                            gccacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACa                                                                                                                      <  5:199404/147‑1 (MQ=255)
                            gccacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACa                                                                                                                      >  6:199404/1‑147 (MQ=255)
                              cacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGc                                                                                                                             >  8:18066/1‑138 (MQ=255)
                              cacagTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGc                                                                                                                             <  7:18066/138‑1 (MQ=255)
                                      tggctgtGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCt                                                                                                            <  6:178650/141‑1 (MQ=255)
                                      tGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTg                                                                                                           >  5:278932/1‑148 (MQ=255)
                                         ccATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAg                                                                                                        >  6:405953/1‑148 (MQ=255)
                                         ccATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACATTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAg                                                                                                        >  2:243349/1‑148 (MQ=255)
                                              gCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGc                                                                                                  <  2:342605/149‑1 (MQ=255)
                                                cACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCct                                                                                                >  1:259126/1‑149 (MQ=255)
                                                 aCGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCctg                                                                                               >  8:145490/1‑149 (MQ=255)
                                                 aCGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCctg                                                                                               <  6:278932/149‑1 (MQ=255)
                                                  cGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCctgc                                                                                              <  5:405953/149‑1 (MQ=255)
                                                   gATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACTGCctgct                                                                                             <  7:156878/149‑1 (MQ=255)
                                                      cACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCa                                                                                          >  1:112864/1‑149 (MQ=255)
                                                      cACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCAGCAGGCCGTCCtgcagcg                                                                                          >  8:217352/1‑148 (MQ=255)
                                                        cccGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCAGCAGGCCGTCCtgcagcggg                                                                                        <  7:217352/149‑4 (MQ=255)
                                                               gTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACa                                                                                                                      <  7:149472/112‑1 (MQ=255)
                                                               gTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACa                                                                                                                      >  8:149472/1‑112 (MQ=255)
                                                                 aCGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTg                                                                                                                                         >  8:290373/1‑92 (MQ=255)
                                                                 aCGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTg                                                                                                                                         <  7:290373/92‑1 (MQ=255)
                                                                    cctgcctgCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGc                                                                            <  3:348643/149‑1 (MQ=255)
                                                                     ctgcctgcCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCt                                                                           <  2:259126/149‑1 (MQ=255)
                                                                        cctgccCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGa                                                                        <  7:132134/149‑1 (MQ=255)
                                                                        cctgccCATTCCTGAACTGCCATGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCg                                                                         <  6:248716/148‑1 (MQ=255)
                                                                             ccATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGcc                                                                                                 <  8:261617/119‑1 (MQ=255)
                                                                                    tGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGAt                                                            >  7:32663/1‑149 (MQ=255)
                                                                                       aCTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGAtgc                                                          <  7:145490/148‑1 (MQ=255)
                                                                                             cTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGAt                                                                       <  7:15418/129‑1 (MQ=255)
                                                                                             cTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGAt                                                                       >  8:15418/1‑129 (MQ=255)
                                                                                              tGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTggg                                                   >  6:3089/1‑148 (MQ=255)
                                                                                              tGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTggg                                                   >  1:256633/1‑148 (MQ=255)
                                                                                                   cAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATgctgct                                                      <  3:66151/140‑1 (MQ=255)
                                                                                                   cAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATgctgct                                                      >  4:66151/1‑140 (MQ=255)
                                                                                                    aaCCGATGAAAAAGACATTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTAGATCACCTGCAGATGCTGCTGGGTACCGGc                                            <  1:243349/149‑1 (MQ=255)
                                                                                                          tGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTc                                      >  2:383465/1‑149 (MQ=255)
                                                                                                               aaGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTcgccgc                                 <  7:286757/149‑1 (MQ=255)
                                                                                                                    tttCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGa                            >  2:294788/1‑149 (MQ=255)
                                                                                                                       cGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCACGGGCCATGAAAGCAACCCGCTGCAGGACGGCATGATGCACCCGATGACCGGGCTCGATCACATCCTGATGATGATGGTTACCGGCGTCAACGCCGCGTTGAcc                          >  1:83245/1‑149 (MQ=255)
                                                                                                                          cTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTc                      >  8:346468/1‑149 (MQ=255)
                                                                                                                                tCTGTTGACGCTGGCACGTCCAGCCTTTGCCCATGCGGGCCATGACAGCAACCCTCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTCGCAAc                 <  4:285131/148‑1 (MQ=255)
                                                                                                                                   gcaaaggCTGGCAGTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGTTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTCGCAACCTGa             >  8:236608/7‑149 (MQ=255)
                                                                                                                                   gTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCAc                                                                    <  7:238018/94‑1 (MQ=255)
                                                                                                                                   gTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCAc                                                                    >  8:238018/1‑94 (MQ=255)
                                                                                                                                    ttGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATgctgct                                                      >  3:370062/1‑107 (MQ=255)
                                                                                                                                    ttGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATgctgct                                                      <  4:370062/107‑1 (MQ=255)
                                                                                                                                          cTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTCGCAACCTGACATTGcc      >  4:222646/1‑149 (MQ=255)
                                                                                                                                              cACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTCGCAACCTGACATTGCCCCTg  <  1:294788/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCTTACCAAAGTGCACCGCATGGCCACAGCCACAGTCATGCCCATGGCCACGATCACCCGTTCGTACGCCTGCCTGCCCATTCCTGAACTGCCCTGGAGCAACCGATGAAAAAGACTTTCGCCTTGTTTCTGTTGATGCTGGCACTGCCAGCCTTTGCCCATCCGGGCCATGACAGCAACCCGCTGCAGGACGGCCTGCTGCACCCGCTGACCGGGCTCGATCACCTGCTGATGCTGCTGGGTACCGGCGTGCTCGCCGCGTTGACCCGTCGCAACCTGACATTGCCCCTG  >  NC_002947/3250511‑3250801

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: