Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 460,938 G→C 7.9% I509M (ATC→ATG pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947460,9380GC 7.9% 198.6 / 6.9 64I509M (ATC→ATGpqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base C (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTA  >  NC_002947/460794‑461084
                                                                                                                                                |                                                                                                                                                  
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                                                                                                                                    >  5:72873/1‑143 (MQ=255)
                  gccatgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCa                                                                                                                                    <  6:72873/143‑1 (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                                                                                                                           <  3:237099/123‑1 (MQ=255)
                                               tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTTCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGc                                                                                                                           >  4:237099/1‑123 (MQ=255)
                                               tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                                 >  5:34893/1‑149 (MQ=255)
                                               tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                                 >  1:168873/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                                                                >  8:219816/1‑149 (MQ=255)
                                                                                gcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGcc                                                                <  3:52373/149‑1 (MQ=255)
                                                                                 caCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCtt                                                                                               >  8:252262/1‑117 (MQ=255)
                                                                                 caCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCtt                                                                                               <  7:252262/117‑1 (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          >  2:50304/1‑149 (MQ=255)
                                                                                                                                      cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa          >  3:17820/1‑149 (MQ=255)
                                                                                                                                      cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCa          >  7:324794/1‑149 (MQ=255)
                                                                                                                                              gggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTTTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTa  <  6:362287/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTA  >  NC_002947/460794‑461084

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: