Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 812,056 G→A 7.9% L67L (CTG→TTG)  PP_0698 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947812,0560GA 7.9% 214.0 / 7.4 77L67L (CTG→TTG) PP_0698LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base A (3/3);  total (34/43)
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.93e-01

ATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCT  >  NC_002947/811910‑812199
                                                                                                                                                  |                                                                                                                                                  
aTGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGGTTTGCGGCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGc                                                                                                                                                 >  2:457994/1‑150 (MQ=255)
     aGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGACGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGAc                                                                                                                                            <  2:1596016/150‑1 (MQ=255)
         cgGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccacccag                                                                                                                                        <  1:1168325/150‑4 (MQ=255)
          ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACc‑‑ggcg                                                                                                                                     <  2:1359059/150‑1 (MQ=255)
          ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTGCAAGGCAcc                                                                                                                                              >  1:1378866/1‑142 (MQ=255)
          ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTGCAAGGCAcc                                                                                                                                              <  2:1378866/143‑2 (MQ=255)
               gcgcGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCgggg                                                                                                                                  >  1:190708/1‑148 (MQ=255)
               gcgcGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCgggg                                                                                                                                  <  2:190708/148‑1 (MQ=255)
               gcgcGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGt                                                                                                                                <  1:1239517/150‑1 (MQ=255)
                cgcgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTg                                                                                                                               <  2:1046089/150‑1 (MQ=255)
                     gATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGaa                                                                                                                         >  1:1407076/1‑150 (MQ=255)
                           ccagcGTATGGCGATGGGCGGTTCGCAAGACAGCACCAGCACATTGGCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGc                                                                                                                   <  1:181272/150‑1 (MQ=255)
                           ccagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGc                                                                                                                   >  2:463935/1‑150 (MQ=255)
                              gcgcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCg                                                                                                                <  2:1038179/150‑1 (MQ=255)
                                gcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccacccag                                                                                                                                        >  1:1430949/1‑124 (MQ=255)
                                gcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccacccag                                                                                                                                        <  2:1430949/127‑4 (MQ=255)
                                  aTGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTg                                                                                                            >  1:1079589/1‑150 (MQ=255)
                                        aTGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTTCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGt                                                                                                      >  2:404322/1‑150 (MQ=255)
                                             cGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGaa                                                                                                 <  2:291642/150‑1 (MQ=255)
                                              ggTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                    >  2:738979/1‑146 (MQ=255)
                                              ggTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                    <  1:738979/146‑1 (MQ=255)
                                                  cGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGACGCTCGGCTGGGTGGTGGCTTGCATGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGtt                                                                                            <  1:610047/150‑1 (MQ=255)
                                                   gCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGcc                                                                                                                  <  2:1471975/127‑1 (MQ=255)
                                                   gCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGcc                                                                                                                  >  1:1471975/1‑127 (MQ=255)
                                                    cAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGCAAAATTGCCGGTTGCGGCGTTCGAACAGTTcc                                                                                          <  1:457994/150‑1 (MQ=255)
                                                     aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGTTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                         >  2:1194242/1‑150 (MQ=255)
                                                     aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCa                                                                                         >  2:1074823/1‑150 (MQ=255)
                                                      aGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCAc                                                                                        <  2:1407076/150‑1 (MQ=255)
                                                          agcaccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATgg                                                                                                                    <  2:1460955/118‑1 (MQ=255)
                                                          agcaccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATgg                                                                                                                    >  1:1460955/1‑118 (MQ=255)
                                                           gcaccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCaa                                                                                   <  2:1079589/150‑1 (MQ=255)
                                                            caccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAg                                                                                  <  1:893872/150‑1 (MQ=255)
                                                             tccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAgg                                                                                 <  1:980903/149‑1 (MQ=255)
                                                             accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGt                                                                                                                                <  2:1572195/104‑1 (MQ=255)
                                                             accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGt                                                                                                                                >  1:1572195/1‑104 (MQ=255)
                                                             accagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAgg                                                                                 <  1:1256499/150‑1 (MQ=255)
                                                               cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                               >  2:596621/1‑150 (MQ=255)
                                                               cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                               >  1:158568/1‑150 (MQ=255)
                                                                 gcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGccc                                                                                                                 >  2:29010/1‑114 (MQ=255)
                                                                 gcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTCCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGccc                                                                                                                 <  1:29010/114‑1 (MQ=255)
                                                                  cacaTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttc                                                                            <  1:748708/150‑1 (MQ=255)
                                                                     aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                         >  1:891053/1‑150 (MQ=255)
                                                                     aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                         >  1:618450/1‑150 (MQ=255)
                                                                     aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCTCTTCGGCACGACC‑‑GCTGCGTGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                         <  1:609620/150‑1 (MQ=255)
                                                                      ttGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATCGCCCGCTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTcc                                                                        <  1:694460/150‑1 (MQ=255)
                                                                          cgcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagca                                                                    <  1:404322/150‑1 (MQ=255)
                                                                                    ttgCTGCAGGCGCTGCACGGCGCCTTCGAGTTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTg                                                          <  1:1194242/150‑1 (MQ=255)
                                                                                      gctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggc                                                        >  2:1403798/1‑150 (MQ=255)
                                                                                      gctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggc                                                        <  2:891053/150‑1 (MQ=255)
                                                                                      gctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggc                                                        <  1:596621/150‑1 (MQ=255)
                                                                                      gctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGAGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggc                                                        <  1:999243/150‑1 (MQ=255)
                                                                                       ctgcAGGCGCTGCACGGCGCCTTCGAGGTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcg                                                       <  1:695685/150‑1 (MQ=255)
                                                                                       ctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcg                                                       >  2:848295/1‑150 (MQ=255)
                                                                                        tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCca                                                                       <  2:360578/133‑1 (MQ=255)
                                                                                        tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCca                                                                       >  1:360578/1‑133 (MQ=255)
                                                                                             gcgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGAc                                                 <  2:1646468/150‑1 (MQ=255)
                                                                                             gcgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGATGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTc                                                                                           >  2:342656/1‑108 (MQ=255)
                                                                                             gcgcTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGATGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTc                                                                                           <  1:342656/108‑1 (MQ=255)
                                                                                                     cGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgt                                         <  1:1074823/150‑1 (MQ=255)
                                                                                                      ggCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtg                                        <  1:848295/150‑1 (MQ=255)
                                                                                                           cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGaaaa                                                                                                                       >  1:267535/1‑66 (MQ=255)
                                                                                                           cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGaaaa                                                                                                                       <  2:267535/66‑1 (MQ=255)
                                                                                                            ttCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGAc                                  >  2:800862/1‑150 (MQ=255)
                                                                                                                 gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAg                                                                                              <  1:546127/85‑1 (MQ=255)
                                                                                                                 gCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAg                                                                                              >  2:546127/1‑85 (MQ=255)
                                                                                                                        gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcag                                                                      >  1:499292/1‑102 (MQ=255)
                                                                                                                        gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcag                                                                      <  2:499292/102‑1 (MQ=255)
                                                                                                                        gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTc                      >  2:491655/1‑150 (MQ=255)
                                                                                                                         aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCca                                                                       <  2:1358729/100‑1 (MQ=255)
                                                                                                                         aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCca                                                                       >  1:1358729/1‑100 (MQ=255)
                                                                                                                           gATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTTAGTGACGTGCAAGTGCTcggc                   >  1:961738/1‑150 (MQ=255)
                                                                                                                              cccgccggtcgtGCCCTGCAAGGCACCACCCAGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTcggcg                  >  1:1632223/7‑150 (MQ=38)
                                                                                                                                  tGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCAGGGGTGANGAAAATGGCNCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGgcgcgc            >  1:723702/1‑150 (MQ=255)
                                                                                                                                   gggTGGTGCCTTGGAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCa           <  1:463935/150‑1 (MQ=255)
                                                                                                                                            cTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGctct  >  2:1342040/1‑150 (MQ=255)
                                                                                                                                               gCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCAc                                                                                        <  2:509970/61‑1 (MQ=255)
                                                                                                                                               gCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCAc                                                                                        >  1:509970/1‑61 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACC‑‑GGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCT  >  NC_002947/811910‑812199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: