Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 931,497 G→T 9.3% L4709I (CTC→ATC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947931,4970GT 9.3% 182.4 / 4.1 60L4709I (CTC→ATC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (24/30);  new base T (2/4);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCT  >  NC_002947/931351‑931634
                                                                                                                                                   |                                                                                                                                         
ggTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAggg                                                                                                                                         >  6:118510/1‑149 (MQ=255)
 gTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                        <  4:181578/149‑1 (MQ=255)
  tGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTTAAGTGTCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTCTCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                        <  1:96541/149‑1 (MQ=255)
  tGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCAATGGCATCGACGGTGGCGCTGCCGATTTCCTAGCCT‑GCCGGGTTGGTCAAGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                        >  8:25705/1‑148 (MQ=255)
   gAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTgt                                                                                                                                      >  3:128868/1‑148 (MQ=255)
       cTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGacac                                                                                                                                   >  1:69395/1‑148 (MQ=255)
       cTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACAcc                                                                                                                                  <  3:324540/149‑1 (MQ=255)
           gcagcTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTGTCACCCTCa                                                                                                                              <  4:128868/149‑1 (MQ=255)
             agcTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCAcc                                                                                                                            <  3:13699/149‑1 (MQ=255)
                tCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTGGCCGGTGAGGGTGACACCATCACCgg                                                                                                                          >  6:268415/1‑148 (MQ=255)
                     gTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACgg                                                                                                                    <  5:118510/149‑1 (MQ=255)
                      ttGCCCCTTGCCGGGGTCAGCGTGTTGGTTATGTGGCCACTTGCATCGACGGTTGCGCTGCCGATTTCCTCGCCT‑GCCGGTTTTGTCACTGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGTTGTTGACACCATCACCGGAAACGGc                                                                                                                   <  7:232850/149‑1 (MQ=255)
                      ttGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGc                                                                                                                   <  3:107178/149‑1 (MQ=255)
                          ccAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTTTCCGGTGATGGTGTCACCCTCACCggccgcgg                                                                                                                    >  4:58278/1‑138 (MQ=38)
                          ccAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTTTCCGGTGATGGTGTCACCCTCACCggccgcgg                                                                                                                    <  3:58278/144‑7 (MQ=38)
                            aGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAggg                                                                                                                                         <  8:59464/121‑1 (MQ=255)
                            aGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAggg                                                                                                                                         >  7:59464/1‑121 (MQ=255)
                                 cGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAAc                                                                                                                      >  6:204703/1‑135 (MQ=255)
                                 cGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAAc                                                                                                                      <  5:204703/135‑1 (MQ=255)
                                      tCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGcc                                                                                                   >  4:223305/1‑149 (MQ=255)
                                          cGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATg                                                                                               >  3:289795/1‑149 (MQ=255)
                                           gTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATgg                                                                                              <  1:91144/149‑1 (MQ=255)
                                           gTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATgg                                                                                              <  4:289795/149‑1 (MQ=255)
                                            ttaTGTTGAAGTGGCCACTGGCATCGACGGTGGCGATGACGATTTCATCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGAGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATgg                                                                                              <  8:243517/146‑1 (MQ=255)
                                              aTGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCgg                                                                                           >  5:106662/1‑149 (MQ=255)
                                                 gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATc                                                                                                                               >  6:135648/1‑110 (MQ=255)
                                                 gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATc                                                                                                                               <  5:135648/110‑1 (MQ=255)
                                                 gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgtgt                                                                                        >  1:144889/1‑149 (MQ=255)
                                                   gAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTgg                                                                                      >  5:101491/1‑149 (MQ=255)
                                                   gAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTgg                                                                                      >  8:169156/1‑149 (MQ=255)
                                                     aGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTa                                                                                    <  6:201658/149‑1 (MQ=255)
                                                       tGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATc                                                                                  >  4:141476/1‑149 (MQ=255)
                                                       tGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATc                                                                                  >  4:142099/1‑149 (MQ=255)
                                                                 cATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTg                                                                        <  7:252930/149‑1 (MQ=255)
                                                                       cGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCAGCGGCCGGTGATGGTGACACCATCACCGGAAACGGCCAGTTCGCCTGGTGCCTATGGCGGTGTGGTATCGGGTTCGGTGAt                                                                      <  4:349190/145‑1 (MQ=255)
                                                                            gcgcTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggc                                                             <  6:101491/149‑1 (MQ=255)
                                                                            gcgcTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTgg                                                              <  1:204860/148‑1 (MQ=255)
                                                                            gcgcTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGCGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggc                                                             <  4:140382/149‑1 (MQ=255)
                                                                              gcTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                           <  1:108842/149‑1 (MQ=255)
                                                                               cTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                           >  6:278175/1‑148 (MQ=255)
                                                                                tGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGac                                                                                                                                     >  5:103481/1‑73 (MQ=255)
                                                                                tGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGac                                                                                                                                     <  6:103481/73‑1 (MQ=255)
                                                                                tGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                           >  6:376585/1‑147 (MQ=255)
                                                                                  ccGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGTGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCa                                                         <  3:223305/147‑1 (MQ=255)
                                                                                  ccGATGTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGATGGCGGCAGa                                                       >  7:336735/1‑149 (MQ=255)
                                                                                          cTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGtt                                                 >  8:23021/1‑147 (MQ=255)
                                                                                          cTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGtt                                                 <  7:23021/147‑1 (MQ=255)
                                                                                          cTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtgc                                               >  6:374677/1‑149 (MQ=255)
                                                                                           tCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGAc                                                      <  3:16305/141‑1 (MQ=255)
                                                                                           tCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGAc                                                      >  4:16305/1‑141 (MQ=255)
                                                                                                    cGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCAt                                      <  3:142099/149‑1 (MQ=255)
                                                                                                            tCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAg                               >  5:5419/1‑148 (MQ=255)
                                                                                                            tCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGc                              >  5:208771/1‑149 (MQ=255)
                                                                                                               cGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGcacc                           <  5:376585/149‑1 (MQ=255)
                                                                                                                         tGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCt                 <  7:25705/149‑1 (MQ=255)
                                                                                                                          gcgcCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTg                >  7:175577/1‑149 (MQ=255)
                                                                                                                           cgcCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTg                <  6:208771/148‑1 (MQ=255)
                                                                                                                            gcCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgtt              <  6:5419/149‑1 (MQ=255)
                                                                                                                                aCTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggcc          <  2:69395/149‑1 (MQ=255)
                                                                                                                                    cgccgcGGCCGGTGAGTGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgtt       <  4:261218/148‑1 (MQ=255)
                                                                                                                                     gccgcGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgttg      <  6:106662/148‑1 (MQ=255)
                                                                                                                                        gcgGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCt  <  5:374677/149‑1 (MQ=255)
                                                                                                                                        gcgGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCt  >  8:256004/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCT‑GCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCT  >  NC_002947/931351‑931634

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