Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,672,934 T→G 5.8% Q25H (CAA→CAC)  PP_4136 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,672,9340TG 5.8% 200.0 / 2.3 64Q25H (CAA→CAC) PP_4136LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCAC  >  NC_002947/4672791‑4673079
                                                                                                                                               |                                                                                                                                                 
gCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTg                                                                                                                                              <  5:125906/149‑1 (MQ=255)
 ccTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGt                                                                                                                                             >  7:27089/1‑149 (MQ=255)
  cTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGtt                                                                                                                                            >  4:108948/1‑149 (MQ=255)
     gATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTc                                                                                                                                           <  1:352786/147‑1 (MQ=255)
       tCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtgga                                                                                                                                       <  1:22890/149‑6 (MQ=255)
                    ccGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGaa                                                                                                                          >  2:222361/1‑149 (MQ=255)
                    ccGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGaa                                                                                                                          <  1:165633/149‑1 (MQ=255)
                        gTAGCCACTTCAAATTTGCGCCGGAAGTTGCATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACTCCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCgg                                                                                                                      <  8:292063/149‑1 (MQ=255)
                         tAGGCAGTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                     <  8:27089/149‑1 (MQ=255)
                         tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                     <  2:258270/149‑1 (MQ=255)
                               cTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCa                                                                                                                                          >  6:237732/1‑122 (MQ=255)
                               cTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCa                                                                                                                                          <  5:237732/122‑1 (MQ=255)
                                ttCAAATTTGCGCGGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCg                                                                                                              >  2:191341/1‑149 (MQ=255)
                                  cAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCGCGTGGGGTCAGGCGCACCGGGGCGTGCGcc                                                                                                            >  2:325197/1‑149 (MQ=255)
                                      tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                        <  3:108948/149‑1 (MQ=255)
                                          cgcAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtg                                                                                                    >  4:249886/1‑149 (MQ=255)
                                                        tACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAg                                                                                                          <  8:172074/129‑1 (MQ=255)
                                                        tACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAg                                                                                                          >  7:172074/1‑129 (MQ=255)
                                                           gTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGt                                                                                   >  2:13358/1‑149 (MQ=255)
                                                            tGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                  <  1:191341/149‑1 (MQ=255)
                                                                    cACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCa                                                                                                                                >  1:60275/1‑95 (MQ=255)
                                                                    cACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCa                                                                                                                                <  2:60275/95‑1 (MQ=255)
                                                                     aCCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCAt                                                                         >  6:290144/1‑149 (MQ=255)
                                                                       cGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGt                                                                       <  3:257765/149‑1 (MQ=255)
                                                                       cGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGt                                                                       >  8:42492/1‑149 (MQ=255)
                                                                         tCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGttt                                                                     >  7:148970/1‑149 (MQ=255)
                                                                         tCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGttt                                                                     <  8:148970/149‑1 (MQ=255)
                                                                                     aTTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                        <  2:314295/102‑1 (MQ=255)
                                                                                     aTTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                        >  1:314295/1‑102 (MQ=255)
                                                                                      ttCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGt                                                        <  1:222361/149‑1 (MQ=255)
                                                                                      ttCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGt                                                        >  8:21420/1‑149 (MQ=255)
                                                                                       tCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGa                                                                                       <  4:21728/117‑1 (MQ=255)
                                                                                       tCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGa                                                                                       >  3:21728/1‑117 (MQ=255)
                                                                                         aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGtt                                                                                                     >  4:96972/1‑101 (MQ=255)
                                                                                         aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGtt                                                                                                     <  3:96972/101‑1 (MQ=255)
                                                                                         aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGAt                                                     >  4:295551/1‑149 (MQ=255)
                                                                                                 ggCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                             >  7:81290/1‑149 (MQ=255)
                                                                                                 ggCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                             >  7:76459/1‑149 (MQ=255)
                                                                                                  gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCgggg                                            >  6:38713/1‑149 (MQ=255)
                                                                                                  gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                             <  8:76459/148‑1 (MQ=255)
                                                                                                  gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                             <  3:249886/148‑1 (MQ=255)
                                                                                                  gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCggg                                             <  8:81290/148‑1 (MQ=255)
                                                                                                     gATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCAt                                         >  5:210580/1‑149 (MQ=255)
                                                                                                           ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATtgct                                     <  7:42492/147‑1 (MQ=255)
                                                                                                               gAGCTCTTGCAATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCAGCCAGTGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGCAGAATTGCTGCCTa                                >  1:304855/1‑148 (MQ=255)
                                                                                                                  ctctTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTg                            >  7:105099/1‑149 (MQ=255)
                                                                                                                   tctTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTg                            >  5:281922/1‑148 (MQ=255)
                                                                                                                    ctTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCgg                                              >  7:98117/1‑129 (MQ=255)
                                                                                                                    ctTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCgg                                              <  8:98117/129‑1 (MQ=255)
                                                                                                                       gCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                          <  5:38713/146‑1 (MQ=255)
                                                                                                                         cATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCAc                     <  3:295551/149‑1 (MQ=255)
                                                                                                                          aTAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCAcc                    >  6:10912/1‑149 (MQ=255)
                                                                                                                                  accacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCAt                                         <  3:170553/120‑1 (MQ=255)
                                                                                                                                  accacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCAt                                         >  4:170553/1‑120 (MQ=255)
                                                                                                                                    gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGGGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca          >  8:247385/4‑149 (MQ=255)
                                                                                                                                    gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGt                        >  6:103985/4‑135 (MQ=255)
                                                                                                                                    gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGt                        <  5:103985/132‑1 (MQ=255)
                                                                                                                                     acaacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGcac         <  5:290144/149‑1 (MQ=255)
                                                                                                                                      caacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTAc                               <  7:208120/126‑1 (MQ=255)
                                                                                                                                      caacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTAc                               >  8:208120/1‑126 (MQ=255)
                                                                                                                                        acaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAgg      <  5:10912/149‑1 (MQ=255)
                                                                                                                                            cACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                          >  3:238361/1‑125 (MQ=255)
                                                                                                                                            cACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                          <  4:238361/125‑1 (MQ=255)
                                                                                                                                            cACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCAc  <  6:210580/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCAC  >  NC_002947/4672791‑4673079

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