Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,137,485 T→A 5.9% intergenic (+29/‑61) PP_5733 → / → czcC metal‑binding protein/CzcC family metal RND transporter outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,137,4850TA 5.9% 210.7 / 2.5 68intergenic (+29/‑61)PP_5733/czcCmetal‑binding protein/CzcC family metal RND transporter outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base A (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GTCAGCATGGTCGAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCG  >  NC_002947/6137343‑6137632
                                                                                                                                              |                                                                                                                                                   
gTCAGCATGGTCGAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcg                                                                                                                                               >  3:88860/1‑149 (MQ=255)
            gAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGAGACCTACTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGa                                                                                                                                   >  5:42907/1‑149 (MQ=255)
             aGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                  <  8:100242/149‑1 (MQ=255)
                ccGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGTTCAGAGCCTTTTATCCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATCGTTCAGCCGCCGAGGCAGAACAATggt                                                                                                                               <  2:112457/149‑3 (MQ=255)
                ccGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                  >  6:53878/1‑146 (MQ=255)
                ccGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGAt                                                                                                                                  <  5:53878/146‑1 (MQ=255)
                  gCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGacac                                                                                                                             >  6:101636/1‑149 (MQ=255)
                    aCTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcggc                                                                                                                                             >  2:106525/1‑131 (MQ=255)
                    aCTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTCGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcggc                                                                                                                                             <  1:106525/131‑1 (MQ=255)
                        aCCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCt                                                                                                                       >  8:149169/1‑149 (MQ=255)
                          cTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCtgt                                                                                                                     <  4:88860/149‑1 (MQ=255)
                               ccACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCt                                                                                                                <  6:180357/149‑1 (MQ=255)
                                cACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTg                                                                                                               <  5:263173/149‑1 (MQ=255)
                                    acacGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGttt                                                                                                           >  5:221100/1‑149 (MQ=255)
                                     cacGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGtttt                                                                                                          >  8:29815/1‑149 (MQ=255)
                                        gtacAGTCGATACCTACGGTCAGAGCCTTTTAACAAGGGCTCCCCCCGATCACTGGAGCCCTCATCGCCCCTGATTCCTTCGTAAAACTCAAACCATTGTTCTGCCTAGTCGGCTGAACGATGACACGCGACTGTGCGCTCGTTTTTTc                                                                                                       >  6:197753/4‑149 (MQ=255)
                                           cAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGa                                                                                                     <  6:55378/148‑1 (MQ=255)
                                           cAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGa                                                                                                     <  7:149169/148‑1 (MQ=255)
                                           cAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAg                                                                                                    <  5:255003/149‑1 (MQ=255)
                                             gTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAgg                                                                                                   >  7:27830/1‑148 (MQ=255)
                                                     cTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCAtctc                                                                                          >  2:8014/1‑149 (MQ=255)
                                                       aCGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCAAGGAATCATCTCCa                                                                                        >  1:35100/1‑149 (MQ=255)
                                                          gTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGa                                                                                     >  5:217827/1‑149 (MQ=255)
                                                             agagCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTc                                                                                  >  4:88882/1‑149 (MQ=255)
                                                              gagCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCATGATTCATTCGTAAAAATCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTcc                                                                                 >  7:288409/1‑149 (MQ=255)
                                                               agCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTccc                                                                                >  3:37210/1‑149 (MQ=255)
                                                                gCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTccc                                                                                <  4:37210/148‑1 (MQ=255)
                                                                   ttttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcggcgg                                                                                                                                           <  7:231737/86‑1 (MQ=255)
                                                                   ttttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcggcgg                                                                                                                                           <  7:236001/86‑1 (MQ=255)
                                                                   ttttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcggcgg                                                                                                                                           >  8:236001/1‑86 (MQ=255)
                                                                   ttttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTcggcgg                                                                                                                                           >  8:231737/1‑86 (MQ=255)
                                                                   ttttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTAt                                                                            >  4:34238/1‑149 (MQ=255)
                                                                    tttAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATg                                                                           >  3:300837/1‑149 (MQ=255)
                                                                          cAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCa                                                                     >  7:217772/1‑149 (MQ=255)
                                                                            gggCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCAt                                                                    >  6:98831/1‑148 (MQ=255)
                                                                                  cccccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGc                                                             >  5:53656/1‑149 (MQ=255)
                                                                                  cccccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCAAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGc                                                             <  2:35100/149‑1 (MQ=255)
                                                                                    cccGATCACTGGCGCCCTCCTCGCGCCTGATTACTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCgg                                                           <  5:197753/149‑1 (MQ=255)
                                                                                     ccGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGa                                                          >  1:94312/1‑149 (MQ=255)
                                                                                           aCTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGAtgctt                                                     <  8:27830/148‑1 (MQ=255)
                                                                                                  cccTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTgccgc                                             <  8:217772/149‑1 (MQ=255)
                                                                                                   cctcctCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGATCCCATTCAGGGCGGATGCTtgctgc                                                >  1:19634/1‑145 (MQ=255)
                                                                                                   cctactCGCTCCTGATTCATTCGTAAAACTCACACCATTGTTCTGCCTCGTCGGCTGAAAGATGACACGCTCCTGTTCGCTTGTTTTTTCGAGGAATAATCTCAATGAATCCCCTATTATACCATTCAGGGCTTAtgctttatgc                                                <  2:19634/145‑6 (MQ=255)
                                                                                                    ctcctcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCt                                           >  1:55227/1‑149 (MQ=255)
                                                                                                       ctcGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTg                                         <  2:55227/148‑1 (MQ=255)
                                                                                                         cgcgCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGt                                       >  7:206281/1‑148 (MQ=255)
                                                                                                         cgcgCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTa                                      <  5:101636/149‑1 (MQ=255)
                                                                                                          gcgcCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTa                                      <  3:34238/148‑1 (MQ=255)
                                                                                                            gcCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGATGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCGTTTGTa                                      <  6:42907/146‑1 (MQ=255)
                                                                                                                gATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCctgc                               >  8:236512/1‑149 (MQ=255)
                                                                                                                    ccTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCtgctgc                            <  6:217827/148‑1 (MQ=255)
                                                                                                                     cTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGAtgc                                                       <  7:233994/120‑1 (MQ=255)
                                                                                                                     cTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGAtgc                                                       >  8:233994/1‑120 (MQ=255)
                                                                                                                     cTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCg                          >  6:237510/1‑149 (MQ=255)
                                                                                                                        cGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGgc                        <  1:8014/148‑1 (MQ=255)
                                                                                                                         gTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGgcgc                      <  5:98831/149‑1 (MQ=255)
                                                                                                                          tAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTCGTACTCCTGCTGCCGgcgc                      >  6:274176/1‑148 (MQ=255)
                                                                                                                            aaaCTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTAt                   <  8:127839/149‑1 (MQ=255)
                                                                                                                             aaCTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATc                  <  8:288409/149‑1 (MQ=255)
                                                                                                                                    acCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCtt                                          >  7:267855/1‑118 (MQ=255)
                                                                                                                                    acCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCtt                                          <  8:267855/118‑1 (MQ=255)
                                                                                                                                       atcgTTCAGCCGCCGAGGCAGAACAATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGATCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCTGCGCTATCGGCGCAGGc         >  3:233800/4‑148 (MQ=255)
                                                                                                                                       aTTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTc                                   <  8:276130/122‑1 (MQ=255)
                                                                                                                                       aTTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGc         >  5:184647/1‑148 (MQ=255)
                                                                                                                                       aTTGTGCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTc                                   >  7:276130/1‑122 (MQ=255)
                                                                                                                                        tcgTTCAGCCGCCGAGGCAGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGa                                                          >  2:299442/3‑98 (MQ=255)
                                                                                                                                        tcgTTCAGCCGCCGAGGCAGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGa                                                          <  1:299442/96‑1 (MQ=255)
                                                                                                                                           ttCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACt    <  4:300837/149‑1 (MQ=255)
                                                                                                                                              tGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCg  >  7:16770/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GTCAGCATGGTCGAGCCCGCACTCACCTTCGCCACTACACGTCCAGTCGATACCTACGGTCAGAGCCTTTTAACCAGGGCTCCCCCCGATCACTGGCGCCCTCCTCGCGCCTGATTCCTTCGTAAAACTCACACCATTGTTCTGCCTCGGCGGCTGAACGATGACACGCGCCTGTGCGCTGGTTTTTTCGAGGAATCATCTCCATGACTCCCCTATGCTCCCATTCAGGGCGGATGCTTGCTGCCGCCTTGGTACTCCTGCTGCCGGCGCTATCGGCGCAGGCGAACTCG  >  NC_002947/6137343‑6137632

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: