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

ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATT  >  NC_000913/4604207‑4604473
                                                                                                                                          |                                                                                                                                  
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGTGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttctt                                                                                                                                  >  1:985378/1‑140 (MQ=255)
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttctt                                                                                                                                  >  1:1097545/1‑140 (MQ=255)
aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATtttttctt                                                                                                                                  >  2:633555/1‑140 (MQ=255)
      cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CAttttttttttttttt                                                                                                                            >  2:952338/1‑137 (MQ=255)
              gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATGCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATgaga                                                                                                                    <  2:899883/140‑1 (MQ=255)
              gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATgaga                                                                                                                    <  2:148204/140‑1 (MQ=255)
                 tCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCCTGAGAAgg                                                                                                                 >  1:575685/1‑139 (MQ=255)
                    gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGAAGTGA‑g                                                                                                             >  2:484572/1‑140 (MQ=255)
                      aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTCTGCCATGAGAAGt                                                                                                                 >  1:664590/1‑135 (MQ=255)
                      aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                 <  2:664590/135‑1 (MQ=255)
                      aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑gcg                                                                                                           <  1:12879/140‑1 (MQ=255)
                        tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCa                                                                                                         <  2:979971/140‑1 (MQ=255)
                        tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCa                                                                                                         <  2:277960/140‑1 (MQ=255)
                        tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCa                                                                                                         <  1:820095/140‑1 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTTCATGAAAAGGGA‑gcgggg                                                                                                        >  2:820095/1‑137 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTTAGAAGAGA‑GCGCag                                                                                                        >  2:1097765/1‑140 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCag                                                                                                        <  2:719559/140‑1 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCag                                                                                                        <  2:1111645/140‑1 (MQ=255)
                         gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGC‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCag                                                                                                        <  1:382335/140‑1 (MQ=255)
                          agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCTGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGTGA‑GGGCaga                                                                                                       >  2:1119944/1‑140 (MQ=255)
                          agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCaga                                                                                                       <  1:1119944/140‑1 (MQ=255)
                           gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCagag                                                                                                      <  1:484572/140‑1 (MQ=255)
                            agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑g                                                                                                             >  1:805376/1‑132 (MQ=255)
                            agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑g                                                                                                             <  2:805376/132‑1 (MQ=255)
                              tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTTCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAAGGA‑GGGCAGAGAgg                                                                                                   >  2:361835/1‑140 (MQ=255)
                              tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGGGA‑GGGCAGAGAgg                                                                                                   >  1:70468/1‑140 (MQ=255)
                              tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCA‑TTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGG‑CCTTTTTTTTTTTTTTTCCCTGAGAAAAGG‑GGGCCGAGAggg                                                                                                  >  1:566980/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAAGGA‑GCGCAGAGAggg                                                                                                  >  1:909959/1‑140 (MQ=255)
                               ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTCTTCCCTGAGAAGTGG‑GGGCagggaggg                                                                                                  >  2:34510/1‑140 (MQ=255)
                                tATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGATTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATTAGAAGTGA‑GCGCAGAGAgggg                                                                                                 >  2:1049179/1‑139 (MQ=255)
                                 aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CCTTTTTTTTTTTTTTTTCTTGAGAAGGGG‑GGGCAGgggggggg                                                                                                >  2:596864/1‑138 (MQ=255)
                                    cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTAt                                                                                             <  1:633555/140‑1 (MQ=255)
                                           tttGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCct                                                                                      <  1:475555/140‑1 (MQ=255)
                                                          tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAt                                                                       >  2:491511/1‑140 (MQ=255)
                                                                cTCCATACCTATACCAGTCAGGAGTCCTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                                                                 <  1:842672/140‑1 (MQ=255)
                                                                   cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCCTGAGAAGGGA‑GCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGttttc                                                              >  2:54434/1‑138 (MQ=255)
                                                                     tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑TCTTTTTTTTTTTTTTACCCTGAGAAGGGA‑ACGCAACGAGGGATTTTCTTGCCCGACTTAAATTGCGGTTAAGTTTcctc                                                            >  2:503274/1‑138 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCATGAGAAGGGA‑GCGCAGAGAGGGGTTTTCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTAcc                                                       >  2:190704/1‑140 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                       <  1:361835/140‑1 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                       <  1:739596/140‑1 (MQ=255)
                                                                          ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc                                                       <  1:1049179/140‑1 (MQ=255)
                                                                           taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc                                                      <  1:503274/140‑1 (MQ=255)
                                                                                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCATTTTTTTTTTTTTTGCC‑TGAGAAGGGA‑GCGGAGAGGGGGGTTATATTTGCTGACTGAACTGGCGATTAAATATCTGCGGAccccccgct                                                >  2:618177/1‑134 (MQ=255)
                                                                                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGAAGGGA‑GCGCAGAGGGGGGTTATCTTGCCCGACTGGACTGGGGATTTAAATTTCGCGTACCCCGCGtt                                                >  2:291580/1‑140 (MQ=255)
                                                                                 tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGAAGGGA‑GCGCAGAGAGGGGTTATTTTTCCTGACTGAACTGGGGGTTAAGTTTCCGCGTGCCCcgcggt                                                >  2:1126628/1‑138 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGAAGTGA‑GCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAATTTTCGCGTACCCCGCGgtt                                               >  2:803856/1‑140 (MQ=255)
                                                                                  cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGTGA‑GCGGAGGGAGGGATTATTTTTTCTGACTGAACTGGCGATTAAGTTTCCGCGGACCCGGCGgtt                                               >  2:82934/1‑140 (MQ=255)
                                                                                    ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                             <  1:30318/140‑1 (MQ=255)
                                                                                     gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGTGA‑GCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGGTTAAGGTTTCGCGGACCCCGCGGTTTGt                                            >  2:634548/1‑140 (MQ=255)
                                                                                        tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAt                                         <  2:1051497/140‑1 (MQ=255)
                                                                                          aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                       <  1:116761/140‑1 (MQ=255)
                                                                                            ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                     <  2:1097545/140‑1 (MQ=255)
                                                                                            ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                     <  2:985378/140‑1 (MQ=255)
                                                                                              cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGga                                   <  1:433633/140‑1 (MQ=255)
                                                                                                 gATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                <  2:831293/140‑1 (MQ=255)
                                                                                                    gCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCCTGAGAAGTGA‑GCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCGGTTTGTTATTGCGgatgatga                             >  2:99396/1‑140 (MQ=255)
                                                                                                    gCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTGCCATGAGAAGGGA‑GCGCAGGGAGGGATTATTTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTTGTTATTTCGggtgatga                             >  2:1057541/1‑140 (MQ=255)
                                                                                                        tGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATAtt                         <  1:16710/140‑1 (MQ=255)
                                                                                                             cgGATCTCGTTTGCGGCGGT‑CATTTTTTTTTTTTTTTCCCTGAGAAGTGA‑GGGCAGAGAGGGGTTATCTTTCCTTACTTAACTGGCGGTTAAATTTTCGCGTACCCCGCGGTTTGTTATTTTGGGTGATGATATTTaaag                    >  2:888540/1‑138 (MQ=255)
                                                                                                             cgGATCTCGTTTGCGGCGGG‑CATTTTTTTTTTTTTTTGCATTTGGAGAGAGGGGCGGGGAGGGA‑TATCTTTCCTTGCTGAAGTGGGGGCTAAGATTACCCCTACCCCGCGTTTGGTTGTTTCGGGTTATGATTTTTaaag                    >  2:535057/1‑138 (MQ=255)
                                                                                                                      tttGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGTGGATGATGATATTGAAGCTCGGCTGAt           <  2:1058281/140‑1 (MQ=255)
                                                                                                                               gT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCTGTGTTTAGTTATTGCGGATGATGATATTGAAGCACGGCTGATTGGTTCAtt  <  2:1190517/140‑1 (MQ=255)
                                                                                                                               gT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt  <  2:665477/140‑1 (MQ=255)
                                                                                                                               gT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt  <  1:1097765/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
ATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGT‑CATTTTTTCTTTTTCTGCCATGAGAAGTGA‑GCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATT  >  NC_000913/4604207‑4604473

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