Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913332,9910GT24.6% 148.1 / 19.1 57L207F (TTG→TTT) pdeLDNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (10/33);  new base T (14/0);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 2.55e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAG  >  NC_000913/332855‑333122
                                                                                                                                        |                                                                                                                                    
tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCa                                                                                                                                   >  1:436912/1‑140 (MQ=255)
 ggCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAg                                                                                                                                  <  2:881244/140‑1 (MQ=255)
   cGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCg                                                                                                                                <  1:352179/140‑1 (MQ=255)
         cATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑ttt                                                                                                                         <  2:1382261/140‑1 (MQ=255)
            ccGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGG‑TTTGaa                                                                                                                      >  2:186095/1‑140 (MQ=255)
            ccGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGA‑TTTGaa                                                                                                                      >  2:732226/1‑140 (MQ=255)
                tCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGA‑TTTGaaaaaa                                                                                                                  >  1:1167714/1‑139 (MQ=255)
                 cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAaga                                                                                                                 <  1:1077187/140‑1 (MQ=255)
                 cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAaga                                                                                                                 <  2:898305/140‑1 (MQ=255)
                 cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAaga                                                                                                                 <  1:1400373/140‑1 (MQ=255)
                 agacaggTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAaga                                                                                                                 <  1:484504/134‑1 (MQ=255)
                    aTTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGC‑TTTGAAAAaaaggg                                                                                                              >  1:1185402/1‑138 (MQ=255)
                    aTTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAgtg                                                                                                              <  2:1167714/140‑1 (MQ=255)
                    aTTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAGCGGGG‑TTTGTAAAaaaggg                                                                                                              >  2:263471/1‑138 (MQ=255)
                     ttGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAgtgt                                                                                                             <  2:582772/140‑1 (MQ=255)
                          aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGG‑TTTTAAAAAGAGGGTCTGaa                                                                                                        >  2:1400373/1‑140 (MQ=255)
                          aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGG‑TTTGAAAAAGAGTGTCTGaa                                                                                                        >  1:911871/1‑140 (MQ=255)
                          aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAaga                                                                                                                 >  1:122687/1‑131 (MQ=255)
                          aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAaga                                                                                                                 <  2:122687/131‑1 (MQ=255)
                          aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGaa                                                                                                        >  1:440079/1‑140 (MQ=255)
                          aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGaa                                                                                                        >  1:197260/1‑140 (MQ=255)
                           tAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAAc                                                                                                       <  2:901758/140‑1 (MQ=255)
                           tAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCCGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAAc                                                                                                       <  2:271807/140‑1 (MQ=255)
                              tGACCCGCCAACTGATGAAACAGTCTGCGGATATTCTGATGCCGGGAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAGCGGGATTTTGAAAAAG‑GGGTCTGGACCTg                                                                                                    >  1:1177737/1‑140 (MQ=255)
                              tGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGG‑TTTTAAAAAGAGGGTCTGAACCTg                                                                                                    >  2:1085173/1‑140 (MQ=255)
                              tGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTg                                                                                                    <  2:440079/140‑1 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGCCAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGCAAAAGAGTGTCTGAACCTGGTta                                                                                                <  1:1177892/140‑1 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGG‑TTTTGAAAAAAGAGGCTGAAACTGGttt                                                                                                >  2:713079/1‑139 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGG‑TTTTAAAAAAAGTGTCTGAAACTGGTta                                                                                                >  2:1177892/1‑140 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGA‑TTTGAAAAAAAGGGTCTGAACCTGGTta                                                                                                >  1:802705/1‑140 (MQ=255)
                                  ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAGCGGGA‑TTTTAAAAAAAGTGTTTGAAccctttct                                                                                                >  2:484504/1‑134 (MQ=255)
                                    gCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGG‑TTTGAAAAAGAGTGTCTGAACCTg                                                                                                    >  2:523452/1‑134 (MQ=255)
                                    gCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTg                                                                                                    <  1:523452/134‑1 (MQ=255)
                                    gCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTtaat                                                                                              <  2:1170737/140‑1 (MQ=255)
                                     ccAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGAGCGGGG‑TTTTAAAAAGAGGGTCTGAAACTGGTtaata                                                                                             >  2:104951/1‑140 (MQ=255)
                                            atgaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAgg                                                                                      >  1:53265/1‑140 (MQ=255)
                                                aaCAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAAc                                                                                  >  1:279656/1‑140 (MQ=255)
                                                   aGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGAt                                                                               <  2:831935/140‑1 (MQ=255)
                                                      cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGG‑TTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATaaa                                                                            >  1:1310221/1‑140 (MQ=255)
                                                      cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGCGCGGGG‑TTTTAAAAAAAGTGTCTGAACCTTGTTAATAAATTAGGTAACGATaaa                                                                            >  2:875875/1‑140 (MQ=255)
                                                         cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGAt                                                                               <  2:544726/134‑1 (MQ=255)
                                                         cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGAt                                                                               >  1:544726/1‑134 (MQ=255)
                                                         cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGG‑TTTTAAAAAAAGGGTCTGAAACTGGTTAATAAAATAGGTAAAGaaaaaaaa                                                                         >  2:698495/1‑138 (MQ=255)
                                                         cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGG‑TTTGAAAAAAAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATaaaaaa                                                                         >  2:1025055/1‑138 (MQ=255)
                                                               ttCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGt                                                                                     >  1:1417678/1‑122 (MQ=255)
                                                               ttCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGt                                                                                     <  2:1417678/122‑1 (MQ=255)
                                                               ttCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGG‑TTTGAAAAAGAGTGTCTGAACCTGGTTTATAAATTTGGGAAAGATAAAAACAAGCTg                                                                   >  2:610076/1‑140 (MQ=255)
                                                                   gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtc                                                               >  1:523549/1‑140 (MQ=255)
                                                                    aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtct                                                              <  2:197260/140‑1 (MQ=255)
                                                                    aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtct                                                              <  1:185171/140‑1 (MQ=255)
                                                                    aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtct                                                              <  1:713079/140‑1 (MQ=255)
                                                                      gCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGGAGCGGGG‑TTTTAAAAAAAGGGTCTGAACCCGGGTAATAAATTAGGTAACGATAAAATCAAAATGGTTCTCg                                                            >  2:1211483/1‑140 (MQ=255)
                                                                       ccGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGa                                                           >  1:779314/1‑140 (MQ=255)
                                                                         ggTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGc                                                         <  2:53265/140‑1 (MQ=255)
                                                                          gTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGA‑TTTGAAAAAAAGGGTCTGAACCTGGGTAATAAATTAGGTAACGATAAAAACAAACTGGTTTTCGAggt                                                        >  2:1190435/1‑138 (MQ=255)
                                                                           tAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGG‑TTTTGAAAAAAGGGTCTGAACCTGGTTTATTAATTAGGTAACGATAAAAACAAACTGGGTCTCGAGGTa                                                       >  2:356080/1‑140 (MQ=255)
                                                                              aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACg                                                    <  1:591924/140‑1 (MQ=255)
                                                                                   ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACg                                               <  1:698495/140‑1 (MQ=255)
                                                                                         ccGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGTGGGG‑TTTGAAAAAAAAAGTTTGAACCTGGTTTATAAATTAGGGAACGCTAAAATTAAGGTGGGTCTTGAGCCAAAGGAAAGTAAccc                                         >  2:366579/1‑140 (MQ=255)
                                                                                           ggACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTa                                       <  2:911871/140‑1 (MQ=255)
                                                                                                         ttGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCa                         <  1:1287307/140‑1 (MQ=255)
                                                                                                           ggCATCAACGTCTCGGCGGGTTTTTTTTTGGCAGCGGGA‑TTTTTAAAAGAGAGTCTGAACCTGGGTAATAAATTAGGTAAAGATAAAATTAAGCTGGGTCTCGAGCTAAAGGAAAGTAACCCTATTCCGGTTAAGcccaa                       >  2:1287307/1‑137 (MQ=255)
                                                                                                           ggCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGA‑TTTTAAAAAAAGGGTCTGAACCTGGGTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGa                       >  2:1356133/1‑140 (MQ=255)
                                                                                                            gCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGGGTCTGAAACTGGGTTATAAATTAGGTAACGATAAAATCAAGCTGGTTTTCGAGCTAACGGAAAGTAACCCTATTCCGGTAACGCCaaaa                      >  2:536443/1‑140 (MQ=255)
                                                                                                               tCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGcc                   <  2:523549/140‑1 (MQ=255)
                                                                                                                   cGTCTCGGCGGGTTGTTTTTTTGCAGCGGGG‑TTTGAAAAAGAGGGGCTGAACCCGGGTTATAAATTAGGGAACGATAAAATCAAGCTGGTTCTCGAGCTAAAGGAACGTAAACCTATTCCGGTAAAACCAGAAGCCagag               >  2:171427/1‑140 (MQ=255)
                                                                                                                   cGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag               <  1:536443/140‑1 (MQ=255)
                                                                                                                     tctcGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCg             <  1:732226/140‑1 (MQ=255)
                                                                                                                      ctcGGCGGGTTGTTTTTTTGCAGCGGGG‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGa            >  1:319111/1‑140 (MQ=255)
                                                                                                                        cggcggGGTGTTTTTTTGCAGCGGGG‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAAAGATAAAATCAAGCTGGGTCTCGAGCTAACGGAACGTAACCCCATTTCGGTAACGCCAGAAGGCAGAGCGata          >  2:628005/1‑140 (MQ=255)
                                                                                                                        cggcggGGTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGGTCTCGAGCTAACGGAACGTAACCCTATTCCGGGAACGCCAGAAACCAGAGCGata          >  2:647031/1‑140 (MQ=255)
                                                                                                                             gggtgttTTTTTGCAGCGGGG‑TTTGAAAAAAAGGGGCTGAACCTGGTTTATAAATTAAGTAAAGATAAAAACAACCTGGGTCTCGCGCTAACCGAACGTAACCCTATTCCGGTAACGCCACAAGCCAGAGCGGTATTTGa     >  2:348707/4‑140 (MQ=255)
                                                                                                                             ggTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa     <  1:356080/140‑1 (MQ=255)
                                                                                                                             ggTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa     <  1:1085173/140‑1 (MQ=255)
                                                                                                                                tgttTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAg  <  2:802705/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGA‑TTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAG  >  NC_000913/332855‑333122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑