Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,635,709 C→T G16D (GGT→GAT)  uvrY ← BarA/UvrY two‑componenent system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,635,7090CT84.9% 168.0 / 20.7 73G16D (GGT→GAT) uvrYBarA/UvrY two‑componenent system response regulator
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base T (31/31);  total (37/36)
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.91e-01

CCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAGAAACGACTGACGT  >  NC_002947/4635563‑4635852
                                                                                                                                                  |                                                                                                                                               
cttCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCCGCACAACTTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCAGCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCg                                                                                                                                               <  8:279225/147‑1 (MQ=255)
   ccAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACTTCCGTGTAATATCTGtt                                                                                                                                            >  1:238588/1‑149 (MQ=255)
    cAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTc                                                                                                                                           <  4:30796/149‑1 (MQ=255)
    cAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTc                                                                                                                                           >  2:110772/1‑149 (MQ=255)
    cAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTc                                                                                                                                           >  4:132968/1‑149 (MQ=255)
        ccgccgATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTAc                                                                                                                                       <  8:130064/149‑1 (MQ=255)
        ccgccgATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTAc                                                                                                                                       <  8:240540/149‑1 (MQ=255)
          gccgATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCa                                                                                                                                     <  4:106998/149‑1 (MQ=255)
               tGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCg                                                                                                                                <  6:199043/149‑1 (MQ=255)
                  cGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGa                                                                                                                             >  1:212246/1‑149 (MQ=255)
                   gggCATTTTCACATCCATCTGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGat                                                                                                                            >  4:99869/1‑149 (MQ=255)
                       aTTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGatcatc                                                                                                                        <  6:44552/149‑1 (MQ=255)
                           tCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGAc                                                                                                                     >  5:335955/1‑148 (MQ=255)
                               atccatcAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACt                                                                                                                 >  3:369739/1‑148 (MQ=255)
                                    tcAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGAcc                                                                                                           >  5:263844/1‑149 (MQ=255)
                                       gCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGAcc                                                                                                                    <  8:247843/137‑1 (MQ=255)
                                       gCACAACGTCCGGTTTCAGCTCCCAGGCCAGTTTCAGTTCCGACTAGCCGGAATATCCCTCACCCCCCACCTGCAGGCCGTCGATGTGGGCCAACATCCGGGTAATTTAGGTTAGTTCCAGCTCGTGATGATCGAcc                                                                                                                    >  7:247843/1‑137 (MQ=255)
                                        cacaACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTaa                                                                                                       >  1:275089/1‑149 (MQ=255)
                                        cacaACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTaa                                                                                                       <  6:263844/149‑1 (MQ=255)
                                        cacaACGTCCGGCTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTaa                                                                                                       <  5:342953/149‑1 (MQ=255)
                                           aaCGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGtt                                                                                                                                            >  7:192469/1‑109 (MQ=255)
                                           aaCGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGtt                                                                                                                                            <  8:192469/109‑1 (MQ=255)
                                             cGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAg                                                                                                  <  2:212246/149‑1 (MQ=255)
                                                 cGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGcaca                                                                                              <  2:207885/149‑1 (MQ=255)
                                                  ggTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGcacac                                                                                             <  4:369739/149‑1 (MQ=255)
                                                   gTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGcacaca                                                                                            >  4:72451/1‑149 (MQ=255)
                                                      tCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAGCATCCGTGTAATATCGGTTCGTACCAGATCGTGATGATCGCCCACTAAGACCCTAATCAAGCACACACCt                                                                                         >  7:206803/1‑149 (MQ=255)
                                                       cAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTg                                                                                                                              <  5:357178/111‑1 (MQ=255)
                                                       cAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTg                                                                                                                              >  6:357178/1‑111 (MQ=255)
                                                          cTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGca                                                                                                >  2:390839/1‑138 (MQ=255)
                                                          cTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGca                                                                                                <  1:390839/138‑1 (MQ=255)
                                                           tCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGCCCACTAAGACCCTAATCAAGCACACACCTCGt                                                                                      >  2:310496/1‑147 (MQ=255)
                                                               gggCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGc                                                                                                 <  5:334241/132‑1 (MQ=255)
                                                               gggCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGc                                                                                                 >  6:334241/1‑132 (MQ=255)
                                                                   cAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACTCTAATCAAGCACACACCTCGTCAGCAGGGc                                                                             >  3:52788/1‑148 (MQ=255)
                                                                          aGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGcaca                                                                                              >  3:59184/1‑124 (MQ=255)
                                                                          aGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGcaca                                                                                              <  4:59184/124‑1 (MQ=255)
                                                                                 aCTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCa                                                              <  6:226546/149‑1 (MQ=255)
                                                                                   tCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATc                                                                                                     <  4:337361/108‑1 (MQ=255)
                                                                                   tCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATc                                                                                                     >  3:337361/1‑108 (MQ=255)
                                                                                                ccTCACCCAGCACCTTCAGGCCGTCGATGGCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAAtt                                                 <  2:238588/147‑1 (MQ=255)
                                                                                                ccTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAg                                                                                                  >  4:47087/1‑98 (MQ=255)
                                                                                                ccTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAg                                                                                                  <  3:47087/98‑1 (MQ=255)
                                                                                                ccTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACaaaa                                                   <  4:146438/145‑1 (MQ=255)
                                                                                                ccTCACCCACCACCTGCAGGCCGTCGATATCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAAtttt                                               >  2:135736/1‑149 (MQ=255)
                                                                                                ccTCACAAACCACCTCCAGGCCGTCGATTTCGGCCAACATCCCGCTAATATCGGTTCGTACCAGAACGTGATCATCGACCACTAAGACCCGAATCAAGCACACACCTCGACAGCAGGGCGAATGGATCCCGCCACCTTAACAAAAtttt                                               >  2:318685/1‑149 (MQ=255)
                                                                                                 cTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGAt                                                                          <  6:38271/121‑1 (MQ=255)
                                                                                                 cTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGAt                                                                          >  5:38271/1‑121 (MQ=255)
                                                                                                 cTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAAttttt                                              <  1:110772/149‑1 (MQ=255)
                                                                                                  tCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAAttttt                                              >  1:39595/1‑148 (MQ=255)
                                                                                                    acccaccACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCAc                                           <  6:335955/149‑1 (MQ=255)
                                                                                                      ccaccaCCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCa                                            <  7:302570/146‑1 (MQ=255)
                                                                                                      ccaccaCCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCa                                            >  8:302570/1‑146 (MQ=255)
                                                                                                           accTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCATGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGc                                    >  7:139019/1‑149 (MQ=255)
                                                                                                             cTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTgg                                      <  8:140462/145‑1 (MQ=255)
                                                                                                             cTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTgg                                      >  7:140462/1‑145 (MQ=255)
                                                                                                             cTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGcc                                   >  6:154253/1‑148 (MQ=255)
                                                                                                               gCAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTaa                                >  5:240191/1‑149 (MQ=255)
                                                                                                                cAGGCCGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTa                                 <  2:39595/147‑1 (MQ=255)
                                                                                                                     cGTCGATTTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCAtt                          <  1:135736/149‑1 (MQ=255)
                                                                                                                     cGTCGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCagt                          >  6:52063/1‑147 (MQ=255)
                                                                                                                       tCGATGTCTGCCAACATCCGTTTAAAATGGGTTCGTACCAGATCGTCATCATCGACCACTAAGACCCTTATCAAGCACACTCCTCGTCAGCAGTGCGATTGGATCCCGCCACCTTAACAAAATTTATCACGGTGGGCCTAACGCATTGc                        <  1:318685/149‑1 (MQ=255)
                                                                                                                        cGATGTCGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGct                       >  4:350120/1‑149 (MQ=255)
                                                                                                                          aTGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGcttc                     >  3:126946/1‑149 (MQ=255)
                                                                                                                              cGGCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGcttcctt                  >  7:68825/1‑148 (MQ=255)
                                                                                                                               ggCCAACATCCGGGTAATACCGGTTCGTACCAGATCGTGATCAGCGACCACTAAGACCCTAATCAAGCACACACCTCGTCTGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCa                <  3:99869/149‑1 (MQ=255)
                                                                                                                                gCCAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAg               <  8:68825/149‑1 (MQ=255)
                                                                                                                                 ccAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAGa              <  3:72451/149‑1 (MQ=255)
                                                                                                                                  cAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACgg                                         >  2:394936/1‑121 (MQ=255)
                                                                                                                                  cAACATCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACgg                                         <  1:394936/121‑1 (MQ=255)
                                                                                                                                      aTCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAGAAACGa         >  3:118953/1‑149 (MQ=255)
                                                                                                                                       tCCGTGTAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAGAAACGAc        <  2:275089/149‑1 (MQ=255)
                                                                                                                                             tAATATCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAGAAACGACTGACGt  >  1:372774/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAGCACACACCTCGTCAGCAGGGCGATTGGATCCCGCCACCTTAACAAAATTTTTCACGGTGGGCCTAACGCATTGCTTCCTTCAGAAACGACTGACGT  >  NC_002947/4635563‑4635852

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: ‑