Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 460,938 G→C 6.7% 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 6.7% 188.4 / 4.0 60I509M (ATC→ATG) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base C (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

AGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGC  >  NC_002947/460799‑461080
                                                                                                                                           |                                                                                                                                              
aGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATCCCCATCAGGGcc                                                                                                                                       <  4:108152/149‑1 (MQ=255)
   cTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCggg                                                                                                                                    <  7:87131/149‑1 (MQ=255)
   cTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCggg                                                                                                                                    <  7:133863/149‑1 (MQ=255)
                   ccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAg                                                                                                                     <  3:332195/148‑1 (MQ=255)
                    caCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCgggg                                                                                                                                   >  1:109014/1‑133 (MQ=255)
                    caCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCgggg                                                                                                                                   <  2:109014/133‑1 (MQ=255)
                         ggTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagc                                                                                                               <  3:88379/148‑1 (MQ=255)
                           tgtgCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagccag                                                                                                            >  3:70697/1‑149 (MQ=255)
                             tgCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagcca                                                                                                             >  1:156968/1‑146 (MQ=255)
                             tgCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagcca                                                                                                             <  2:156968/146‑1 (MQ=255)
                               cAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCaa                                                                                                         >  4:258687/1‑148 (MQ=255)
                               cAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCaa                                                                                                         <  3:258687/148‑1 (MQ=255)
                               cAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAAc                                                                                                        >  8:100791/1‑149 (MQ=255)
                               cAGCCACAAACTGTCGAGTTGCTGCTGAATCAGCATACAGGCAAACACCGCTGGGGGCCGGCTACCCCGTAACGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAAc                                                                                                        >  8:55405/1‑149 (MQ=255)
                                   cacaGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagccag                                                                                                            <  2:350923/141‑1 (MQ=255)
                                   cacaGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagccag                                                                                                            >  1:350923/1‑141 (MQ=255)
                                   cacaGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGc                                                                                                    <  4:70697/149‑1 (MQ=255)
                                    acaGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCa                                                                                                   <  2:92038/149‑1 (MQ=255)
                                     ttctatGTCGAGTTGCTGCTGGGTCAGCGGGCCGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAg                                                                                                  <  6:78048/144‑1 (MQ=255)
                                        gCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTg                                                                                               >  3:112854/1‑149 (MQ=255)
                                          tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                             <  4:203496/149‑1 (MQ=255)
                                              gAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGa                                                                                         <  7:100791/149‑1 (MQ=255)
                                              gAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGa                                                                                         <  5:310755/149‑1 (MQ=255)
                                                      gctgGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAg                                                                                        <  6:360362/142‑1 (MQ=255)
                                                      gctgGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAg                                                                                        >  5:360362/1‑142 (MQ=255)
                                                          ggTCAGCGTTCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggc                                                                               <  8:142313/147‑1 (MQ=255)
                                                          ggTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAg                                                                                                                                           <  8:230366/87‑1 (MQ=255)
                                                          ggTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAg                                                                                                                                           >  7:230366/1‑87 (MQ=255)
                                                              aGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                         <  4:112854/149‑1 (MQ=255)
                                                               gCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCt                                                                        >  6:206196/1‑149 (MQ=255)
                                                                         aggcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTg                                                              <  4:109460/149‑1 (MQ=255)
                                                                             aCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                         <  8:23123/134‑1 (MQ=255)
                                                                             aCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                         >  7:23123/1‑134 (MQ=255)
                                                                               cGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCgca                                                        <  2:296243/149‑1 (MQ=255)
                                                                                 cTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggc                                                                               <  7:282057/124‑1 (MQ=255)
                                                                                 cTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggc                                                                               >  8:282057/1‑124 (MQ=255)
                                                                                         aGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaagg                                                                                <  8:116767/115‑1 (MQ=255)
                                                                                           gCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                         >  2:288980/1‑120 (MQ=255)
                                                                                           gCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                         <  1:288980/120‑1 (MQ=255)
                                                                                            cTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTgg                                           >  6:274199/1‑149 (MQ=255)
                                                                                               ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                        <  5:206196/149‑1 (MQ=255)
                                                                                                  cGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAg                                                                                        <  3:243226/98‑1 (MQ=255)
                                                                                                  cGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAg                                                                                        >  4:243226/1‑98 (MQ=255)
                                                                                                  cGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTAGCAGGcc                                     >  3:29712/1‑149 (MQ=255)
                                                                                                      ccGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                 >  5:106802/1‑149 (MQ=255)
                                                                                                        gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                 >  4:114375/1‑147 (MQ=255)
                                                                                                          gTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                             >  2:122321/1‑149 (MQ=255)
                                                                                                               gCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAAc                        >  4:345854/1‑149 (MQ=255)
                                                                                                                  gTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCATGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                     <  8:211160/149‑1 (MQ=255)
                                                                                                                      cTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGc                           >  7:167551/1‑139 (MQ=255)
                                                                                                                      cTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGc                           <  8:167551/139‑1 (MQ=255)
                                                                                                                      cTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAg                  >  4:324416/1‑148 (MQ=255)
                                                                                                                       ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                >  7:62963/1‑149 (MQ=255)
                                                                                                                            cccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCCGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc           >  1:36807/4‑149 (MQ=255)
                                                                                                                            cccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                             <  5:276821/128‑1 (MQ=255)
                                                                                                                            cccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                             >  6:276821/4‑131 (MQ=255)
                                                                                                                            cAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTAGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc           <  4:29712/149‑1 (MQ=255)
                                                                                                                                      aTGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActgc  <  6:106802/148‑1 (MQ=255)
                                                                                                                                       tGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc           <  7:54745/138‑1 (MQ=255)
                                                                                                                                       tGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc           >  8:54745/1‑138 (MQ=255)
                                                                                                                                           |                                                                                                                                              
AGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGC  >  NC_002947/460799‑461080

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