Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,727 G→A 5.2% P97S (CCG→TCG)  mgsA ← methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,7270GA 5.2% 252.3 / 2.6 77P97S (CCG→TCG) mgsAmethylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base A (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTA  >  NC_000913/1026583‑1026875
                                                                                                                                                |                                                                                                                                                    
gCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcggcacg                                                                                                                                                  >  6:273228/1‑149 (MQ=255)
        tCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAg                                                                                                                                          <  4:328498/149‑1 (MQ=255)
          ggggATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGCAGTCGGCCGGTGCCCCGTGGGTGGCGACCGGAATGGTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGGCAGGATCGTGCGGCACGGCATTTAGTg                                                                                                                                        <  5:22478/149‑1 (MQ=255)
          ggggATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTg                                                                                                                                        <  8:432140/149‑1 (MQ=255)
              aTCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATc                                                                                                                                    <  3:89909/149‑1 (MQ=255)
               tCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATcc                                                                                                                                   >  1:23198/1‑149 (MQ=255)
               tCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAATATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATcc                                                                                                                                   <  6:153580/149‑1 (MQ=255)
                cAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATcc                                                                                                                                   >  3:249682/1‑148 (MQ=255)
                   aaTATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCag                                                                                                                                >  4:71492/1‑148 (MQ=255)
                          acagcgTCGTTGAAATTCTGCGACTGGATTATGAAGTCTTCCGTAGCCACTATGGTGGCGACCGGAATTAACCATACCGTAGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAAt                                                                                                                        >  2:182111/4‑149 (MQ=255)
                           ccGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATc                                                                                                                       <  8:49891/149‑1 (MQ=255)
                             gcgTCGTTGAAATGCGGCGACTGGATTATGCAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCaa                                                                                                                     <  3:235163/149‑1 (MQ=255)
                             gcgTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCaa                                                                                                                     <  5:423260/149‑1 (MQ=255)
                                   ttGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACAt                                                                                                                >  5:272883/1‑148 (MQ=255)
                                            ggcgACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTccc                                                                                                      <  1:342334/149‑1 (MQ=255)
                                            ggcgACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTccc                                                                                                      >  4:16050/1‑149 (MQ=255)
                                                aCTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCt                                                                                                  >  8:238541/1‑149 (MQ=255)
                                                aCTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCt                                                                                                  <  6:272883/149‑1 (MQ=255)
                                                  tGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCt                                                                                                  >  6:391354/1‑147 (MQ=255)
                                                  tGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCt                                                                                                  >  4:213143/1‑147 (MQ=255)
                                                    gATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTgaga                                                                                              <  2:315053/149‑1 (MQ=255)
                                                    gATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTga                                                                                                <  6:379062/147‑1 (MQ=255)
                                                    gATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTga                                                                                                <  2:292453/147‑1 (MQ=255)
                                                    gATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGACACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTga                                                                                                <  7:28023/147‑1 (MQ=255)
                                                        aTGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCaa                                                                                          <  3:71492/149‑1 (MQ=255)
                                                         tGAAGTCTGCCGTTGCCCCGTGGGGGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCaa                                                                                          <  4:249682/148‑1 (MQ=255)
                                                         tGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAAt                                                                                         >  8:50594/1‑149 (MQ=255)
                                                           aaGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTgg                                                                                                                                       >  7:232943/1‑101 (MQ=255)
                                                           aaGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTgg                                                                                                                                       <  8:232943/101‑1 (MQ=255)
                                                             gTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATgcgc                                                                                     <  5:273228/149‑1 (MQ=255)
                                                              tCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTg                                                                                                 <  4:402120/136‑1 (MQ=255)
                                                              tCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTg                                                                                                 >  3:402120/1‑136 (MQ=255)
                                                              tCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATgcgc                                                                                     >  3:348316/1‑148 (MQ=255)
                                                                     ttGCCACGTTGGTGGCGACCGGAATGTTCCAAACCGTCGCCAGACGCAGCAAGGCGTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACctgc                                                                             <  2:205105/149‑1 (MQ=255)
                                                                        ccACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGACACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTgg                                                                          >  1:211346/1‑149 (MQ=255)
                                                                            gTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAc                                                                      <  1:388919/149‑1 (MQ=255)
                                                                              tggtggCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAccc                                                                    >  3:273842/1‑149 (MQ=255)
                                                                                gtggCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAccccc                                                                  >  3:379103/1‑149 (MQ=255)
                                                                                     gACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTg                                                             >  8:85938/1‑149 (MQ=255)
                                                                                       ccGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTggg                                                           <  7:238541/149‑1 (MQ=255)
                                                                                          gAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTga                                                                                                >  7:403048/1‑109 (MQ=255)
                                                                                          gAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTga                                                                                                <  8:403048/109‑1 (MQ=255)
                                                                                          gAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCa                                                        >  4:190023/1‑149 (MQ=255)
                                                                                             tGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTc                                                     <  4:273842/149‑1 (MQ=255)
                                                                                               ttCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGCAGAAAATCACTACATCAATTTTCCCTTCTGAGATCAATGCGCCACCCTGCTGGTCACCCCCCATTGGGCCACTCaa                                                   >  3:18700/1‑149 (MQ=255)
                                                                                               ttCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCACCCTGCTGGTCACCCCCCATTGGGCCACTCaa                                                   >  7:383995/1‑149 (MQ=255)
                                                                                                   aTACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATc                                               <  8:144499/149‑1 (MQ=255)
                                                                                                     aCCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATcgc                                             >  6:204011/1‑149 (MQ=255)
                                                                                                       cGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACctgc                                                                             <  3:157987/115‑1 (MQ=255)
                                                                                                       cGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACctgc                                                                             >  4:157987/1‑115 (MQ=255)
                                                                                                       cGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACctgc                                                                             >  4:163920/1‑115 (MQ=255)
                                                                                                       cGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACctgc                                                                             <  3:163920/115‑1 (MQ=255)
                                                                                                       cGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGt                                           <  4:379103/149‑1 (MQ=255)
                                                                                                          cGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGa                                        <  3:16050/149‑1 (MQ=255)
                                                                                                              aGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGtt                                    >  3:182698/1‑149 (MQ=255)
                                                                                                              aGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGtt                                    >  1:6076/1‑149 (MQ=255)
                                                                                                               gACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTc                                   <  6:188061/149‑1 (MQ=255)
                                                                                                                aCGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCa                                  <  3:237375/149‑1 (MQ=255)
                                                                                                                aCGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCa                                  <  4:348316/149‑1 (MQ=255)
                                                                                                                aCGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCa                                  <  5:364943/149‑1 (MQ=255)
                                                                                                                 cGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCAGTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCAt                                 <  3:204114/149‑1 (MQ=255)
                                                                                                                  gcagcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATg                                >  7:418712/1‑149 (MQ=255)
                                                                                                                    agcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCAtt                                                              <  8:46361/117‑1 (MQ=255)
                                                                                                                    agcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCAtt                                                              >  7:46361/1‑117 (MQ=255)
                                                                                                                     gcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCg                             <  8:385221/149‑1 (MQ=255)
                                                                                                                      caAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCGCCCCCCATTGGGCCACTCAACCTCGCGAGGACGTTCATGCCgg                            >  8:158290/1‑149 (MQ=255)
                                                                                                                          gCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgc                        >  8:117830/1‑149 (MQ=255)
                                                                                                                                  gTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                >  3:158930/1‑149 (MQ=255)
                                                                                                                                  gTCAGGATCGTGCGACACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATc                                               >  5:84334/1‑118 (MQ=255)
                                                                                                                                  gTCAGGATCGTGCGACACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATc                                               <  6:84334/118‑1 (MQ=255)
                                                                                                                                   tCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAAt               >  6:401689/1‑149 (MQ=255)
                                                                                                                                     aGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAAtt              >  6:231489/1‑148 (MQ=255)
                                                                                                                                      ggATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTaa            <  2:321454/149‑1 (MQ=255)
                                                                                                                                      ggATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTaa            <  5:204011/149‑1 (MQ=255)
                                                                                                                                       gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAg           >  4:363486/1‑149 (MQ=255)
                                                                                                                                       gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAg           <  7:50594/149‑1 (MQ=255)
                                                                                                                                       gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAg           <  7:241093/149‑1 (MQ=255)
                                                                                                                                         tCGTGCGGCACGGCATTTAGTGGATCCCAGATGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGtt         <  5:231607/149‑1 (MQ=255)
                                                                                                                                              cggcacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCgg    <  2:23198/149‑1 (MQ=255)
                                                                                                                                               ggcacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCACTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCgg    >  6:348103/1‑148 (MQ=255)
                                                                                                                                                gcacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTa  <  4:158930/149‑1 (MQ=255)
                                                                                                                                                gcacggcaGTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTa  <  1:385030/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTA  >  NC_000913/1026583‑1026875

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