Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,0000TG51.2% 81.9 / 105.5 87Y557S (TAC→TCC) maeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base T (17/25);  new base G (44/0);  total (62/25)
Fisher's exact test for biased strand distribution p-value = 8.28e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-09

TCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTAC  >  NC_000913/1553851‑1554148
                                                                                                                                                     |                                                                                                                                                    
tCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTa                                                                                                                                                     <  2:714132/151‑1 (MQ=255)
  ggATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGAtt                                                                                                                                                   >  1:784245/1‑151 (MQ=255)
       atgTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCgg                                                                                                                                                >  1:974194/1‑149 (MQ=255)
       atgTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGc                                                                                                                                               <  2:51225/150‑1 (MQ=255)
       atgTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGc                                                                                                                                               <  1:327246/150‑1 (MQ=255)
       atgTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTggggt                                                                                                                                              >  2:151257/1‑149 (MQ=255)
         gTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTGCGGCGGTAGTCGCGGGATTCGGCTTg                                                                                                                                            >  2:1753983/1‑151 (MQ=255)
         gTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGGTTg                                                                                                                                            >  2:1006943/1‑151 (MQ=255)
              gCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTggccga                                                                                                                                       >  2:786561/1‑146 (MQ=255)
              gCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGGCAGa                                                                                                                                       >  2:1775433/1‑151 (MQ=255)
              gCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTTGGCTTGCCAGa                                                                                                                                       >  2:726558/1‑151 (MQ=255)
                 tAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTTCCAGAAAt                                                                                                                                    >  2:1676838/1‑151 (MQ=255)
                 tAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTTGCAGAAAt                                                                                                                                    >  1:823584/1‑151 (MQ=255)
                 tAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAAt                                                                                                                                    >  2:684690/1‑151 (MQ=255)
                 tAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAt                                                                                                                                    >  2:515797/1‑151 (MQ=255)
                 tAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGCCTGCCAGgaaa                                                                                                                                    >  1:1667028/1‑150 (MQ=255)
                   gCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGGTTGCCCGaaaatt                                                                                                                                  >  2:103950/1‑148 (MQ=255)
                    cAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTa                                                                                                                                  <  1:436839/150‑1 (MQ=255)
                             cAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATgg                                                                                                                        <  1:684387/151‑1 (MQ=255)
                               aaTGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGc                                                                                                                       >  2:1602503/1‑150 (MQ=255)
                                aTGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt                                                                                                                     <  1:1161501/151‑1 (MQ=255)
                                  gagaTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCCATGGCCCGt                                                                                                                   >  1:1760331/1‑151 (MQ=255)
                                    gaTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGttt                                                                                                                 >  2:1214355/1‑150 (MQ=255)
                                    gaTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCCTGCCCGAAATTTTCTTCAATTGCCTGTTg                                                                                                                 >  2:1085792/1‑151 (MQ=255)
                                    gaTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCCGTTg                                                                                                                 >  2:1558058/1‑151 (MQ=255)
                                     aTTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGc                                                                                                                <  1:151211/151‑1 (MQ=255)
                                       ttATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAg                                                                                                              <  2:372727/151‑1 (MQ=255)
                                       ttATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGGTTGCCAGAAAATATCGGCAAAGGCCTGTTGCAg                                                                                                              >  2:1502791/1‑151 (MQ=255)
                                         aTTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAggg                                                                                                            <  2:278822/151‑1 (MQ=255)
                                         aTTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAggg                                                                                                            <  2:204124/151‑1 (MQ=255)
                                         aTTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAggg                                                                                                            >  2:1320067/1‑151 (MQ=255)
                                         aTTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCCGAAAATATCGTCAAAGGCCTGTTGCAggg                                                                                                            >  2:947831/1‑151 (MQ=255)
                                          ttCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGc                                                                                                           <  1:947317/151‑1 (MQ=255)
                                              ccACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTTTCGTCAATGGGCTGTTGTAGGGCtttg                                                                                                       >  1:1656433/1‑149 (MQ=255)
                                              ccACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTGGCAGAAATTATCGTCCATGGCCTGTTGCAGGGCTTCg                                                                                                       >  2:277586/1‑151 (MQ=255)
                                                aCTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGc                                                                                                     <  2:1660162/151‑1 (MQ=255)
                                                aCTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCCGTTGCAGGGCTTCgg                                                                                                      >  1:640850/1‑150 (MQ=255)
                                                   aCCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGc                                                                                                     <  1:1302323/148‑1 (MQ=255)
                                                   aCCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTAATCGGCTTTCCAGAAATTATCCTCAATGGCCTGTTGCAGGGCTTCGGCaga                                                                                                  >  2:1142894/1‑151 (MQ=255)
                                                   aCCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGTTTTGGCTTGCCAGCAAATATTGTCAATGGCCTGTTGCAGGGCTTCGGCagg                                                                                                  >  2:1719162/1‑150 (MQ=255)
                                                   aCCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGCAGGGCTTCGGCaga                                                                                                  >  1:1098254/1‑151 (MQ=255)
                                                   aCCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGGTTCGGc                                                                                                     >  2:1302914/1‑148 (MQ=255)
                                                     ccAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAgg                                                                                                <  1:1671765/151‑1 (MQ=255)
                                                        gTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGGTTCGGCCGAGGgtt                                                                                             >  2:1298503/1‑151 (MQ=255)
                                                        gTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGgtt                                                                                             >  2:60752/1‑151 (MQ=255)
                                                        gTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGGTTGCCAGAAATTATCGTCAAAGGGCTGTTGCAGGGCTTCGGGAGAGGgtt                                                                                             >  1:693141/1‑151 (MQ=255)
                                                          aTGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCCGTTGCAGGGCTTCGGCAGAGGttttt                                                                                           >  1:701902/1‑150 (MQ=255)
                                                             gATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCAc                                                                                         <  1:1771000/150‑1 (MQ=255)
                                                             gATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCAc                                                                                         <  2:1295395/150‑1 (MQ=255)
                                                             gATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAATCGCGGGATTCGGCCTGCCAGAAAATATCGTCAATGGCCTGTTGCAGGGCTTCGGCCGAGGGTTTCAc                                                                                         >  2:1600553/1‑150 (MQ=255)
                                                             gATGAGATCTGAAAAACGAAGAGGGAAAAGGCCCGGTAGACTTCACTACCGGGCGCAGGCTTAGATGGAGGTCCGGGGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCAc                                                                                         <  1:1469182/150‑1 (MQ=255)
                                                                     cTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGCAGGGCGTCGGCAGAGGGTTTTACCGCCACGc                                                                                 >  1:386877/1‑150 (MQ=255)
                                                                       gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCCGAAATTATCGTCAATGGCCTGTTTCAGGGCTTCGGCAGGGGGTTTTACCCCCCCGcc                                                                                >  2:240866/1‑149 (MQ=255)
                                                                       gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTTCCAGAAATTTTCGTCAATGGCCTGTTTCAGGGCTTCGGCAGAGGCTTTAACCGCCCCGCCtt                                                                              >  1:1718906/1‑151 (MQ=255)
                                                                       gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGCAGGGGTTCGGCAGAGGTTTTTAccgccccgcc                                                                                >  2:47084/1‑149 (MQ=255)
                                                                       gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTCGGCTTGCCCGAAATTATTCTCAATGGCCTGTTGCAGGGCTTCGGGAGAGGGTTTCACCGccccccc                                                                                >  1:1187052/1‑144 (MQ=255)
                                                                              ggAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAGGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTtgctgc                                                                         <  2:205138/149‑1 (MQ=255)
                                                                                   gggAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTGCCAGCAAATATCGTCGATAGACCGATTCAGGGACTCGGCAAAGGGTTTTACCGCCCCGACTTGATGCTgcggc                                                                   >  1:1023138/1‑148 (MQ=255)
                                                                                        aTAGCCCGGTAGCCTTCACTACCAGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCCGAAAATATCGTCAATGGCCTGTTGCAGGGCTTTGGGAGAGGGTTTTACCGCCACGCCCTGCTGCTGCGCCCTTTTg                                                             >  1:19065/1‑151 (MQ=255)
                                                                                           gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAAAGGCCTGTTTCAGGGGTTCGGCAGAGGGTTTTACCGCCCCCCCCTTCTGCTGCGCCATTTTGCCa                                                          >  2:1751472/1‑151 (MQ=255)
                                                                                            ccgggTACCCTTCCCTACCGGGGGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCa                                                          <  2:380882/147‑1 (MQ=255)
                                                                                                gTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCCCGCCTTGCTGCTGCGCCCTTTTTCCAACCg                                                      >  2:44329/1‑150 (MQ=255)
                                                                                                     cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGGCCGAAATTATCGTCAATGGCCTGTTGCAGGGCCTCGGCAGAGGGTTTTACCGCCCCCCCTTTCTGCTGCGCCATTTTTCCAACCGCAAACg                                                >  1:1286337/1‑151 (MQ=255)
                                                                                                     cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCCGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCCCGCCTTGCTGCTGCGCCCTTTTTCCAACCGCCAAc                                                 >  2:1279515/1‑150 (MQ=255)
                                                                                                     cTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGACGCGGGCTTAGGCTTGCCCGGAATTTTTGTCCATGGGCTGTTGCAGGGCCTCGGCCGAGGGTTTTACCGCCCCGCCTTGCTTCTGCGCCATTTTTCCAACCGCAAACg                                                >  1:580013/1‑151 (MQ=255)
                                                                                                       tCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTCTGCGCGGAATTTGGGTTGCCAGAAATTTTCGACAAAGGCCTGTTGCAGGGCTTCGGCCGAGGGTTTCCCCGCCACGCCCTTCTGCCCCGCCATTTTGCCAaacgcaaacgca                                              >  2:1434513/1‑151 (MQ=255)
                                                                                                       tCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTCAATGGCCTGTTGCAGGGGGTCGGCAGAGGGTTTCACCGCCACGCCCTTCTGCTGCGCCATTTTTCCAACCGCCAACGCa                                              >  1:1285404/1‑151 (MQ=255)
                                                                                                       tCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCCTTTTGCCAACCGCAAACGc                                               >  2:1389415/1‑150 (MQ=255)
                                                                                                           tACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTg                                          <  1:98358/151‑1 (MQ=255)
                                                                                                              cGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGTTTCGGCAGAGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTgcgc                                       <  1:60732/151‑1 (MQ=255)
                                                                                                               gggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTTGCTGCTGCGCCCTTTTGCCAACCGCAAACGCCATTGcgcg                                      >  2:1240335/1‑151 (MQ=255)
                                                                                                                      ggCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTTGGCTTGCCAGAAATTTTCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCCCGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCCATTGCGCGGggga                                 >  2:1278021/1‑149 (MQ=255)
                                                                                                                      ggCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATTGTCAATGGCCTGTTGGAGGGGTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTTCCAACCGCAAAAGCAATTTCGCGGggga                                 >  2:1120040/1‑149 (MQ=255)
                                                                                                                      ggCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGCCGGGCTTCGGCAGAGGGTTTTACCGCCCCGCCTTTCTGCTGCGCCCTTTTTCCAACCCCAAACCCAATTGCGCGGggga                                 >  2:35296/1‑149 (MQ=255)
                                                                                                                      ggCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTTGTGTTTCCTGAAATTTTTGTCAAGGGGCTGTTGCATGGCTTCGGGCGAGGTTTTTACCGCCCCGCCTCTCCGCTTCGCCATTTTTCCCACCGCAAAAGCCATTTcgcgtggga                                 >  1:191421/1‑144 (MQ=255)
                                                                                                                      ggCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTTCTGCTGCGCCCTTTTGCCAACCGCAAACGCAATTGCGCGGggga                                 >  2:1685891/1‑149 (MQ=255)
                                                                                                                      ggCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCGGCTTTCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCCCCGCCCCGCCTTTCTGCTGCGCCACTTTGCCAACCGCCAACGCCATTGCGCGGggga                                 >  2:1169460/1‑149 (MQ=255)
                                                                                                                        cTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGGTTGCCAGAAATTATCGTCAATGGCCTGTTTCAGGGCTTCGGCAGAGGGTTTCACCCCCCCGCCTTGCCGCTGCGCCATTTTGCCAACCGCCAACGCAATTGCGCGGGAGACttt                             >  2:1246423/1‑151 (MQ=255)
                                                                                                                            gATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCCCGCCTTGCTGCTGCGCCCTTTTGCCAACCGCAAACGCAATTTCGCGGGGGACTTTTTGa                         >  2:435013/1‑151 (MQ=255)
                                                                                                                             aTGGAGGTACGGCGGTAGTCGTGGTTTTCGGCTCGCCAAAAATTGTCGTCAAGGGCTTGTTGCAGGGTTTCGGAAGAGTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGaa                        <  1:495986/151‑1 (MQ=255)
                                                                                                                             aTGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTg                          <  2:538534/149‑1 (MQ=255)
                                                                                                                              tGGAGGTAGGGGGGTAGGCGGGGTATCCGGCTTGCCAGAAATTATCGCCAATGGCCTGTTGCGGGGTTTCGGCAGAGGTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGaa                        <  1:1430781/150‑1 (MQ=255)
                                                                                                                                gAGGTACGGCGGTAGTCGCGGGATTTGGCTTGCCAGAAATTATTGTCAATGGCCCGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCCTTCTGCTGCGCCATTTTGCCCAACGCCAACGCAATTGCGCCGGGGACTTTTTGGatat                     >  1:511355/1‑151 (MQ=255)
                                                                                                                                gAGGTACGGCGGTAGGCGGGGTATTCGGCTTTCCAGAAATTATCGTCAATGGCCTGTTGCGGGGCTTCGGCAGAGGGTTTCACCGCCCCCCCTTGCTGCTGCGCCATTTTGACAACCGCCAACGCAATTTCGCGGGAGAGTTTCTTAat                       >  1:1073763/1‑149 (MQ=255)
                                                                                                                                    tgcggcggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCGGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCtt                  <  2:517051/148‑1 (MQ=255)
                                                                                                                                                    gTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTAc  >  1:960753/1‑150 (MQ=255)
                                                                                                                                                    gTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTAc  >  1:773758/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCGCCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGCAAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTAC  >  NC_000913/1553851‑1554148

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