Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 460,935 G→C 7.8% P510P (CCC→CCG) ‡ pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947460,9350GC 7.8% 201.0 / 6.7 64P510P (CCC→CCG) ‡pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base C (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

GGTATCCTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCA  >  NC_002947/460788‑461076
                                                                                                                                                   |                                                                                                                                             
ggTATCCTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATgg                                                                                                                                              <  2:253693/149‑1 (MQ=255)
      cTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGAt                                                                                                                                          >  5:154022/1‑147 (MQ=255)
       tCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCCGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAg                                                                                                                                       >  3:124508/1‑149 (MQ=255)
           aGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATCCCCATCAGGGcc                                                                                                                                   <  7:282129/149‑1 (MQ=255)
              cTGCGGCCTTGCCATGCCACTCGGCGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCGGGGGGCAGGCTACCCCGTCCCGCGTCTGGCAGTTGCTTGCGCAGGGCCCTGATGGGGATCAGGGCCggg                                                                                                                                <  7:277804/149‑1 (MQ=255)
                    ccatgccatgCAACTAGGTGTGAAGCCACAGGATGTCGAGTTGCTGCTTGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAgc                                                                                                                          >  3:154819/1‑149 (MQ=255)
                        gccatgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCa                                                                                                                            <  6:72873/143‑1 (MQ=255)
                        gccatgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCa                                                                                                                            >  5:72873/1‑143 (MQ=255)
                              ccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAgg                                                                                                                <  6:154022/149‑1 (MQ=255)
                              ccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAgg                                                                                                                <  4:124508/149‑1 (MQ=255)
                              ccaCTCGGGGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGGCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAg                                                                                                                 <  7:84290/148‑1 (MQ=255)
                               caCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAggg                                                                                                               >  5:290426/1‑149 (MQ=255)
                                 cTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCgggg                                                                                                                               <  6:86737/131‑1 (MQ=255)
                                 cTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCgggg                                                                                                                               >  5:86737/1‑131 (MQ=255)
                                    ggTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCAGCGCTGGGGGCAAGCAACCCCCGACCGCCTCAGGGAATTTCATTAGCAGGGCCCTGAAGGGGCTAAAGGCCGGGGGCCACGCCCTGCACGGcacc                                                                                                           >  1:92404/1‑146 (MQ=255)
                                             ccACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTgg                                                                                                 >  7:153474/1‑149 (MQ=255)
                                                 aGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAgg                                                                                             >  7:88352/1‑149 (MQ=255)
                                                  ggCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGt                                                                                            >  2:182835/1‑149 (MQ=255)
                                                     tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTTCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGc                                                                                                                   >  4:237099/1‑123 (MQ=255)
                                                     tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTTCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGc                                                                                                                   <  3:237099/123‑1 (MQ=255)
                                                     tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                         >  1:168873/1‑149 (MQ=255)
                                                     tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                         >  5:34893/1‑149 (MQ=255)
                                                       tCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCtt                                                                                       >  8:247774/1‑149 (MQ=255)
                                                          aGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAg                                                                                    <  2:92404/149‑1 (MQ=255)
                                                              gctgctGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCAt                                                                                 <  4:154819/148‑1 (MQ=255)
                                                                tgctgGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATc                                                                                >  1:300031/1‑147 (MQ=255)
                                                                  ctgGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGt                                                                                            <  6:286506/133‑1 (MQ=255)
                                                                  ctgGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGt                                                                                            >  5:286506/1‑133 (MQ=255)
                                                                    gggTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggc                                                                           >  4:203109/1‑148 (MQ=255)
                                                                       tCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTg                                                                                      >  6:298103/1‑134 (MQ=255)
                                                                       tCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTg                                                                                      <  5:298103/134‑1 (MQ=255)
                                                                       tCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaag                                                                       <  3:203109/149‑1 (MQ=255)
                                                                           cGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCtt                                                                   <  2:300031/149‑1 (MQ=255)
                                                                               caggcaggCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGtt                                                                <  8:153474/148‑1 (MQ=255)
                                                                               caggcaggCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTc                                                               <  1:182835/149‑1 (MQ=255)
                                                                                aggcaggcACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGtt                                                                >  2:340197/1‑147 (MQ=255)
                                                                                aggcaggcACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGtt                                                                <  1:340197/147‑1 (MQ=255)
                                                                                   caggcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACt                                                           >  4:52373/1‑149 (MQ=255)
                                                                                    aggcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCACTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCCCTg                                                          >  6:362687/1‑149 (MQ=255)
                                                                                      gcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGcc                                                        <  3:52373/149‑1 (MQ=255)
                                                                                      gcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGcc                                                        >  8:219816/1‑149 (MQ=255)
                                                                                       caCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCtt                                                                                       <  7:252262/117‑1 (MQ=255)
                                                                                       caCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCtt                                                                                       >  8:252262/1‑117 (MQ=255)
                                                                                                gggCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActg                                              <  6:34893/149‑1 (MQ=255)
                                                                                                        tACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaagg                                                                            >  2:135213/1‑111 (MQ=255)
                                                                                                                 ccGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                             <  7:219816/149‑1 (MQ=255)
                                                                                                                     gTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                         <  6:290426/149‑1 (MQ=255)
                                                                                                                     gTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                         <  2:168873/149‑1 (MQ=255)
                                                                                                                        tGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                     <  6:262175/102‑1 (MQ=255)
                                                                                                                        tGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                     >  5:262175/1‑102 (MQ=255)
                                                                                                                          gCAGTTGCTTGAGCAGGGCCCTGATGGGAATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                    >  1:197991/1‑149 (MQ=255)
                                                                                                                             gTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                 <  8:88352/149‑1 (MQ=255)
                                                                                                                                gCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa               <  2:197991/148‑1 (MQ=255)
                                                                                                                                    gAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa               <  2:294746/144‑1 (MQ=255)
                                                                                                                                    gAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa               >  1:294746/1‑144 (MQ=255)
                                                                                                                                        ccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAg                          <  6:184656/127‑1 (MQ=255)
                                                                                                                                        ccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAg                          >  5:184656/3‑129 (MQ=255)
                                                                                                                                          ggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGc                                                                                   >  6:284933/1‑70 (MQ=255)
                                                                                                                                          ggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGc                                                                                   <  5:284933/70‑1 (MQ=255)
                                                                                                                                          ggCCCTGATGGGGATCAGGGCCGGGGGCAGCCCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCgg    >  7:168257/1‑149 (MQ=255)
                                                                                                                                          ggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCa                                                             >  6:147270/1‑92 (MQ=255)
                                                                                                                                          ggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCa                                                             <  5:147270/92‑1 (MQ=255)
                                                                                                                                            cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa  >  7:324794/1‑149 (MQ=255)
                                                                                                                                            cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa  >  3:17820/1‑149 (MQ=255)
                                                                                                                                            cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa  >  2:50304/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GGTATCCTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCA  >  NC_002947/460788‑461076

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: