Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,177 G→C 6.6% W176C (TGG→TGC)  ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,1770GC 6.6% 246.6 / 5.9 76W176C (TGG→TGC) ybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (40/31);  new base C (2/3);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCGGCTGGCCT  >  NC_002947/4479030‑4479325
                                                                                                                                                   |                                                                                                                                                    
cGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTgg                                                                                                                                                      <  1:324146/148‑1 (MQ=255)
 ggATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGc                                                                                                                                                     >  5:106344/1‑148 (MQ=255)
 ggATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGc                                                                                                                                                     >  2:100870/1‑148 (MQ=255)
  gATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTTAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCaa                                                                                                                                                   <  1:100870/149‑1 (MQ=255)
       aGCAATCTGTATGGAGGGATCCCTGGAGGCCAAGATGATTACGTTAGGGTGCCCAAAGTCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGCCAAGCGGt                                                                                                                                              <  6:145997/149‑1 (MQ=255)
        gCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGCCAAGCGGtg                                                                                                                                             <  1:124785/149‑2 (MQ=255)
         ccATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAg                                                                                                                                                  <  4:132564/143‑1 (MQ=255)
         ccATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAg                                                                                                                                                  >  3:132564/1‑143 (MQ=255)
                  tGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGcc                                                                                                                                    >  7:189727/1‑148 (MQ=255)
                        gATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTg                                                                                                                             >  3:224785/1‑149 (MQ=255)
                            ccTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCg                                                                                                                          >  1:6435/1‑148 (MQ=255)
                            ccTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCg                                                                                                                          <  4:224785/148‑1 (MQ=255)
                              tGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGCCAAGCGGtg                                                                                                                                             >  8:98759/1‑126 (MQ=255)
                              tGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGCCAAGCGGtg                                                                                                                                             <  7:98759/127‑2 (MQ=255)
                              tGGAGGCCAAGCTGATTACGATAGGGTGCCCAAAGGCAATGTCGGCCCATTCAAAGACCCGACGACCCTCGCGGATGACAAACTGCTGTTCCTATCCGACCTCCTTCACACCGCTTGGCAAGCGGTCTTCAATGCACACAATGGCGAAg                                                                                                                       >  4:203333/1‑149 (MQ=255)
                                    ccAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCAt                                                                                                                 >  7:92593/1‑149 (MQ=255)
                                      aaGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAg                                                                                                                       >  8:283704/1‑141 (MQ=255)
                                      aaGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAg                                                                                                                       <  7:283704/141‑1 (MQ=255)
                                      aaGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCg                                                                                                               >  5:219423/1‑149 (MQ=255)
                                      aaGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCg                                                                                                               >  2:82652/1‑149 (MQ=255)
                                        gCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCgg                                                                                                              >  6:104368/1‑148 (MQ=255)
                                                cGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGc                                                                                                                     >  2:157815/1‑133 (MQ=255)
                                                cGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGc                                                                                                                     <  1:157815/133‑1 (MQ=255)
                                                       gTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGgcgc                                                                                              >  1:65082/1‑149 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGaa                                                                                                                        >  3:108897/1‑119 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGaa                                                                                                                        >  3:113578/1‑119 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGaa                                                                                                                        <  4:108897/119‑1 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGaa                                                                                                                        <  4:113578/119‑1 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgggg                                                                                          >  8:109464/1‑149 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgggg                                                                                          >  6:320226/1‑149 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgggg                                                                                          <  4:52242/149‑1 (MQ=255)
                                                           ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgggg                                                                                          >  6:148289/1‑149 (MQ=255)
                                                                 gCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCggg                                                                                           <  1:309670/142‑1 (MQ=255)
                                                                 gCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCggg                                                                                           >  2:309670/1‑142 (MQ=255)
                                                                 gCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTg                                                                                    <  5:104368/149‑1 (MQ=255)
                                                                  cAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTgg                                                                                   <  1:82652/149‑1 (MQ=255)
                                                                   aaTGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTggg                                                                                  >  2:135176/1‑149 (MQ=255)
                                                                       tCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTg                                                                              >  1:43335/1‑149 (MQ=255)
                                                                         gggCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGtt                                                                            <  2:6435/149‑1 (MQ=255)
                                                                              aTTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTg                                                                       <  5:320226/149‑1 (MQ=255)
                                                                              aTTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTg                                                                       >  4:168443/1‑149 (MQ=255)
                                                                                  aaaGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTgccgc                                                                   >  5:182447/1‑149 (MQ=255)
                                                                                      tCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgg                                                                                            >  7:53381/1‑120 (MQ=255)
                                                                                      tCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgg                                                                                            <  8:53381/120‑1 (MQ=255)
                                                                                        ccGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTgcgc                                                             >  2:216468/1‑149 (MQ=255)
                                                                                          gacgacCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCAc                                                           <  5:148289/149‑1 (MQ=255)
                                                                                                ccTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGcc                                                                                                                                    >  6:223350/1‑70 (MQ=255)
                                                                                                ccTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGcc                                                                                                                                    <  5:223350/70‑1 (MQ=255)
                                                                                                     cgGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgc                                                <  2:65082/149‑1 (MQ=255)
                                                                                                      ggATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTg                                                                              <  5:7051/118‑1 (MQ=255)
                                                                                                      ggATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTg                                                                              >  6:7051/1‑118 (MQ=255)
                                                                                                      ggATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgcg                                               <  6:182447/149‑1 (MQ=255)
                                                                                                       gATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgcgc                                              <  8:189727/149‑1 (MQ=255)
                                                                                                             aaaGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGt                                                                                     >  5:224979/1‑104 (MQ=255)
                                                                                                             aaaGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGt                                                                                     <  6:224979/104‑1 (MQ=255)
                                                                                                             aaaGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCacac                                        <  2:43335/149‑1 (MQ=255)
                                                                                                              aaGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACAcc                                       >  4:6225/1‑149 (MQ=255)
                                                                                                                          tATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggtg                           >  6:16837/1‑149 (MQ=255)
                                                                                                                            tCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATTCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACCCCGTCTTCAtggtggt                         >  4:22877/1‑148 (MQ=255)
                                                                                                                            tCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggtgg                          <  1:216468/148‑1 (MQ=255)
                                                                                                                            tCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGa                         <  3:168443/149‑1 (MQ=255)
                                                                                                                              cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgcgcg                                             >  4:3881/1‑127 (MQ=255)
                                                                                                                              cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgcgcg                                             <  3:3881/127‑1 (MQ=255)
                                                                                                                              cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACAcc                                       >  4:135057/1‑133 (MQ=255)
                                                                                                                              cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACAcc                                       <  3:135057/133‑1 (MQ=255)
                                                                                                                                aCATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGGTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGa                      >  4:235139/1‑148 (MQ=255)
                                                                                                                                 cATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggtgg                          <  6:251470/143‑1 (MQ=255)
                                                                                                                                 cATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggtgg                          >  5:251470/1‑143 (MQ=255)
                                                                                                                                   tCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTg                                                                                                                             >  4:75910/1‑42 (MQ=255)
                                                                                                                                      ttCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATCGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGCTGCTGGGCGCGCCCCCCGTCTTCATGGTGGATGACAGCGAc               >  8:267204/1‑149 (MQ=255)
                                                                                                                                         ccACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACaa                   >  7:55149/1‑142 (MQ=255)
                                                                                                                                         ccACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACaa                   <  8:55149/142‑1 (MQ=255)
                                                                                                                                         ccACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTAc            <  6:219423/149‑1 (MQ=255)
                                                                                                                                           aCCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACa                    <  6:86500/139‑1 (MQ=255)
                                                                                                                                           aCCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACa                    >  5:86500/1‑139 (MQ=255)
                                                                                                                                           aCCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCg          <  5:16837/149‑1 (MQ=255)
                                                                                                                                              gCTTGCCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACTGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCgg         >  8:276450/1‑147 (MQ=255)
                                                                                                                                                   gCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCGGCTGGCCt  <  1:135176/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCGGCTGGCCT  >  NC_002947/4479030‑4479325

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