Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT39.6% 62.6 / 32.3 53S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (4/28);  new base T (21/0);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 3.98e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NC_000913/4604207‑4604481
                                                                                                                                         |                                                                                                                                             
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttt                                                                                                                                             >  2:115896/1‑140 (MQ=255)
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctt                                                                                                                                             >  1:46550/1‑140 (MQ=255)
     tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                        <  1:307583/140‑1 (MQ=255)
     tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                        <  2:812847/140‑1 (MQ=255)
     tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                        <  1:519392/140‑1 (MQ=255)
     tCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                        <  2:595964/140‑1 (MQ=255)
       aTGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTGACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                      <  2:105723/140‑1 (MQ=255)
                 tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGGCATgcacagt                                                                                                                            >  2:926396/1‑134 (MQ=255)
                 tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                            <  2:191527/140‑1 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACTAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGc                                                                                                                                      >  2:294874/1‑127 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACTAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                      <  1:294874/127‑1 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAg                                                                                                                         >  2:343828/1‑140 (MQ=255)
                      aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGGgcg                                                                                                                       >  2:208312/1‑140 (MQ=255)
                          agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGATAAGTGAGCGC‑‑aga                                                                                                                 <  2:46550/140‑1 (MQ=255)
                          agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑aga                                                                                                                 <  1:906596/140‑1 (MQ=255)
                           gagTTTATTCAGAAAGTTTGCGTGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑agag                                                                                                                <  1:149952/140‑1 (MQ=255)
                           gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGTCCTGAGAAGGGAGCGC‑‑aggg                                                                                                                >  1:254997/1‑138 (MQ=255)
                           gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑agag                                                                                                                <  1:554260/140‑1 (MQ=255)
                              tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGAGAGCGC‑‑AGAGAgg                                                                                                             >  1:391679/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAggg                                                                                                            <  1:627766/140‑1 (MQ=255)
                                   tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTa                                                                                                        >  1:286548/1‑140 (MQ=255)
                                       aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGC‑‑agag                                                                                                                >  1:639294/1‑128 (MQ=255)
                                       aaaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑agag                                                                                                                <  2:639294/128‑1 (MQ=255)
                                        aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTg                                                                                                   <  1:947500/140‑1 (MQ=255)
                                        aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCATTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTg                                                                                                   <  2:375951/140‑1 (MQ=255)
                                                    aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACt                                                                                       <  2:286548/140‑1 (MQ=255)
                                                                     tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCAGAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTcc                                                                      >  2:362820/1‑140 (MQ=255)
                                                                     tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGGTTTcgc                                                                    >  2:351034/1‑140 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGAGAGCGC‑‑AGAGAGGGGGTATCTTTCCTTACTGAACTCGGGGTTAAAGTTTCCCGTAccc                                                              >  2:551707/1‑140 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCATGAGAAGTGAGCGG‑‑GGAGAGGGGTTTTCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTAccc                                                              >  2:676070/1‑140 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                              <  1:933953/140‑1 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                              <  1:884795/140‑1 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                              <  1:275051/140‑1 (MQ=255)
                                                                             ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGC‑‑AGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGGAccccg                                                            >  2:277614/1‑138 (MQ=255)
                                                                             ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAGGGC‑‑AGAGAGGGATTATTTTTCCTGACTTAACTGGCGGTTTAATTTTCGCGTAccccg                                                            >  2:731214/1‑138 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCAAAAGAAGTGAGCAC‑‑GGAGAAGGATTATCTTGCCTGGCTGAATCCGGGGTTTAGGTTACGCGTACCCCGCGt                                                         >  2:902217/1‑140 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGC‑‑AGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTTAGTTTTCGCGTACCCGGCgg                                                         >  2:24863/1‑139 (MQ=255)
                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGC‑‑AGAGAGGGGTTATTTTTCCTGACTGAACTGGCGGTTTAGTTTTCGCGTACCCGGCgg                                                         >  2:16331/1‑139 (MQ=255)
                                                                                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTCTTTCATGAGAAGCGAGTGC‑‑AGAAGGGGAGTTTCTTTCCTGACTGAAATGGCGATTAAATTTCCGCGCACCCGGCggt                                                        >  2:386218/1‑138 (MQ=255)
                                                                                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGGGGGC‑‑GGAGGGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAAGTTTCCCGGACCCgggggt                                                        >  2:414590/1‑136 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGC‑‑AGAGGGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCCGCGttt                                                       >  2:933953/1‑140 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGTGAGCGC‑‑AGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGgtt                                                       >  2:509586/1‑140 (MQ=255)
                                                                                     gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGGGC‑‑AGAGAGGGGTTATTTTTCCTGACTGAACTGGGGATTAAGTTTTCGCGTACCCCGCGGTTTGt                                                    >  2:41597/1‑140 (MQ=255)
                                                                                       gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGG‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                                  <  1:676070/140‑1 (MQ=255)
                                                                                        tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                                 <  1:731214/140‑1 (MQ=255)
                                                                                          aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                               <  1:926396/140‑1 (MQ=255)
                                                                                          aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                               <  1:41597/140‑1 (MQ=255)
                                                                                            ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                             <  1:208312/140‑1 (MQ=255)
                                                                                              cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                           <  1:509586/140‑1 (MQ=255)
                                                                                              cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                           <  1:351034/140‑1 (MQ=255)
                                                                                                ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGC‑‑AGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTTAGTTTCCGCGTACCCCGCGTTTTGTTATTTCGgatg                                         >  2:809006/1‑140 (MQ=255)
                                                                                                 gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGGGAAGGGAGGGC‑‑AGAGAGGGGTTATCTTGCCTTACTGAACTGGCGGTTAAGTTTTCGCGGACCCCGGGGTTAGTTTTTTCGgatga                                        >  2:100467/1‑140 (MQ=255)
                                                                                                 gATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                        <  1:551707/140‑1 (MQ=255)
                                                                                                  aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatgat                                       <  1:414590/140‑1 (MQ=255)
                                                                                                      cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGGTTAGTTATTTGGGATGATGata                                   >  2:652307/1‑140 (MQ=255)
                                                                                                           cgcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                              <  1:652307/140‑1 (MQ=255)
                                                                                                             cgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGG‑‑CGAGAGGGGTTATCTTGCCTGACTGAACTGGGGGTTTAGTTTTCGCGTACCCGGCGTTTTGTTTTTGCGGGTGATGATATTTAAgg                            >  2:817545/1‑139 (MQ=255)
                                                                                                                aTCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                           <  1:362820/140‑1 (MQ=255)
                                                                                                                     gTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGa                    <  1:809006/140‑1 (MQ=255)
                                                                                                                       ttGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGCGG‑‑AGAGAGGGGTTATTTTTGCTGACTGAACTTGGGGTTTAGTTTTCCCGTGCCCGGCGGTTTGTTTTTTCGGGTGATGATATTTAAGCTCGGCTGGtt                  >  2:478641/1‑140 (MQ=255)
                                                                                                                               gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt          <  1:343828/140‑1 (MQ=255)
                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  1:817545/140‑1 (MQ=255)
                                                                                                                                       ttctttttctGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  1:14418/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                             
ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGC‑‑AGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NC_000913/4604207‑4604481

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