Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,326,4500TG43.8% 97.6 / 18.4 65G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/33);  new base G (28/0);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.02e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326305‑2326598
                                                                                                                                                 |                                                                                                                                                    
cgcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                                  >  1:690473/1‑150 (MQ=255)
 gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                                <  2:414595/151‑1 (MQ=255)
   gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaaaa                                                                                                                                              <  2:186436/151‑1 (MQ=255)
   gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                                  >  1:849079/1‑147 (MQ=255)
     cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                                             <  1:700961/150‑1 (MQ=255)
      aGGCACTGTTAAAAATCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAata                                                                                                                                            >  2:114401/1‑150 (MQ=255)
      aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAatat                                                                                                                                           <  2:978688/151‑1 (MQ=255)
        gCACTGTTAAAAAGCAGTTTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTGCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCGGAGCATTGTGGCGGGGGGTATGGAAAATATGa                                                                                                                                         <  1:708297/151‑1 (MQ=255)
           cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGATCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGtaaagagtgg                                                                                                                                        >  1:101302/1‑143 (MQ=255)
             gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                                      <  1:286842/149‑1 (MQ=255)
             gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTa                                                                                                                                    <  1:691393/151‑1 (MQ=255)
               tAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGttt                                                                                                                                     >  2:516606/1‑148 (MQ=255)
                aagaaGCGGGCTGGCAGAAACGGAGTGCGGATTCACGGTCAATAAAGTATGCGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGGCAGGCGCAGAGCATTGTGGCGGGGGGGGGGGAAAAAATGAGtg                                                                                                                                      >  1:867614/4‑145 (MQ=255)
                 aaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGccc                                                                                                                                <  1:275106/151‑1 (MQ=255)
                 aaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATGAGTTTAGccc                                                                                                                                >  2:769530/1‑151 (MQ=255)
                 aaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGccc                                                                                                                                >  2:1138120/1‑151 (MQ=255)
                      cgggcgggCAGAAACGGTGTGCGGATTCACGGTCAATAACGGATGTGGTTCGGGTCGTAAAAGTGTGGCGCTTGCCGGCCAGGCCATTCAGGCAGGTCAGGGGCAGAGCATTGTGGCGGAGGGGACTGGAAATATGAgggtgacaaccttc                                                                                                                           >  1:153329/1‑138 (MQ=255)
                       gggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtact                                                                                                                          <  2:103425/151‑1 (MQ=255)
                       gggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtact                                                                                                                          <  1:282707/151‑1 (MQ=255)
                       gggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTTGAAAATATGGGTTTATtccccccc                                                                                                                           >  2:687950/1‑147 (MQ=255)
                         gCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTTGCCCCCtactt                                                                                                                         >  2:126036/1‑150 (MQ=255)
                            ggCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtacttact                                                                                                                      <  2:543000/150‑1 (MQ=255)
                             gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCtactt                                                                                                                         >  1:959812/1‑146 (MQ=255)
                              cAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTc                                                                                                                     >  2:862630/1‑149 (MQ=255)
                               aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGCGTGGAAAAAATTAGATTACTCCCCtacctttttcat                                                                                                                  >  2:318717/1‑142 (MQ=255)
                               aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGACCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGAt                                                                                                                  >  2:180349/1‑151 (MQ=255)
                                      gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAg                                                                                                           <  2:35848/151‑1 (MQ=255)
                                         tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                                                                                         <  1:712836/150‑1 (MQ=255)
                                         tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTCTTGCCCCCCACTTATTCAGTGCTAAAg                                                                                                           >  1:622571/1‑148 (MQ=255)
                                           cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGGGTTTTGCCCCCCACTTACTCGATGCAAAAGCACGc                                                                                                      >  1:425735/1‑151 (MQ=255)
                                                    ggTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTAt                                                                                             <  1:1064795/151‑1 (MQ=255)
                                                         aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGGAATATGAGTTTATCCCCCTACTTACTCGCTGCAAAAGCACGCTCCGGGTATCGt                                                                                          >  2:1048282/1‑149 (MQ=255)
                                                           aaaGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTg                                                                                      <  2:182289/151‑1 (MQ=255)
                                                           aaaGTATGTGGTTCGAGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTg                                                                                      <  1:957138/151‑1 (MQ=255)
                                                              gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCtttg                                                                                     >  2:704716/1‑147 (MQ=255)
                                                                 tgtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGgag                                                                                   >  2:1151973/1‑148 (MQ=255)
                                                                  gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGTTTAGCCCCCTACCTACTCGATGCAAAAGCACCCTCTGGGTTTCGTCTTTGGGAc                                                                                 >  1:362355/1‑149 (MQ=255)
                                                                   tgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGAAAAAATAGGTTTTGCCCCATTCATCTTTCATGCAAAAGCACACCCTTGTTTATGTCTTGGAGAc                                                                                 >  2:798586/1‑148 (MQ=255)
                                                                     gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACa                                                                            <  2:670794/151‑1 (MQ=255)
                                                                      ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGAAgg                                                                               >  2:390063/1‑147 (MQ=255)
                                                                        cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACAAGCTCTGGTTATCGTCTTGGGGACGGAAAgg                                                                          >  2:636659/1‑150 (MQ=255)
                                                                          ggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGttt                                                                       <  1:205278/151‑1 (MQ=255)
                                                                                aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAATATTAGGTTTGCCCCCCCCTTAATCGCTGCAAAAACACGCCCTGGGTATTTTCTTTGGGAGGGGCAGGgtt                                                                       >  1:1024078/1‑145 (MQ=255)
                                                                                 aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGTTTTGTCCCCCACTTTATTCATGAAAAAACAAGCCCTGGGTATCGGCTTTGGGACGGGCAGGTTTAAGa                                                                   >  1:269225/1‑148 (MQ=255)
                                                                                    gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAt                                                             <  1:1006955/151‑1 (MQ=255)
                                                                                    gtgtGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAAt                                                             <  2:577444/151‑1 (MQ=255)
                                                                                         gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGCAATATGAGTTTAGCCCCCTACTTAATCGATGCAAAAGCACGCTCTGGTTTTCGTCTTGGTGAACGGCAGGTTTATGACGTAATcc                                                           >  2:37122/1‑148 (MQ=255)
                                                                                           gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                      <  2:410778/151‑1 (MQ=255)
                                                                                           gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                      <  1:574428/151‑1 (MQ=255)
                                                                                                              caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGGCGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGcc                                    <  1:269433/150‑1 (MQ=255)
                                                                                                              caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGcc                                    <  2:342437/150‑1 (MQ=255)
                                                                                                               aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCAc                                  <  2:742982/151‑1 (MQ=255)
                                                                                                                  caggTCAGGCGCAGAGCATTGTGGCGGGGGGTGTGGTAAATATGTGGTTAGCCCCCTACTTAATCGCTGCAAAAGCACGCCCTGGGTTTTGGCTTTGGGGCGGACAGGGTTATGACGGAATCCTGCGCGATGGGCTGATGTGCGCCAccc                                >  1:97045/1‑150 (MQ=255)
                                                                                                                            gcAGAGCATTGTGGCGGGGGGGGAGGAAAAAATTAGTTTAGCCCCCCCCTTACTCGATGCCAAAGGACGCTCTGGTTTTCGTCTTGGGGACGGGGAGGGTTATGAAGTAAACCTCCGCGGTGGGCCGAAGTGCGCCCCCCCTgggt                          >  1:1134223/1‑144 (MQ=255)
                                                                                                                             cAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGGCAGGGTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCat                     >  2:893884/1‑150 (MQ=255)
                                                                                                                              agagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCata                    <  2:255958/150‑1 (MQ=255)
                                                                                                                                  cATTATGGCGGGGGGTATGGAAAATATGAGTTTAGCTCCCTAGTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTAGAGACGGACAGGTTTATAACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCGATGGTTATCATATg                  <  1:1004833/148‑1 (MQ=255)
                                                                                                                                   aTTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGa              <  1:126032/151‑1 (MQ=255)
                                                                                                                                      gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTa           <  2:635513/151‑1 (MQ=255)
                                                                                                                                      gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTa           <  2:952032/151‑1 (MQ=255)
                                                                                                                                         gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccg        <  1:237416/151‑1 (MQ=255)
                                                                                                                                         gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccg        >  2:269974/1‑151 (MQ=255)
                                                                                                                                           ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc      >  2:1238352/1‑151 (MQ=255)
                                                                                                                                            gggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCTTGATGTGCGCCACCCATGGTTATCATATGGGGATTCccgccg     <  1:159643/151‑1 (MQ=255)
                                                                                                                                               ggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaa  <  1:114400/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAA  >  NC_000913/2326305‑2326598

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