Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,915,148 A→C 6.8% Q39H (CAA→CAC PP_3455 → multidrug RND transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,915,1480AC 6.8% 187.3 / 4.1 59Q39H (CAA→CACPP_3455multidrug RND transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base C (2/2);  total (30/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.61e-01

GCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTACCGCGAAGGCAGCGACGTCAA  >  NC_002947/3915001‑3915294
                                                                                                                                                   |                                                                                                                                                   
gCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGTGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCaa                                                                                                                                                     <  4:329077/148‑1 (MQ=255)
gCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCaa                                                                                                                                                     <  5:344486/148‑1 (MQ=255)
gCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAg                                                                                                                                                    >  8:90994/1‑149 (MQ=255)
          ggTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCACTTGcgcggtgg                                                                                                                                          <  5:296283/149‑5 (MQ=255)
            tcatcaTCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGaaa                                                                                                                                        <  5:125739/149‑1 (MQ=255)
                 atcGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCt                                                                                                                                    <  3:220748/148‑1 (MQ=255)
                 atcGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTg                                                                                                                                   <  3:260455/149‑1 (MQ=255)
                               aTGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCt                                                                                                                                    >  6:372499/1‑134 (MQ=255)
                               aTGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCt                                                                                                                                    <  5:372499/134‑1 (MQ=255)
                                 gCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTgg                                                                                                                   >  4:145751/1‑149 (MQ=255)
                                         cccTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATcag                                                                                                            >  5:20528/1‑148 (MQ=255)
                                              tcCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGt                                                                                                        <  6:20528/147‑1 (MQ=255)
                                                cccGTGCCTGCGCCCGTTGGCGCTCATGGCGTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTg                                                                                                                                   >  1:31242/1‑118 (MQ=255)
                                                cccGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTg                                                                                                                                   <  2:31242/118‑1 (MQ=255)
                                                cccGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTg                                                                                                       <  2:206112/146‑1 (MQ=255)
                                                 ccGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACt                                                                                                   >  2:104045/1‑149 (MQ=255)
                                                 ccGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACt                                                                                                   >  4:239948/1‑149 (MQ=255)
                                                   gTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCCGCCCCTGGCCATCAGCAGTGAACt                                                                                                   >  2:235687/1‑147 (MQ=255)
                                                   gTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCa                                                                                                 <  1:353044/149‑1 (MQ=255)
                                                      ccTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCg                                                                                              >  2:294052/1‑149 (MQ=255)
                                                                ttGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCt                                                                                    >  5:275230/1‑149 (MQ=255)
                                                                   gcgcTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGcc                                                                                           <  7:217359/139‑1 (MQ=255)
                                                                   gcgcTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGcc                                                                                           >  8:217359/1‑139 (MQ=255)
                                                                        cATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATcag                                                                                                            >  1:157914/1‑117 (MQ=255)
                                                                        cATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATcag                                                                                                            <  2:157914/117‑1 (MQ=255)
                                                                           ggCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGAAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCa                                                                         >  5:281724/1‑149 (MQ=255)
                                                                                     aGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGaa                                                                 <  6:185368/147‑1 (MQ=255)
                                                                                       ccTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAgtgc                                                             >  6:305104/1‑149 (MQ=255)
                                                                                         tGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAgtgcgt                                                           >  6:158214/1‑149 (MQ=255)
                                                                                          gTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAgtgcgtg                                                          >  6:21820/1‑149 (MQ=255)
                                                                                              cTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCACTTGcgcggtg                                                                                                                                           <  1:18388/64‑4 (MQ=255)
                                                                                              cTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCACTTGcgcggtg                                                                                                                                           >  2:18388/1‑61 (MQ=255)
                                                                                              cTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGcgcg                                                      <  4:359593/149‑1 (MQ=255)
                                                                                                ggCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAgt                                                               <  8:290372/138‑1 (MQ=255)
                                                                                                ggCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAgt                                                               >  7:290372/1‑138 (MQ=255)
                                                                                                  ccGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTc                                                  >  7:190993/1‑149 (MQ=255)
                                                                                                     gCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCa                                               <  1:148039/149‑1 (MQ=255)
                                                                                                      cTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAg                                              <  8:190993/149‑1 (MQ=255)
                                                                                                      cTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAg                                              <  5:305104/149‑1 (MQ=255)
                                                                                                         cAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAgg                                             <  3:239948/147‑1 (MQ=255)
                                                                                                         cAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCg                                           <  1:104045/149‑1 (MQ=255)
                                                                                                           gCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCgtg                                         <  3:145751/149‑1 (MQ=255)
                                                                                                           gCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCgtg                                         <  5:158214/149‑1 (MQ=255)
                                                                                                            cGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCgtgg                                        <  7:38548/149‑1 (MQ=255)
                                                                                                             gACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCgtggt                                       >  3:252740/1‑149 (MQ=255)
                                                                                                              aCAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCgtggtg                                      >  6:174756/1‑149 (MQ=255)
                                                                                                               cAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGc                                     <  1:294052/149‑1 (MQ=255)
                                                                                                                   gcCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTc                                  >  1:176031/1‑147 (MQ=255)
                                                                                                                   gcCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGAGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGAGTGCGCGGGTCGCAGGCGTGGTGCTc                                  >  4:235554/1‑147 (MQ=255)
                                                                                                                         ccAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGgt                          >  3:193058/1‑149 (MQ=255)
                                                                                                                         ccAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGgt                          <  4:252740/149‑1 (MQ=255)
                                                                                                                            gACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTa                       >  5:104093/1‑149 (MQ=255)
                                                                                                                             aCGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTAc                      <  6:104093/149‑1 (MQ=255)
                                                                                                                               gAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTACcg                    <  6:281724/149‑1 (MQ=255)
                                                                                                                                    ggCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCCCGGGTCGCAGGCGTGGTGCTCAAGCGGGAGTACCGCGAAg               >  2:311408/1‑149 (MQ=255)
                                                                                                                                       ccccACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTACCGCGAAGGCa            <  6:275230/149‑1 (MQ=255)
                                                                                                                                       ccccACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGCGTACCGCGAAGGCa            >  4:252986/1‑149 (MQ=255)
                                                                                                                                           tacgcgCACTTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTACCGCGAAGGCAGCGa        >  3:217075/3‑149 (MQ=255)
                                                                                                                                               cgcAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGgt                          >  7:273203/1‑127 (MQ=255)
                                                                                                                                               cgcAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGgt                          <  8:273203/127‑1 (MQ=255)
                                                                                                                                                gcAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTACCGCGAAGGCAGCGACGTc    >  4:169756/1‑148 (MQ=255)
                                                                                                                                                  aaGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCCAAGCGGGTGTACCGCGAAGGCAGCGACGTCaa  <  1:235687/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCATATACACGGTCATCATCGGAAGATCCGCATGCCTACTACCCTCTCCCCGTGCCTGCGCCCGTTGGCGCTCATGGCCTTGCTGAGCCTGTCCCTGGCCGGCTGCAGCGACAGCGCCGAGCCAGACGAACAGGCCCCCACCGCGCAAGTGCGCGTAGAAACCCTGCAACTGCAGCCACTGGCCATCAGCAGTGAACTCAGCGGCCGTATCCTCGCGCCACGCACCGCCGAAGTGCGTGCGCGGGTCGCAGGCGTGGTGCTCAAGCGGGTGTACCGCGAAGGCAGCGACGTCAA  >  NC_002947/3915001‑3915294

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