Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,7350AT42.0% 56.5 / 37.4 50Y2* (TAT→TAA) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (4/25);  new base T (21/0);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 3.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413588‑1413885
                                                                                                                                                      |                                                                                                                                                           
cAGTTTTACTTTCTTCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAGTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTtata                                                                                                                                                           <  2:51810/150‑1 (MQ=255)
 aGTTTTACTTTCTGCGCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTtttat                                                                                                                                                          >  2:536762/1‑149 (MQ=255)
       aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAac                                                                                                                                                    <  1:1048150/150‑1 (MQ=255)
           tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑A‑AAA‑TTTATTCATTAACACCT‑c                                                                                                                                              >  1:266547/1‑151 (MQ=255)
              gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTT‑AAT‑TTTTTTCATGAA‑ACCT‑cc                                                                                                                                             >  2:645493/1‑149 (MQ=255)
              gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCa                                                                                                                                            <  2:649203/150‑1 (MQ=255)
                tctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑TAA‑TTTTTTCATTAGAACCT‑CCctt                                                                                                                                          >  2:687367/1‑150 (MQ=255)
                   ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCatta                                                                                                                                         <  1:794903/148‑1 (MQ=255)
                     ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCatta‑ttatt                                                                                                                                   <  1:513312/151‑1 (MQ=255)
                     ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCatta‑ttat                                                                                                                                    <  2:670075/150‑1 (MQ=255)
                            ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTg                                                                                                                           <  1:421651/151‑1 (MQ=255)
                             cTGCTTTTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑g                                                                                                                             <  2:874763/148‑1 (MQ=255)
                                 ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTACATGAACACCT‑CCCTtt‑tttttt                                                                                                                                  >  1:817548/1‑140 (MQ=255)
                                 ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAAAAAAAACCT‑TCACTA‑TTTTTTTCA‑ATgggt                                                                                                                        >  1:697499/1‑147 (MQ=255)
                                     tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAAT‑TTTTAAAATGAAAACCC‑CCCTTT‑TTTTTTCCG‑GGGGtgggtt                                                                                                                    >  1:904546/1‑145 (MQ=255)
                                     tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑ttt                                                                                                                  <  1:731829/150‑1 (MQ=255)
                                      ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAA‑TTTTTTTATGAAAAACC‑CCCCTA‑TTAttttttgtgggttgtgt                                                                                                                    >  2:1109947/1‑135 (MQ=255)
                                        ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTAtt                                                                                                               <  2:670095/150‑1 (MQ=255)
                                         tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTATATGAACACCT‑CCCTTA‑TTTTTTTCatcgggttgtttttt                                                                                                                 >  1:17719/1‑135 (MQ=255)
                                           ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTACATGTACACCC‑CCCTTA‑TTTTTTTCA‑ATGGGTTG‑tttttt                                                                                                               >  1:380391/1‑147 (MQ=255)
                                           ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCAt                                                                                                           <  1:1102343/151‑1 (MQ=255)
                                            cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTATATTAAAAACC‑CCCTTT‑TTTTTTTCA‑GTGGGTCG‑ttttttt                                                                                                              >  2:542734/1‑146 (MQ=255)
                                            cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAATTAAAACCT‑CCCTTT‑TTATTTTCC‑GTGGGTCG‑tttttttc                                                                                                             >  1:1037128/1‑146 (MQ=255)
                                                 agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTTAATAAACACCTACCTTTT‑TTTTTTCAT‑GTGGTTGT‑TTTTTTCcttctt                                                                                                        >  2:1073302/1‑149 (MQ=255)
                                                   cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAg                                                                                                   <  1:2306/151‑1 (MQ=255)
                                                      tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTg                                                                                                 <  1:444150/150‑1 (MQ=255)
                                                       ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTTAAAAAACACCC‑CCCTTT‑TTTTTTTCT‑GTGGTTTG‑TTTTTTTCCTTTTTTGAGt                                                                                                  >  1:1067449/1‑148 (MQ=255)
                                                           cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGcttct                                                                                            <  1:772348/150‑1 (MQ=255)
                                                              aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTT‑AAT‑TTTTTA‑ATTAACACCC‑CCCTTT‑TTTTTTTCT‑GTGGTTTG‑TTTTTTTCATTTTTTTTGTGCtttttttt                                                                                        >  2:249596/1‑151 (MQ=255)
                                                              aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑TAT‑TTTATTAATGAACACAT‑CCCTTATTTTTTTCCT‑GGGGTTGG‑TTTTTTctttcttttggaggtttttttt                                                                                         >  2:616191/1‑137 (MQ=255)
                                                                tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAA‑TTTTTATATGTAAAACC‑CCCTTA‑TTTTTTTTATGGGGTTGG‑TTATTTCCTTTTTTC‑AGTGCTTCttttt                                                                                        >  1:364178/1‑149 (MQ=255)
                                                                   gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AATATTTTTATAAGAAAACAC‑CATTTA‑TTTTTTCTA‑‑TGGTTTGTTTTATTCTA‑CTTTCTAGTGCCTttttttttttt                                                                                   >  2:1046190/1‑146 (MQ=255)
                                                                    aTTTGGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATAAATAAAAACAT‑CCATTT‑TTTTTTTCA‑GTGGGTGG‑GTTTTTTTATCTTTTGAGTGGTTGtttttttctt                                                                                   >  2:379832/1‑145 (MQ=255)
                                                                       tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑G‑AAT‑TTTTTAAAAGAAAACCC‑CCACTA‑TTTTTTTCC‑ATGGGTCG‑TTTTTTTCTTTTTTTTAGTGCGTttttttttctt                                                                                   >  2:753933/1‑142 (MQ=255)
                                                                          ttatttttcaTTTTGAATATCTTGTCCATTTAGTAGAAGTGCCACCTTCCTTTTCAATAGTGGCGC‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GGGGTTCT‑TTTATTCCATCTTTCGATTGCTTCTTTTTCATTTccac                                                                               <  2:880934/139‑1 (MQ=255)
                                                                               ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AAcc                                                                       <  1:1022052/150‑1 (MQ=255)
                                                                                 aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCgg                                                                     <  2:904546/150‑1 (MQ=255)
                                                                                 aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGt                                                                    <  1:458344/151‑1 (MQ=255)
                                                                                   aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑CAT‑TTTTTATAAATAAACCC‑CCACTT‑TTATTTTCC‑CTAGGTCG‑TTTTTTTCCTCTTTTGTCGGCTTTTTTTTTTCTTCCACCCT‑AAcccg                                                                     >  1:276474/1‑146 (MQ=255)
                                                                                    ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCCCTTCCACCAT‑AACCg                                                                      <  2:747797/146‑1 (MQ=255)
                                                                                      ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCgg                                                                <  1:981593/150‑1 (MQ=255)
                                                                                      ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCgg                                                                <  1:300861/150‑1 (MQ=255)
                                                                                                tGTCCAGTTAGTAGGAGGGCCACCTTCCTTTTCAATAGTGGCGG‑T‑GAT‑TTTTTTATTTAATAACT‑CCCTTC‑TTTTTTTTT‑GTTGTTCG‑GTTTTTTCTTTTTTCCCATGCTTCTTTTTTTCTTTCTCCCT‑AATCGGCTTggggtgaccct                                                       >  2:962340/1‑141 (MQ=255)
                                                                                                 gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATcc                                                     <  2:141344/150‑1 (MQ=255)
                                                                                                  tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTACATGAACAACC‑CCACTATTTTTTTCAA‑GTGG‑TCG‑TTTTTTTTATCTTTTTTGTTCTTTTTTTTTTCTTTCACCCC‑AACCCGGTCGGGGTTCTCATcc                                                     >  2:205122/1‑149 (MQ=255)
                                                                                                  tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAATTAACACCC‑CCCTTT‑TTATTTTCA‑GTGGGTCG‑TTTTTTTCCTCTTTCTTGGGCTTCTTTTTTTCTTTCCCCAT‑AACCAGTTTGGGGTTCGCATcc                                                     >  1:925962/1‑149 (MQ=255)
                                                                                                            agGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAAAGAACACCC‑CCCCTT‑TTATTTTCA‑GTGGGTTG‑TTTTTTTCTTTTTTTGAGGGCTTTTTTTTTTTTTCCCCCCTAAACCGGTTGGGGGTT‑GCATCCCTTACACtt                                            >  1:1170449/1‑148 (MQ=255)
                                                                                                              gAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCt                                        <  2:328860/150‑1 (MQ=255)
                                                                                                                gTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg                                     <  1:502659/151‑1 (MQ=255)
                                                                                                                  gCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTa                                   <  2:855017/151‑1 (MQ=255)
                                                                                                                             tttCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGa                          <  2:774537/149‑1 (MQ=255)
                                                                                                                                        gCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc             >  1:382050/1‑151 (MQ=255)
                                                                                                                                                      ataCATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  1:760778/151‑1 (MQ=255)
                                                                                                                                                      ataCATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   <  2:704889/150‑1 (MQ=255)
                                                                                                                                                      ataCATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   <  1:1209361/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                           
CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413588‑1413885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑