Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,800,435 G→A G97D (GGT→GAT)  ygaM → DUF883 domain‑containing protein YgaM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,4350GA96.4% 164.9 / ‑6.7 55G97D (GGT→GAT) ygaMDUF883 domain‑containing protein YgaM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/28);  minor base T (2/0);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGATATCCAGAATGATGTCAATCAATTAGCTGACAGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGAC  >  NC_000913/2800188‑2800677
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cAGGATATCCAGAATGATGCCAATCAATTAGCTGACAGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGatat                                                                                                                                                                                                                                                 <  2:746851/251‑1 (MQ=255)
             atgatgTCAATCAATTAGCTGACAGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCa                                                                                                                                                                                                                                    <  1:139662/251‑1 (MQ=255)
               gatgTCAATCAATTAGCTGACAGCCTCGAATCTGTATTGAAGTCCTGGGTCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCAc                                                                                                                                                                                                                                   <  2:525028/250‑1 (MQ=255)
                 tgTCAATCAATTAGCTGACAGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACACTGGTTGAGCGGGGGTCCAGCAGCTGCGGTAGTTATTTTTATCGGCGCACTgg                                                                                                                                                                                                                                >  2:30008/1‑250 (MQ=255)
                    caatcaatTAGCTGACAGCCTAGAATCTGTATTGAAGTCCTGGGGCAGTGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGATGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGATGAGCGGGTTCTTTTGTTTGTGAAAGCCCCTGTGGTAGCGTGTGAACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTa                                                                                                                                                                                                                              >  2:421596/1‑250 (MQ=255)
                                    aGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGt                                                                                                                                                                                                              >  1:404146/1‑250 (MQ=255)
                                    aGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTa                                                                                                                                                                                                             <  1:77493/251‑1 (MQ=255)
                                      ccTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGt                                                                                                                                                                                                           <  1:716193/251‑1 (MQ=255)
                                             tCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGTAGCCGCACGCACTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGAGCGGATTCTTTTGTTCGTGCAAGACCCTGGTGTAGCGTGGGTACCGCAGCTGCGGTAGATCTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGtgaggcaaaa                                                                                                                                                                                                     >  2:2800/1‑241 (MQ=255)
                                             tCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCaaaa                                                                                                                                                                                                     >  2:110180/1‑250 (MQ=255)
                                                gTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGCAGCAGCACGACGTACAGCTCAGGCTTTGCTGAAAGAAGACAGGGCGCGAATGCCTGGCCGTACTCGCGTACAGCAAGACGCGCGAGATGCCGTTGACTACGCGGGTTCTTTTGCTCGTGAAAGACCCTGGGGTAGCGTGGGTGCAGCCGCTGCGGGAGACATTTTTGTCGGCGCACTGTTAAGCATgccaaactcggggtgagaaaatgg                                                                                                                                                                                                 >  2:247132/1‑229 (MQ=255)
                                                gTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGATGCGGTGGTTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGCAGTGccataatgt                                                                                                                                                                                                 >  2:701315/1‑243 (MQ=255)
                                                    tGAAGTCCTGGGGCAGCGACGCCACAGGGGAAGCTGAAGCAGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGAGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCATGGTGTAGCGTGGGTACAACAGCTGCGGGAGATAGTTTTATGGGCGAACTGTTAAGCATGCTCAAATAGTAGTGCAAAAATGataa                                                                                                                                                                                              >  2:657439/1‑250 (MQ=255)
                                                     gAAGTCCTGGAGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCTTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTACGCATGCGCAAATCGTAGTGCAAAAATGataat                                                                                                                                                                                             >  2:815992/1‑250 (MQ=255)
                                                              gggCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGAGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAACTACGCGt                                                                                                                                                                                   >  2:859818/1‑251 (MQ=255)
                                                              gggCAACGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGt                                                                                                                                                                                   <  1:177676/251‑1 (MQ=255)
                                                                gCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCt                                                                                                                                                                                 >  1:267784/1‑251 (MQ=255)
                                                                  aggggcGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCAGAAAGAAACCAGGGCGCGAATGCAGGGTCGTACGCGCGTCCAGCAAGCCGCGCGCGATTCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTTCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCttt                                                                                                                                                                               <  2:29974/246‑1 (MQ=255)
                                                                     gACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGCGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGAc                                                                                                                                                                            >  1:377439/1‑251 (MQ=255)
                                                                            aaGGGGAAGCTGAAGCCGCACGCAGTAAAGATCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAACGACCCTGGAGTAGCGGGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCCCTGTTAAGCATGCGCAAATCATAGTGCAAAAATGATAATAAATACGAGTCTTTGAcccaga                                                                                                                                                                       >  2:134496/1‑246 (MQ=255)
                                                                               gggAAGCTGAAGCCGCACGCAGGAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                                                                                                  <  1:2800/251‑1 (MQ=255)
                                                                                       gAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                                                                                            <  1:59030/249‑1 (MQ=255)
                                                                                           ccGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGGGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                                                                                      <  2:297915/251‑1 (MQ=255)
                                                                                              cacgcaGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCt                                                                                                                                                   <  2:802728/251‑1 (MQ=255)
                                                                                                  caGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                                                                               >  1:502587/1‑251 (MQ=255)
                                                                                                          cTCAGGCATTTCTGAAAGAAACCAGGGCGCGCATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                       >  1:247311/1‑251 (MQ=255)
                                                                                                            cAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                       <  1:77191/249‑1 (MQ=255)
                                                                                                            cAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCAGGGGCGTACTCGCGGCCAGCAAGCCGCGCGCGATGCCGTTGGCGGCGCGGAGTCTGTTGTTCGTGAAAGACCCTGGTGTAGCGGGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACGGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAAGAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                       <  1:112840/249‑1 (MQ=255)
                                                                                                                           gaaaCCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGgc                                                                                                                                                                                                                                       <  1:415445/138‑1 (MQ=255)
                                                                                                                              aCCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGCCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                                   <  1:776109/251‑1 (MQ=255)
                                                                                                                                 aGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                                                                                >  2:692930/1‑251 (MQ=255)
                                                                                                                                  gggCTCGAATGCAGGGGCGTACGCTTGTCCAGCAAGCCGCGCGCGAGGGCGGTGGCTGCGCGGATGTTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTCCAGCATCTGCGGTAGATATTTTTATCGGCGCATTGTTAAGCAGGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGACGCATGTCTTCGGGGTTTCGTTTTGCCTTGTGAATCACAAAAATCCCCCTAcccct                                                                                                               <  2:662814/251‑2 (MQ=255)
                                                                                                                                          aTGCATGGTCGTACTCGAGTCCAGCAAGCCGCGCGCGAGGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                                                                                       <  1:859818/251‑1 (MQ=255)
                                                                                                                                                agTCGTACTAGCGGCCAGAAAGCCGCGCGAGATGCAGTTGGGTGCGCGGCTTCTTTTGTTCGTGGAAGACCCTGGTGCAGCGTGGGTACAGCAGCTGCGGAAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATACATACGCGTCTTTGACACCGAAGCCTGACTTCGGGGTTTCTTTTTGCCTGGCGAATCACAAAAATCCCCCACCCCCGTCACGATCatattc                                                                                                 >  2:499478/2‑249 (MQ=255)
                                                                                                                                                   cGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAgg                                                                                              <  1:692930/251‑1 (MQ=255)
                                                                                                                                                          gcgTACAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAAttt                                                                                       <  1:565632/251‑1 (MQ=255)
                                                                                                                                                          ccggccAGCAAGCCGCGCGGGATGCCGTTTGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGCCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCGGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAAttt                                                                                       <  1:622804/247‑1 (MQ=255)
                                                                                                                                                              ccagcacgactcgcgcgATACCGTCGACTGCGCGGAATCTTTTGCTAGTGAAAGACCATGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATATTTATCGAAGACCTGTTAAGAATGCGCAAATCGTAGGGCAAAAATGATAATAAATACGCGTCTTTGACCCCGCAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGCATCACGAAAATCCCCCTACCCCGTCACGCTAATATCCCGGGTATTTTcacc                                                                                   >  2:397610/12‑248 (MQ=255)
                                                                                                                                                              ccAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                                   >  1:366948/1‑251 (MQ=255)
                                                                                                                                                                                          tGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTg                                                        >  2:161310/1‑250 (MQ=255)
                                                                                                                                                                                            cgcgGATGCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                                     <  1:805430/251‑1 (MQ=255)
                                                                                                                                                                                                 aTTCTTTGGTTCGTGAAAGACCCTGGTGGAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCGATAGATTGTGTGGTTGAATCtttt                                                <  2:267677/251‑1 (MQ=255)
                                                                                                                                                                                                       ttGTTCGTGAAAGACCCTGGTGTAGCGGGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                                           <  1:596603/250‑1 (MQ=255)
                                                                                                                                                                                                        tGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACt                                         >  1:598163/1‑251 (MQ=255)
                                                                                                                                                                                                             gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGAAATTTCGACCACTATTTGCTATATCTTGTGTGGTTGAATGTTTTTTCAACTACATc                                    >  1:654000/1‑251 (MQ=255)
                                                                                                                                                                                                              tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCt                                   >  1:308958/1‑251 (MQ=255)
                                                                                                                                                                                                                aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                                 <  1:293531/251‑1 (MQ=255)
                                                                                                                                                                                                                   gACCCTGGTTTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCACAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAt                              <  2:3840/251‑1 (MQ=255)
                                                                                                                                                                                                                   gACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAt                              >  1:50707/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTGGGTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGCTAATACATACGCTTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTATTTTTGCCTGGTGTATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCGGGGTAATTTCGACGACTATTTGCTATATATTGTGTAGTTGCATCTTTTTTCAAATCCATCTAGTATCTCTGTATCAACa                 >  2:682464/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gggTACAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACagag              <  2:303984/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aCAGCAGCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCTGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAgagaga           <  2:759512/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cagcagCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                                                                                            <  2:791834/104‑1 (MQ=255)
                                                                                                                                                                                                                                        cagcagCTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                                                                                            >  1:791834/1‑104 (MQ=255)
                                                                                                                                                                                                                                         agcagcTGCGGTAGATATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACa         >  1:856052/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ggCGGTAGATATTTTTATCGGCTCATAGTTAAGCATGCGAAAATCGTAGTGCAAAAATGATAATAAATACGCGGCTTTGACCCCGAAGCCTGTCTACGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCTCTATTTGCTATATATGGTGTGGTTGAATCTTGTTTCAACTACATCTAGTATCTCTGTATCAACAGAGATACAACCCGAc  <  2:751380/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CAGGATATCCAGAATGATGTCAATCAATTAGCTGACAGCCTCGAATCTGTATTGAAGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGAC  >  NC_000913/2800188‑2800677

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