Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,698,163 C→G 6.1% G103A (GGC→GCC)  ybhN ← phospholipid modification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,698,1630CG 6.1% 254.1 / 3.8 77G103A (GGC→GCC) ybhNphospholipid modification enzyme
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base G (2/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCGACCAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAACGCGATCGCC  >  NC_002947/3698017‑3698309
                                                                                                                                                  |                                                                                                                                                  
tctcGACCAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGcc                                                                                                                                                   >  2:397856/1‑148 (MQ=255)
    gACCAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTAcc                                                                                                                                               >  3:108744/1‑148 (MQ=255)
    gACCAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCACATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCAAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTaac                                                                                                                                               >  8:158192/1‑146 (MQ=255)
     aCCAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCa                                                                                                                                             <  1:339870/149‑1 (MQ=255)
      ccAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAg                                                                                                                                            <  1:393404/149‑1 (MQ=255)
       cAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAgg                                                                                                                                           >  1:89588/1‑149 (MQ=255)
         ggCCCAGGGCAATTACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTa                                                                                                                                         <  5:156487/149‑1 (MQ=255)
            ccAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGt                                                                                                                                       >  3:100113/1‑148 (MQ=255)
            ccAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTg                                                                                                                                      <  2:89588/149‑1 (MQ=255)
               gggCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGaaa                                                                                                                                   >  2:269701/1‑149 (MQ=255)
                 gCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAAtt                                                                                                                                 >  5:19718/1‑149 (MQ=255)
                   aaaTACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTCCCCAGGTAGTGAAAtt                                                                                                                                 >  8:285119/1‑147 (MQ=255)
                    aaTACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAg                                                                                                                              <  7:122216/149‑1 (MQ=255)
                     aTACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAgg                                                                                                                             <  4:100113/149‑1 (MQ=255)
                             ccGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGc                                                                                                                     <  1:259666/149‑1 (MQ=255)
                                 ccAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAg                                                                                                                 <  8:248122/149‑1 (MQ=255)
                                     cagcaTATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAgg                                                                                                                                           >  2:332188/1‑119 (MQ=255)
                                     cagcaTATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAgg                                                                                                                                           <  1:332188/119‑1 (MQ=255)
                                        caTATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAgg                                                                                                                             >  5:246764/1‑130 (MQ=255)
                                        caTATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAgg                                                                                                                             <  6:246764/130‑1 (MQ=255)
                                        caTATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACg                                                                                                           >  2:205119/1‑148 (MQ=255)
                                          tataGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAg                                                                                                                              <  2:356562/127‑1 (MQ=255)
                                          tataGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAg                                                                                                                              >  1:356562/1‑127 (MQ=255)
                                           ataGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAg                                                                                                                              >  2:96791/1‑126 (MQ=255)
                                           ataGCCCATCCAGTTGGTCAACAGACCCAAGGTGAGGATTCGTGTAATGGTGGATGAATCCAGGCACAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTACTGAAATTGAGGTTAAAGGCGTAGCACACGAAcg                                                                                                       >  5:395486/1‑149 (MQ=255)
                                               cccATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTa                                                                                                                                         <  5:4018/111‑1 (MQ=255)
                                               cccATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTa                                                                                                                                         >  6:4018/1‑111 (MQ=255)
                                               cccATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGAc                                                                                                   >  8:331116/1‑149 (MQ=255)
                                                ccATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACg                                                                                                  <  1:397856/149‑1 (MQ=255)
                                                   tCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAg                                                                                                                                            >  6:337257/1‑104 (MQ=255)
                                                   tCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAg                                                                                                                                            <  5:337257/104‑1 (MQ=255)
                                                    ccAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGc                                                                                               >  1:275649/1‑148 (MQ=255)
                                                       gTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAAt                                                                                            >  8:392116/1‑148 (MQ=255)
                                                              aaCAGCCCCTAGGTGAGGATTCGTGTATTGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCATTACCTGGc                                                                                     <  2:324201/148‑1 (MQ=255)
                                                                cAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCggg                                                                                  >  1:42086/1‑149 (MQ=255)
                                                                 aGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGc                                                                                 <  2:275649/149‑1 (MQ=255)
                                                                         ggTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTa                                                                                                                  <  3:373406/108‑1 (MQ=255)
                                                                         ggTGAGGATTCGTGTAATGGTCGCGGTATCCAGCCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCTCCTACCCAGGTACTCAAATTGAGGTTAAACGAGTa                                                                                                                  >  4:373406/1‑108 (MQ=255)
                                                                           tGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACg                                                                                                           <  7:295441/113‑1 (MQ=255)
                                                                           tGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACg                                                                                                           >  8:295441/1‑113 (MQ=255)
                                                                               ctttCGGGTAATTGTCCCGGTTACCACGCCCCCCCGCCCATAAATGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTTTTTAAATTTTGGTTAAAGGCTTAGCACACGAACGCGACGGGCAAAACCTGGCGGGCAGGTAGCGTATg                                                                     <  7:61685/145‑1 (MQ=255)
                                                                                 ttCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCg                                                                 <  1:205119/149‑1 (MQ=255)
                                                                                   cGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGGGTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTCGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTATCGTATGGCCgg                                                                <  6:395486/148‑1 (MQ=255)
                                                                                    gtgtAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGgt                                                               <  7:331116/148‑1 (MQ=255)
                                                                                    gtgtAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGgt                                                               <  7:392116/148‑1 (MQ=255)
                                                                                     tgtAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAAtt                                                                                                                                 <  6:94905/81‑1 (MQ=255)
                                                                                     tgtAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAAtt                                                                                                                                 >  5:94905/1‑81 (MQ=255)
                                                                                        aaTGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGca                                                                                                               <  3:211940/96‑1 (MQ=255)
                                                                                        aaTGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGca                                                                                                               >  4:211940/1‑96 (MQ=255)
                                                                                        aaTGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAg                                                          >  2:143447/1‑149 (MQ=255)
                                                                                                cgGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACGCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCTGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGcc                                                  <  6:19718/149‑1 (MQ=255)
                                                                                                  gTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCaa                                                >  8:222809/1‑149 (MQ=255)
                                                                                                   tATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAAc                                               <  5:247008/149‑1 (MQ=255)
                                                                                                     tCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAg                                             >  4:407889/1‑149 (MQ=255)
                                                                                                        aGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTc                                          >  1:259413/1‑149 (MQ=255)
                                                                                                            ccAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGt                                        >  3:44547/1‑147 (MQ=255)
                                                                                                                ccgccCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCATGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTgg                                  <  7:158192/149‑1 (MQ=255)
                                                                                                                ccgccCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTgg                                  >  3:188551/1‑149 (MQ=255)
                                                                                                                  gccCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTc                                          <  8:285435/139‑1 (MQ=255)
                                                                                                                  gccCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTc                                          >  7:285435/1‑139 (MQ=255)
                                                                                                                          aGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGaa                        <  4:108744/149‑1 (MQ=255)
                                                                                                                            gCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACa                                              <  7:217499/125‑1 (MQ=255)
                                                                                                                            gCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACa                                              >  8:217499/1‑125 (MQ=255)
                                                                                                                                   cgcAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGca               <  4:44547/149‑1 (MQ=255)
                                                                                                                                    gcAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATgg                                                                    >  6:274920/1‑95 (MQ=255)
                                                                                                                                    gcAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATgg                                                                    <  5:274920/95‑1 (MQ=255)
                                                                                                                                    gcAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGca               >  7:216334/1‑148 (MQ=255)
                                                                                                                                        ggCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAAcg          <  4:188551/149‑1 (MQ=255)
                                                                                                                                         gCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAAcgc         >  5:386714/1‑149 (MQ=255)
                                                                                                                                            aCCCCGGCTACCCAGGTAGTGATATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAAc                                     <  1:334052/118‑1 (MQ=255)
                                                                                                                                            aCCCCGGCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAAc                                     >  2:338455/1‑118 (MQ=255)
                                                                                                                                            aCCCCGGCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAAc                                     >  2:334052/1‑118 (MQ=255)
                                                                                                                                            aCCCCGGCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAAc                                     <  1:338455/118‑1 (MQ=255)
                                                                                                                                            aCCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAACGCGAt      >  6:1936/1‑149 (MQ=255)
                                                                                                                                             ccccGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTCAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTCGCTGCACAACGCGATc     >  2:246205/1‑149 (MQ=255)
                                                                                                                                                cGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAACGCGATCGcc  <  3:407889/149‑1 (MQ=255)
                                                                                                                                                cGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAACGCGATCGcc  <  5:1936/149‑1 (MQ=255)
                                                                                                                                                 gCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAACGCGATCGcc  >  6:166112/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCTCGACCAGGCCCAGGGCAAATACGGTGCCGGCCAGCAGCATATAGCCCATCCAGTTGGTCAACAGCCCCAAGGTGAGGATTCGTGTAATGGTCGCGGTATCCAGGCCCAGCCGCCCATAAAGGCGGTAGCGCAAGGCAACCCCGCCTACCCAGGTAGTGAAATTGAGGTTAAAGGCGTAGCACACGAACGCGACGGGCAATACCTGGCGGGCAGGTAGCGTATGGCCGGTGTAAGCCCGGGCCAACAGGTCGTAACTGGCAAACGTGAAGTAGCTGCACAACGCGATCGCC  >  NC_002947/3698017‑3698309

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: