Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,401,650 T→C 6.5% R215R (CGT→CGC)  queE → 7‑carboxy‑7‑deazaguanine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,401,6500TC 6.5% 176.9 / 2.8 61R215R (CGT→CGC) queE7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base T (32/25);  new base C (2/2);  total (34/27)
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.55e-01

GCCGAGCGTGCCGGTGAGGTGTTGTTCTCACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGACTACGGCCA  >  NC_002947/1401504‑1401786
                                                                                                                                                  |                                                                                                                                        
gCCGAGCGTGCCGGTGAGGTGTTGTTCTCACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGtt                                                                                                                                         <  2:194978/148‑1 (MQ=255)
     gCGTGCCGGTGAGGTGTTGTTCTCACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGttgatt                                                                                                                                     >  7:6600/1‑147 (MQ=255)
                aGGTGTTGTTCTCACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTa                                                                                                                         >  5:370159/1‑148 (MQ=255)
                       gttCTCACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACACCCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGa                                                                                                                  >  4:331475/1‑148 (MQ=255)
                           tcACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTgag                                                                                                                                  >  1:116850/1‑128 (MQ=255)
                           tcACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTgag                                                                                                                                  <  2:116850/128‑1 (MQ=255)
                           tcACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGa                                                                                                              >  3:312532/1‑148 (MQ=255)
                             aCCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGaaa                                                                                                            >  3:211243/1‑148 (MQ=255)
                             aCCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGaaa                                                                                                            >  2:332258/1‑148 (MQ=255)
                               ccAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGTAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACa                                                                                                                <  6:114772/142‑1 (MQ=255)
                               ccAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGTAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACa                                                                                                                >  5:114772/1‑142 (MQ=255)
                               ccAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAAc                                                                                                         <  4:211243/149‑1 (MQ=255)
                                 aGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAAcgc                                                                                                       <  8:6600/149‑1 (MQ=255)
                                     atcatcAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCGTGTACGCTTCCAACTGCAGATGAACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAg                                                                                                   >  8:297996/1‑149 (MQ=255)
                                     atcatcAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGCTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAg                                                                                                   >  1:16186/1‑149 (MQ=255)
                                     atcatcAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGCTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAg                                                                                                   >  1:22546/1‑149 (MQ=255)
                                             gTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCtg                                                                                           >  1:321703/1‑149 (MQ=255)
                                              tGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCtgt                                                                                          >  5:121449/1‑149 (MQ=255)
                                                     cAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCcgg                                                                                   >  1:365516/1‑149 (MQ=255)
                                                     cAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCcgg                                                                                   >  3:310721/1‑149 (MQ=255)
                                                       gTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCcggcg                                                                                 >  2:246738/1‑149 (MQ=255)
                                                       gTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCcggcg                                                                                 >  3:324018/1‑149 (MQ=255)
                                                            cTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTgagga                                                                                                                                >  3:193168/1‑97 (MQ=255)
                                                            cTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTgagga                                                                                                                                <  4:193168/97‑1 (MQ=255)
                                                               gCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTgg                                                                           <  7:240450/147‑1 (MQ=255)
                                                               gCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGAt                                                                         <  6:199534/149‑1 (MQ=255)
                                                                 tGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTc                                                                       >  2:239795/1‑149 (MQ=255)
                                                                   aCTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCg                                                                      <  1:239795/148‑1 (MQ=255)
                                                                   aCTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGCTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCg                                                                      <  2:16186/148‑1 (MQ=255)
                                                                   aCTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGCTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCg                                                                      <  2:22546/148‑1 (MQ=255)
                                                                    cTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGc                                                                    <  1:189228/149‑1 (MQ=255)
                                                                     tGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGc                                                                                                                              <  6:99584/90‑1 (MQ=255)
                                                                     tGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGc                                                                                                                              >  5:99584/1‑90 (MQ=255)
                                                                     tGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCtgtt                                                                                         <  2:380785/127‑1 (MQ=255)
                                                                     tGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCtgtt                                                                                         >  1:380785/1‑127 (MQ=255)
                                                                       gATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCAc                                                                 >  2:216051/1‑149 (MQ=255)
                                                                        aTCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCAcc                                                                >  6:214773/1‑149 (MQ=255)
                                                                           gTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGCAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTg                                                             <  1:332258/149‑1 (MQ=255)
                                                                            ttGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGGCTGGATTCGGCCACCGTgg                                                            <  6:121449/149‑1 (MQ=255)
                                                                            ttGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTgg                                                            <  6:370159/149‑1 (MQ=255)
                                                                               cTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTg                                                         <  4:324018/149‑1 (MQ=255)
                                                                               cTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTg                                                         >  3:291421/1‑149 (MQ=255)
                                                                                     aCCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATgg                                                   <  4:291421/149‑1 (MQ=255)
                                                                                          cctgTACGCTTCCTACTGCAGCTGCTCAAGTTGCTTTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATTACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCAGCGTGGTTGCCATGGCCAAg                                              <  4:310721/149‑1 (MQ=255)
                                                                                          cctgTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCaa                                               <  2:365516/148‑1 (MQ=255)
                                                                                            tgTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGc                                            >  7:326204/1‑149 (MQ=255)
                                                                                            tgTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGc                                            >  2:26298/1‑149 (MQ=255)
                                                                                                cGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGaa                                        <  5:400074/149‑1 (MQ=255)
                                                                                                cGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGaa                                        <  4:312532/149‑1 (MQ=255)
                                                                                                   ttCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGc                                     >  4:265771/1‑149 (MQ=255)
                                                                                                     ccAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCt                                    <  5:214773/148‑1 (MQ=255)
                                                                                                       aaCTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTGGCCATGGCCAAGGCCGAAGGCTa                                   >  6:268870/1‑147 (MQ=255)
                                                                                                          tGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAgct                              >  2:196773/1‑149 (MQ=255)
                                                                                                           gCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAgctg                             >  8:383651/1‑149 (MQ=255)
                                                                                                             aGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCg                                                                      <  6:27843/106‑1 (MQ=255)
                                                                                                             aGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCg                                                                      >  5:27843/1‑106 (MQ=255)
                                                                                                                  cacaAGTTGCTGTGGAACGACGAGCCCGGACGTTGAATGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACAcc                      <  3:331475/149‑1 (MQ=255)
                                                                                                                        ttGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGTGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGc                >  1:65962/1‑149 (MQ=255)
                                                                                                                             gtgGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGa           >  7:193355/1‑149 (MQ=255)
                                                                                                                                gAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGACTa        >  2:252189/1‑149 (MQ=255)
                                                                                                                                 aaCGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGACTAc       <  1:246738/149‑1 (MQ=255)
                                                                                                                                 aaCGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGACTAc       <  1:216051/149‑1 (MQ=255)
                                                                                                                                      cgaGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGACTACGGCCa  >  8:363704/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GCCGAGCGTGCCGGTGAGGTGTTGTTCTCACCCAGCCATCATCAGGTGAACGCCAGTGACCTGGCTGACTGGATCGTTGCTGACAACCTGCCTGTACGCTTCCAACTGCAGCTGCACAAGTTGCTGTGGAACGACGAGCCCGGACGTTGATTGAGGAGCACGTACACATGACAGAAAAACGCGCAGTAATCCTGTTGTCCGGCGGCCTGGATTCGGCCACCGTGGTTGCCATGGCCAAGGCCGAAGGCTACAGCTGCTACACCATGAGCTTCGACTACGGCCA  >  NC_002947/1401504‑1401786

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