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% 216.9 / 3.5 69Q25H (CAA→CACPP_4136LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base G (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCACC  >  NC_002947/4672787‑4673080
                                                                                                                                                   |                                                                                                                                                  
aaTGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACtt                                                                                                                                                   >  5:253245/1‑149 (MQ=255)
       ttGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGTTGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttg                                                                                                                                            >  1:251885/1‑148 (MQ=255)
           tCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtggt                                                                                                                                        <  2:251885/149‑6 (MQ=255)
              cagcagCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCATACCACAACAGCACTTGCTGTTCACGTGGc                                                                                                                                     <  1:21879/149‑1 (MQ=255)
              cagcagCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGc                                                                                                                                     <  1:27506/149‑1 (MQ=255)
                       cccGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCg                                                                                                                             >  7:320328/1‑148 (MQ=255)
                            gTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCgg                                                                                                                       <  3:337494/149‑1 (MQ=255)
                            gTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCgg                                                                                                                       <  8:53678/149‑1 (MQ=255)
                             tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      <  4:321308/149‑1 (MQ=255)
                             tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      <  6:87609/149‑1 (MQ=255)
                             tAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTAGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      <  2:176356/149‑1 (MQ=255)
                                 cACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGc                                                                                                                                     <  8:137357/130‑1 (MQ=255)
                                 cACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGc                                                                                                                                     >  7:137357/1‑130 (MQ=255)
                                 cACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGt                                                                                                                  >  1:290884/1‑149 (MQ=255)
                                    ttaaaaTTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTg                                                                                                                                       <  2:230158/122‑1 (MQ=255)
                                    ttCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCa                                                                                                                                 >  6:257837/1‑131 (MQ=255)
                                    ttCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCa                                                                                                                                 <  5:257837/131‑1 (MQ=255)
                                          tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                         <  2:290884/149‑1 (MQ=255)
                                              cgcAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtg                                                                                                     >  3:271661/1‑149 (MQ=255)
                                                   gATGTTGGATACGTGAAAATTCACCGTCGACACCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtgctgc                                                                                                 >  1:219966/1‑148 (MQ=255)
                                                    aTGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGa                                                                                               <  2:219966/149‑1 (MQ=255)
                                                              cGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCg                                                                                     >  2:157598/1‑149 (MQ=255)
                                                                tGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                   >  5:115266/1‑149 (MQ=255)
                                                                 gAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTg                                                                                                                                               >  6:211118/1‑88 (MQ=255)
                                                                 gAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTg                                                                                                                                               <  5:211118/88‑1 (MQ=255)
                                                                      ttCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCAGTTGCCAAAACCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTgtcagt                                                                             >  1:268692/1‑149 (MQ=255)
                                                                       tCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTgtcagt                                                                             >  4:337184/1‑148 (MQ=255)
                                                                           cGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTc                                                                                          >  6:61148/1‑131 (MQ=255)
                                                                           cGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTc                                                                                          <  5:61148/131‑1 (MQ=255)
                                                                            gTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGtt                                                                       >  2:147336/1‑149 (MQ=255)
                                                                               gACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTg                                                                     <  4:271661/148‑1 (MQ=255)
                                                                                    cTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTc                                                                >  4:124214/1‑148 (MQ=255)
                                                                                     ttGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCa                                                              >  4:108472/1‑149 (MQ=255)
                                                                                     ttGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCa                                                              >  5:17986/1‑149 (MQ=255)
                                                                                       gCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCAc                                                             <  1:157598/148‑1 (MQ=255)
                                                                                        cATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACgg                                                           <  2:268692/149‑1 (MQ=255)
                                                                                         aTTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACggg                                                          >  8:7658/1‑149 (MQ=255)
                                                                                           tCAAGGATCTGGCCGATCTCCCATGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGaa                                                                   <  6:222975/138‑1 (MQ=255)
                                                                                           tCAAGGATCTGGCCGATCTCCCATGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGaa                                                                   >  5:222975/1‑138 (MQ=255)
                                                                                             aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGttt                                                                      >  7:142749/1‑133 (MQ=255)
                                                                                             aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGAt                                                      >  3:7918/1‑149 (MQ=255)
                                                                                                    tGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCgg                                               >  4:332551/1‑149 (MQ=255)
                                                                                                      gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGc                                                 <  4:7918/145‑1 (MQ=255)
                                                                                                      gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCgggg                                             <  8:320328/149‑1 (MQ=255)
                                                                                                               ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCAtt                                         <  6:253245/144‑1 (MQ=255)
                                                                                                               ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCAt                                          <  3:124214/143‑1 (MQ=255)
                                                                                                               ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATtgctgc                                    >  1:246382/1‑149 (MQ=255)
                                                                                                                 aCGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                   >  3:223652/1‑100 (MQ=255)
                                                                                                                 aCGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                   <  4:223652/100‑1 (MQ=255)
                                                                                                                    aGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTAcc                               <  3:337184/149‑1 (MQ=255)
                                                                                                                         ttGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      <  3:192270/57‑1 (MQ=255)
                                                                                                                         ttGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      <  3:196874/57‑1 (MQ=255)
                                                                                                                         ttGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      >  4:192270/1‑57 (MQ=255)
                                                                                                                         ttGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCggg                                                                                                                      >  4:196874/1‑57 (MQ=255)
                                                                                                                           gCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                           <  6:17986/146‑1 (MQ=255)
                                                                                                                             cATAAGCACACCACAACAGCACTTGCTGTTCACGGGGCGGCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCAc                      <  4:56384/149‑1 (MQ=255)
                                                                                                                                 aGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCtt                   <  2:246382/148‑1 (MQ=255)
                                                                                                                                      acgtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTAc                 >  7:281146/6‑145 (MQ=255)
                                                                                                                                      acgtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTAc                 <  8:281146/140‑1 (MQ=255)
                                                                                                                                        cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca           >  6:220937/1‑149 (MQ=255)
                                                                                                                                         acaacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca           >  8:152075/1‑148 (MQ=255)
                                                                                                                                           aacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAg        <  3:332551/149‑1 (MQ=255)
                                                                                                                                            acaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGaa              >  7:4284/1‑142 (MQ=255)
                                                                                                                                            acaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGaa              <  8:4284/142‑1 (MQ=255)
                                                                                                                                            acaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAgg       <  7:152075/149‑1 (MQ=255)
                                                                                                                                             caGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                           <  6:56981/128‑1 (MQ=255)
                                                                                                                                             caGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCa                           >  5:56981/1‑128 (MQ=255)
                                                                                                                                             caGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGCGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc      <  6:115266/149‑1 (MQ=255)
                                                                                                                                              aGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc      >  1:350623/1‑148 (MQ=255)
                                                                                                                                                 aCTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAAGTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCAcc  >  2:17390/1‑149 (MQ=255)
                                                                                                                                                   ttGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc      <  7:131831/143‑1 (MQ=255)
                                                                                                                                                   ttGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc      >  8:131831/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTCATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCACC  >  NC_002947/4672787‑4673080

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: