Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,512,951 C→G 9.2% A180G (GCC→GGC)  urtD → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,512,9510CG 9.2% 176.6 / 6.5 61A180G (GCC→GGC) urtDABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base G (3/3);  total (32/29)
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.90e-01

CTGACCGTGTTCGAGAACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGAATTTACCGCCGAGCTGTTC  >  NC_002947/5512815‑5513099
                                                                                                                                        |                                                                                                                                                    
cTGACCGTGTTCGAGAACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGcc                                                                                                                                         >  2:550724/1‑150 (MQ=255)
      gtgtTCGAGAACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTg                                                                                                                                   <  2:1088309/150‑1 (MQ=255)
        gtTCGAGAACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCt                                                                                                                                 <  1:1226430/150‑1 (MQ=255)
        gtTCGAGAACCTGGAGATGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGTCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCt                                                                                                                                 >  1:102174/1‑150 (MQ=255)
         ttCGAGAACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTg                                                                                                                                >  1:130903/1‑150 (MQ=255)
            gagaACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCg                                                                                                                             >  1:987727/1‑150 (MQ=255)
            gagaACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCg                                                                                                                             >  2:911759/1‑150 (MQ=255)
            gagaACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCg                                                                                                                             >  1:975106/1‑150 (MQ=255)
                 ccTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATgg                                                                                                                        >  1:1266556/1‑150 (MQ=255)
                 ccTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATgg                                                                                                                        <  2:130903/150‑1 (MQ=255)
                  cTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATggg                                                                                                                       <  1:658588/150‑1 (MQ=255)
                  cTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATggg                                                                                                                       <  2:855285/150‑1 (MQ=255)
                       gatGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGaa                                                                                                                  <  2:102174/148‑1 (MQ=255)
                          ggCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGTTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCa                                                                                                               <  1:976269/150‑1 (MQ=255)
                               tGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTc                                                                                                          <  2:987727/150‑1 (MQ=255)
                                 aaGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTg                                                                                                        <  1:741236/150‑1 (MQ=255)
                                 aaGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTACCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTg                                                                                                        <  1:98755/150‑1 (MQ=255)
                                   gACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGga                                                                                                      <  1:470850/150‑1 (MQ=255)
                                    aCCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGgag                                                                                                     <  1:1014057/150‑1 (MQ=255)
                                    aCCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAACGTGGTCCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGgag                                                                                                     >  1:1396767/1‑150 (MQ=255)
                                      cGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGAt                                                                                                   >  2:816737/1‑150 (MQ=255)
                                          aaGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGc                                                                                               >  2:382497/1‑150 (MQ=255)
                                            gTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCAt                                                                                             >  1:1108615/1‑150 (MQ=255)
                                            gTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCAt                                                                                             <  2:813778/150‑1 (MQ=255)
                                            gTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCACCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCAt                                                                                             <  2:813862/150‑1 (MQ=255)
                                              cGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGc                                                                                                                                                      >  1:545536/1‑91 (MQ=255)
                                              cGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGc                                                                                                                                                      <  2:545536/91‑1 (MQ=255)
                                                    gggCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTgg                                                                                     <  1:816737/150‑1 (MQ=255)
                                                     ggCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGt                                                                                    >  1:159137/1‑150 (MQ=255)
                                                      gCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTg                                                                                   >  2:384817/1‑150 (MQ=255)
                                                       ccAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTg                                                                                                        >  1:500410/1‑128 (MQ=255)
                                                       ccAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTg                                                                                                        <  2:500410/128‑1 (MQ=255)
                                                           ccTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGa                                                                              <  2:1102153/150‑1 (MQ=255)
                                                             tGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGGTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGc                                                                            >  1:379793/1‑150 (MQ=255)
                                                                 cgcgcgCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCCGCGCCAGGCCGGCTTGCTGTCGCATGGGTAGGAGCATTGTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCa                                                                        >  1:673525/1‑150 (MQ=255)
                                                                      gcCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgc                                                                   >  1:1252077/1‑150 (MQ=255)
                                                                        cTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTTCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgctg                                                                 >  1:1056363/1‑150 (MQ=255)
                                                                         tGAGCGTTGACCAGCGACAGCGCATCGAGGACGTGGTGCCCACTTTGCGCCTGCTGCCGCTGGCCCAGAGCCAGGCCGGCTTGCGGTCTCATGGGCAGAAGCAGTTTCTGGAGATCGACCAGCTGCTGGGGCAGGAGCCGCATctgcatc                                                                >  2:95053/1‑147 (MQ=255)
                                                                         tGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgctgc                                                                <  2:797184/150‑1 (MQ=255)
                                                                               gTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGATATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCg                                                          <  2:1032377/150‑1 (MQ=255)
                                                                               gTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGc                                                                                               <  1:336174/113‑1 (MQ=255)
                                                                               gTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTGCTGGAGATCGGc                                                                                               >  2:336174/1‑113 (MQ=255)
                                                                                  aCCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGAATCCTGAAAATCGGCCAGCGCCAGGCCGGCCTTCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATg                                                       >  1:27142/1‑150 (MQ=255)
                                                                                           agcgcATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTgg                                              >  2:1243346/1‑150 (MQ=255)
                                                                                            gcgcATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGc                                             >  1:758599/1‑150 (MQ=255)
                                                                                            gcgcATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGc                                             >  2:187279/1‑150 (MQ=255)
                                                                                                      gACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGAcc                                   >  2:926315/1‑150 (MQ=255)
                                                                                                       aCGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCg                                  >  1:1267223/1‑150 (MQ=255)
                                                                                                            cTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGcc                             >  1:447294/1‑150 (MQ=255)
                                                                                                              gACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCga                           <  1:550724/150‑1 (MQ=255)
                                                                                                               accaccTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgctgc                                                                <  2:1041242/112‑1 (MQ=255)
                                                                                                               accaccTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgctgc                                                                >  1:1041242/1‑112 (MQ=255)
                                                                                                                 caccTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGAc                        >  1:1081423/1‑150 (MQ=255)
                                                                                                                     ttGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGaa                    <  2:773029/150‑1 (MQ=255)
                                                                                                                        gccggcctggCGCTGGGCCAGCGGCA‑GCAGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGAATTTa                <  2:221240/144‑1 (MQ=255)
                                                                                                                          ccTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGAATTTAc               <  2:1081423/150‑1 (MQ=255)
                                                                                                                              ctggcGCTGGGCCAGCGGCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgctg                                                                 <  2:354368/92‑1 (MQ=255)
                                                                                                                              ctggCGCTGGGCCAGCGGCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActgctg                                                                 >  1:354368/1‑96 (MQ=255)
                                                                                                                               tggcGCTGGGCCAGCGGCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActg                                                                    <  2:1382238/89‑1 (MQ=255)
                                                                                                                               tggcGCTGGGCCAGCGGCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAActg                                                                    >  1:1382238/4‑92 (MQ=255)
                                                                                                                                gcCGGTGCCCCGGGGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGAATTTACCGCCGa         <  2:159137/150‑1 (MQ=255)
                                                                                                                                     tGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGAATTTACCGCCGAGCTGt    <  2:1267223/150‑1 (MQ=255)
                                                                                                                                       gCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAAAAGGCGCGCGGTGGAACAGATACCCAAAGCGAAGGTAGCGCCGAGCTGTTc  >  2:1084612/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
CTGACCGTGTTCGAGAACCTGGAGCTGGCGCTGAAGACCGACAAGTCGGTGTGGGCCAGCCTGGCCGCGCGCCTGAGCGGTGACCAGCGCCAGCGCATCGAGGACGTGCTGACCACCTTGCGCCTGCTGCCGCTGGCCCAGCGCCAGGCCGGCTTGCTGTCGCATGGGCAGAAGCAGTTTCTGGAGATCGGCATGCTGCTGGTGCAGGAGCCGCAACTGCTGCTGCTCGATGAACCGGTGGCCGGCATGACCGATGCCGAGACCGAATTTACCGCCGAGCTGTTC  >  NC_002947/5512815‑5513099

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