Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 47,142 C→G 8.0% A745P (GCG→CCG) ‡ czcA‑I ← cation efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294747,1420CG 8.0% 195.5 / 5.7 62A745P (GCG→CCG) ‡czcA‑Ication efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/3);  total (31/31)
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.83e-01

CGGGATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGAG  >  NC_002947/46996‑47286
                                                                                                                                                  |                                                                                                                                                 
cGGGATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGa                                                                                                                                                 <  5:59415/149‑1 (MQ=255)
   gATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGt                                                                                                                                              >  6:184469/1‑149 (MQ=255)
    aTCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGtt                                                                                                                                             <  3:118427/149‑1 (MQ=255)
       aGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCGCGGATGTTCa                                                                                                                                           <  3:37860/148‑1 (MQ=255)
        gCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCGCGGATGTTCAga                                                                                                                                         <  1:359097/149‑2 (MQ=255)
                   aTACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCAGCGATGTTCAGGCCGTAGCGtgc                                                                                                                              >  6:350449/1‑149 (MQ=255)
                                    gCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCaa                                                                                                            >  5:182441/1‑149 (MQ=255)
                                      gAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATg                                                                                                          <  7:219932/149‑1 (MQ=255)
                                       aaCGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGt                                                                                                         <  4:204318/149‑1 (MQ=255)
                                           gTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGAt                                                                                                     <  5:184469/149‑1 (MQ=255)
                                              tctGGTAGGCGTACTACCAAGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATgg                                                                                                   >  3:164563/1‑148 (MQ=255)
                                                 ggTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTc                                                                                                                    <  7:297921/128‑1 (MQ=255)
                                                 ggTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTc                                                                                                                    >  8:297921/1‑128 (MQ=255)
                                                  gTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCa                                                                                                <  2:7655/147‑1 (MQ=255)
                                                            tacCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGc                                                                                    >  2:136934/1‑149 (MQ=255)
                                                             acCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGAAGGcc                                                                                   >  6:169016/1‑149 (MQ=255)
                                                              ccATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCg                                                                                  >  7:205825/1‑149 (MQ=255)
                                                              ccATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCg                                                                                  >  7:53835/1‑149 (MQ=255)
                                                              ccATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCg                                                                                  >  6:312107/1‑149 (MQ=255)
                                                              ccATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGCCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCTCGGTCAATGATGATGGTCAGCACCGGCAGGCCg                                                                                  >  4:322218/1‑149 (MQ=255)
                                                               cATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCg                                                                                  >  8:148091/1‑148 (MQ=255)
                                                                aTGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGa                                                                                >  3:127546/1‑149 (MQ=255)
                                                                     gAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAgg                                                                                                                                         <  1:131156/88‑1 (MQ=255)
                                                                     gAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAgg                                                                                                                                         >  2:131156/1‑88 (MQ=255)
                                                                      aaCCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTg                                                                          >  7:33044/1‑149 (MQ=255)
                                                                         cGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTc                                                                       <  4:127546/149‑1 (MQ=255)
                                                                               tCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACttt                                                                 <  4:118234/149‑1 (MQ=255)
                                                                                 gCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTc                                                                <  1:136934/148‑1 (MQ=255)
                                                                                 gCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTc                                                                <  5:350449/148‑1 (MQ=255)
                                                                                  cccTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACtt                                                                  >  6:17703/1‑145 (MQ=255)
                                                                                  cccTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACtt                                                                  <  5:17703/145‑1 (MQ=255)
                                                                                   ccTCATACAGCGTGCCGGCCTGACGGCCCCCCACGGCGATAGCGATCGAGTTCTGAACATCCGCCATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTGAATGTTGATGGTCAGCACCGGCAGGCCGGATGACTGCTCAACTTTCACt                                                             >  8:114798/1‑149 (MQ=255)
                                                                                    cTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACtt                                                            <  8:205825/149‑1 (MQ=255)
                                                                                       aTACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTc                                                           <  7:148091/147‑1 (MQ=255)
                                                                                        tACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGa                                                        <  5:6668/149‑1 (MQ=255)
                                                                                        tACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTAAGGCCGTAGCGTGCTGCT‑TTTTCGCCGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTAAACTTTAACTTCCGa                                                        >  8:155972/1‑149 (MQ=255)
                                                                                          cAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGAt                                                       <  6:182441/148‑1 (MQ=255)
                                                                                                gCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCg                                                >  4:67638/1‑149 (MQ=255)
                                                                                                 ccGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCTCGGGCAATGTTGATGGTCAGCTCCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCg                                                <  5:140306/148‑1 (MQ=255)
                                                                                                  cGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTc                                                                                                 >  8:212632/1‑98 (MQ=255)
                                                                                                  cGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTc                                                                                                 <  7:212632/98‑1 (MQ=255)
                                                                                                    gCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTc                                                                >  3:100139/1‑129 (MQ=255)
                                                                                                    gCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTc                                                                <  4:100139/129‑1 (MQ=255)
                                                                                                    gCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGAc                                            >  4:254192/1‑149 (MQ=255)
                                                                                                          cGGCCACCCACGTCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTATCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCTGAAGGCCGGATGTCTGCTCAACTTTCAGTTCCGATGATCCCGGGACCGCttt                                      <  4:164563/149‑1 (MQ=255)
                                                                                                          cGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCttt                                      <  5:312107/149‑1 (MQ=255)
                                                                                                                 ccacGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCt                               >  7:190458/1‑149 (MQ=255)
                                                                                                                 ccacGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCt                               >  4:176280/1‑149 (MQ=255)
                                                                                                                 ccacGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGAACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGATCCCGGGACCGCTTTGAGCGCt                               >  5:277573/1‑149 (MQ=255)
                                                                                                                  cacGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTg                              <  8:190458/149‑1 (MQ=255)
                                                                                                                    cGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTGAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGcc                            <  8:33044/149‑1 (MQ=255)
                                                                                                                           aGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCt                     >  5:78739/1‑149 (MQ=255)
                                                                                                                               aTCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCttgtt                 >  6:35678/1‑149 (MQ=255)
                                                                                                                                       cTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGgtg          <  3:67638/148‑1 (MQ=255)
                                                                                                                                        tGAACATCGCGGATGTTCAGGCCGTAGCGTNCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGtgtt        >  1:30185/1‑149 (MQ=255)
                                                                                                                                        tGAACATCGCGGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCttgttg                <  4:230373/141‑1 (MQ=255)
                                                                                                                                        tGAACATCGCGGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCttgttg                >  3:230373/1‑141 (MQ=255)
                                                                                                                                          aaCATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGtgttgt      >  2:36822/1‑149 (MQ=255)
                                                                                                                                          aaCATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGAGGtgtt        <  5:169016/147‑1 (MQ=255)
                                                                                                                                             aTCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGa   <  5:35678/149‑1 (MQ=255)
                                                                                                                                             aTCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGa   <  3:176280/149‑1 (MQ=255)
                                                                                                                                               ccGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGAg  >  4:132431/1‑148 (MQ=255)
                                                                                                                                                cgcgATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGa                                    >  1:12441/1‑113 (MQ=255)
                                                                                                                                                cgcgATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGa                                    <  2:12441/113‑1 (MQ=255)
                                                                                                                                                cgcgATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTgg               <  3:4526/134‑1 (MQ=255)
                                                                                                                                                cgcgATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTgg               >  4:4526/1‑134 (MQ=255)
                                                                                                                                                 gcgATGTTCAGGCCGTAGCGTGCTGCTATTTTCGCGGTAAATGTTGATGTTCAGAACCGGCAGGCCGGATGTCTGATCAACTTTAAATTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGAt                       >  1:205667/1‑126 (MQ=255)
                                                                                                                                                 gcgATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGAt                       <  2:205667/125‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGGGATCAGCAAGCTGGACATACCCGCTACGTCGGTGCGAACGGTCTCTGGTAGGCGTACTACCATGTCGAACCGACGGTCGCCCTCATACAGCGTGCCGGCCTGACGGCCACCCACGGCGATAGCGATCGAGTTCTGAACATCCGCGATGTTCAGGCCGTAGCGTGCTGCT‑TTTTCGCGGTCAATGTTGATGGTCAGCACCGGCAGGCCGGATGTCTGCTCAACTTTCACTTCCGATGAGCCCGGGACCGCTTTGAGCGCTGCCGCGATCTTGTTGGCGGTGTTGTTGAG  >  NC_002947/46996‑47286

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