Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,857 C→T 10.1% intergenic (‑99/+16) yecJ ← / ← azuC DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,8570CT10.1% 158.7 / 5.2 58intergenic (‑99/+16)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base T (4/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

TTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCG  >  NC_000913/1987558‑1988151
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                       
ccgatctaTTTCCTGCATGGTGCCCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCg                                                                                                                                                                                                                                                                                                        <  2:15884/294‑1 (MQ=255)
   ttGGATTGCCCGTATGGTGACCGGGTCTGCGGCGGCTGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAg                                                                                                                                                                                                                                                                                                     <  2:281971/301‑1 (MQ=255)
               gtgggTGGCCGGTTCGGGGGGGGCGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTGGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAACCTGCTCCGGCAGGtttttttgtg                                                                                                                                                                                                                                                                                      <  1:300243/298‑5 (MQ=255)
                  ggTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGATTTTTGTTaa                                                                                                                                                                                                                                                                                     >  2:14744/1‑300 (MQ=255)
                    tggccGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCTATATTCTTCAGCTTACTGGACATTTTTTCACGAACTTCAACGGGCGCGAGTACATCGAGAACATCCAGAATTTGTTTGCTAACCAACTGGCCAGCAACCACATCATAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAACAAACCTGCCGGAGCATGTTTTTTGTTATCg                                                                                                                                                                                                                                                                                   <  1:141983/298‑1 (MQ=255)
                     gACCGGGGCTGCGGCTGCGGGATGGTTATTTAAGTCGATTTTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCgg                                                                                                                                                                                                                                                                                  <  1:14744/301‑1 (MQ=255)
                            tCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATAtt                                                                                                                                                                                                                                                                          <  2:278579/301‑1 (MQ=255)
                                            ggTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACACAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATAGCCTGGCGGTACGtt                                                                                                                                                                                                                                                           >  1:161019/1‑300 (MQ=255)
                                             gTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAACCTGCTCCGGCAGGATTTTTTTTGTccagatcggaagagcaccatctaaactc                                                                                                                                                                                                                                                         >  2:242573/1‑274 (MQ=255)
                                                     aGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAACAAACCTGCCGGAACAGGATTTTTGTTATCGGAACCTACTGCCTGGCGGTACGGCTGTGAACGt                                                                                                                                                                                                                                                  >  1:341765/1‑301 (MQ=255)
                                                     aGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGAGGTACATCTTTGAACga                                                                                                                                                                                                                                                  >  2:40632/1‑300 (MQ=255)
                                                                        aatgggcgttttttCACTAGCTTCAGCGGGAGGGGTGACTGCGCGAACATCTAGAAGTTGTTGTATAACCAACGGGCAAGCAACCACATCAGGAACCAATGCCGGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATATTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGGGCAGGTTTTTGGTTATCGGAACATATTGCCTTGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAAc                                                                                                                                                                                                                               <  1:595307/293‑1 (MQ=255)
                                                                            ggagttcagacgtgtgctcttccgatctcgATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATAc                                                                                                                                                                                                                           <  2:541247/273‑1 (MQ=255)
                                                                             cgacgctcttccGATCTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACt                                                                                                                                                                                                                          <  1:122849/290‑1 (MQ=255)
                                                                              cATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTATTGCAATAGTTATAGAACATACtt                                                                                                                                                                                                                         >  1:143943/1‑301 (MQ=255)
                                                                                  ttttCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTATCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTAAACATACTTTTCa                                                                                                                                                                                                                     >  2:449567/1‑301 (MQ=255)
                                                                                     tCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATAATGCCTGGCGGTACGTCTTTAAACGTATTACAATAGTAATTGAACAAACTTTTCAggg                                                                                                                                                                                                                  >  1:451865/1‑300 (MQ=255)
                                                                                        cGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAAACTGCCGGAGCAGGTTTTTTGTTATCGGGACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTAAGGAttt                                                                                                                                                                                                               >  2:330328/1‑301 (MQ=255)
                                                                                         gAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAACAGGTATTTTGTTAACGGAACATCTTGCCTGCCTGTACGTCTTTGAACATCTTGCAATAGTTATTGAAACTACTTTTCAGTAtttt                                                                                                                                                                                                              >  2:554417/1‑301 (MQ=255)
                                                                                           aCTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTgc                                                                                                                                                                                                            <  2:229775/301‑1 (MQ=255)
                                                                                              tCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACtagatcggaagagcacaa                                                                                                                                                                                                         >  2:122973/1‑284 (MQ=255)
                                                                                               cAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCCACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCCTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAACCCTGCCGGACCAGGTTTTTTGTTATCGGAACATATTGCCTGACGATACGTATTTGAACGTCTAGCAATAGTTATTGAACATACTTTCCAGGAttataaccaa                                                                                                                                                                                                        >  1:204944/1‑293 (MQ=255)
                                                                                                    ggCGCGATTACTGTGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCa                                                                                                                                                                                                   <  2:452429/301‑1 (MQ=255)
                                                                                                        cgATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATAcagatcggaagagcgtcgtgtagggaacg                                                                                                                                                                                               >  1:540574/1‑273 (MQ=255)
                                                                                                          tttACATCGAGAACATCCATAATTTGTTTGATCACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGTTCGGCATTCAGCGGCTGGGACATCGTAAACTCCTTATAGCATTTGGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTGTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCgg                                                                                                                                                                                             <  2:342228/300‑1 (MQ=255)
                                                                                                              cATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGACACt                                                                                                                                                                                         >  2:165424/1‑301 (MQ=255)
                                                                                                              cATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCACGGCAca                                                                                                                                                                                         >  2:336965/1‑300 (MQ=255)
                                                                                                                   agaACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                                                                                                                                                                    >  2:169062/1‑301 (MQ=255)
                                                                                                                         tCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTg                                                                                                                                                                              <  2:144069/301‑1 (MQ=255)
                                                                                                                          ccAGTATTTTTTTGATAACCAACTGGAAATCATCCACATCAGATACCAATTCCTGATCGGCATTCAGCGGCTGGGACTTCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGt                                                                                                                                                                             <  1:165272/301‑1 (MQ=255)
                                                                                                                               ttttgtttgATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATAc                                                                                                                                                                        <  2:301803/300‑1 (MQ=255)
                                                                                                                                    ttttATAACCAACTGGCAAGCAACCACATCAGGAACCATTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtg                                                                                                                                                                   <  1:168865/301‑1 (MQ=255)
                                                                                                                                       gATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGAACTCCGACCATTTGTTATACAGGTGATa                                                                                                                                                                >  2:465907/1‑301 (MQ=255)
                                                                                                                                            ccAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTc                                                                                                                                                           >  2:118256/1‑301 (MQ=255)
                                                                                                                                                ccggCAAGTAACCACATCAGAAACCATTGCCTTATCGGCTTTCGTCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGGTATAGTAGCGACTTCACATCTGCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTGTTGTTATCGGGACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTg                                                                                                                                                       <  1:329924/299‑1 (MQ=255)
                                                                                                                                                         aaCCACATCGTAAACCAATTCCGTATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGTATTTTGAAAGCCGTTATAGTAGCGACTTCACAGCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCTGGTTTTTTGTTATCGGAACATACTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTataa                                                                                                                                              <  1:465343/301‑1 (MQ=255)
                                                                                                                                                           ccACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTataata                                                                                                                                            >  2:47372/1‑301 (MQ=255)
                                                                                                                                                            cacaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTATTATTATCTTTGCGCTTATAatat                                                                                                                                           >  1:123015/1‑301 (MQ=255)
                                                                                                                                                               aTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGCTTTTGCGCAGTTTCACCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCCATTATAATATTAc                                                                                                                                        >  1:179408/1‑301 (MQ=255)
                                                                                                                                                                   ttaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAAcc                                                                                                                                    <  2:139462/299‑1 (MQ=255)
                                                                                                                                                                     aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCaa                                                                                                                                  >  1:139337/1‑300 (MQ=255)
                                                                                                                                                                           ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACaaa                                                                                                                            >  2:301616/1‑301 (MQ=255)
                                                                                                                                                                                    ggCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCCTCGCTGCACTCCGACCATTTGTTATACAGGTTTTAATGTCTTTGCGCTTATAATATGACAACCATCACAAAaatcaatca                                                                                                                   >  2:266468/1‑300 (MQ=255)
                                                                                                                                                                                           aGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTAtgtg                                                                                                            <  2:161170/301‑1 (MQ=255)
                                                                                                                                                                                           aGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTAtgta                                                                                                            >  2:16205/1‑300 (MQ=255)
                                                                                                                                                                                             cggcgtggATATCGTAATATCCTGGTAGTATTTTGAAAGGTGTTTTAGTAGCGACGTCACATCTTCAGCGATAGTCTCATCCACCGTCATCAGTACACAAATAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGAt                                                                                                          <  1:336560/295‑1 (MQ=255)
                                                                                                                                                                                                       aTCTTAAACTCCTGATAGCATTTTGAAATCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATcaca                                                                                                <  2:404093/301‑1 (MQ=255)
                                                                                                                                                                                                        tCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACAt                                                                                               <  1:301212/301‑1 (MQ=255)
                                                                                                                                                                                                                                           agtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAACCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTAAAAAAATGATAATCaa                                                            >  1:372842/1‑301 (MQ=255)
                                                                                                                                                                                                                                                ctaCTTCACATCTTCAGCTATTGTCACATCTACCGTCATCAGGACACAAAACAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAAc                                                       <  1:47319/299‑1 (MQ=255)
                                                                                                                                                                                                                                                        caTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCCAAAAAATCAATCTTTATGTAATACAAATCAACAAAATAACCCTTTAATGTTATAAAAATGATAATCAAAAAACACcccccc                                               >  1:42600/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   ttAGTCACATCCACCTTGTGCAGGACAAAAATAATCCTGCTGGAGTAGGTTTTTTGTTATCGGAACATATTGCCTGTGGGTACGTCTTTGAACGTCGTGAAATAGTGAGTGAACATACTTGTTAGGATTTTGCGCAGTTTCAGTGCGGCACTCAGACCCTGTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGAAAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGaca                                    <  1:449003/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                         acttgccccGTCATCAGGACACAAAAAAACCTGCCGGTGCAGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACa                              <  1:118132/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       aaGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCATTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATa                <  2:179605/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              acaaaaaaCCTGCTCCGGCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATCTGACAACCATCACAAAAATCAATCTTTATGTGATACCAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAACAAGCCCCCCTATTTCTGACACCTACAGATGGACAAAAATAgcccat          >  1:583811/3‑297 (MQ=255)
                                                                                                                                                                                                                                                                                               caaaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGc        >  2:603205/2‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                 aaaaaCCTGCTCCGGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATCATATGACAACCATCACACAAATCAATCTTTATGTGATCCAAATCACATAAATACCCCTTTACTGTTATAACAATGATAATCAAAAAACAACCCCCATATTTCTGACACCTACAGAAGGCAAGAAAAAACGccacaca      >  1:246716/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                                     aCCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCg  >  2:475428/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                       
TTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGGTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCG  >  NC_000913/1987558‑1988151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: