Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,842,9140AG41.4% 68.9 / 16.2 59M182V (ATG→GTG) gdhAglutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base A (3/31);  new base G (24/0);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 2.28e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTT  >  NC_000913/1842769‑1843058
                                                                                                                                                 |                                                                                                                                                      
cGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatgat                                                                                                                                                    <  2:369642/150‑1 (MQ=255)
     cGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGAAGgaaa                                                                                                                                               >  2:1313351/1‑150 (MQ=255)
          ggAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGctct                                                                                                                                         <  1:1349601/151‑1 (MQ=255)
            aaaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAGACCTcc                                                                                                                                       >  1:1309566/1‑151 (MQ=255)
            aaaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCtctc                                                                                                                                        <  1:1044984/150‑1 (MQ=255)
             aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGaaaaaatactcc                                                                                                                                       >  1:205433/1‑143 (MQ=255)
             aaaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCca                                                                                                                                      <  1:584891/151‑1 (MQ=255)
              aaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGaaaaaaaacccc                                                                                                                                       >  2:333530/1‑142 (MQ=255)
              aaaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGaaaaaacaatccc                                                                                                                                      >  1:989577/1‑142 (MQ=255)
                aGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCcaaca                                                                                                                                   <  1:817787/151‑1 (MQ=255)
                  cGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGaaaaaaacctcccacca                                                                                                                                  >  2:331392/1‑138 (MQ=255)
                   gAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAACTCTCcaacaa                                                                                                                                  >  2:1019332/1‑149 (MQ=255)
                    aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGCAAATGAAAAAACTCTCca                                                                                                                                      >  2:757479/1‑144 (MQ=255)
                    aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATAc                                                                                                                               <  1:719131/151‑1 (MQ=255)
                      ggTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCg                                                                                                                             <  2:717966/151‑1 (MQ=255)
                               aTGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGGAAAAACACACcaacaaaaa                                                                                                                               >  1:1045080/1‑137 (MQ=255)
                                   gTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGAAAAAGCTCTCCCACAATACCGCCTGCGTC‑TT‑C‑Ac                                                                                                             >  2:1165773/1‑151 (MQ=255)
                                   gTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑a                                                                                                              <  2:17566/150‑1 (MQ=255)
                                      tttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAAATCTCCCACCAAAACCCCCGCGTC‑TT‑c‑ccc                                                                                                            >  2:518357/1‑149 (MQ=255)
                                      tttGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑g                                                                                                          <  2:862416/150‑1 (MQ=255)
                                          ccAGGCGCTGATGACTGACCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGt                                                                                                        <  1:729276/148‑1 (MQ=255)
                                          ccAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAAAACTCCCACAAAAACGCCCGCGCC‑TT‑C‑Acccg                                                                                                          >  1:863670/1‑145 (MQ=255)
                                           cAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGTGAAAAAACACTCCCACAAAAACCCCTGCGTC‑TT‑T‑ACC‑gggcagg                                                                                                    >  2:1006665/1‑146 (MQ=255)
                                           cAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAACCTCCCACAAAAACCCCCGCGTC‑TT‑C‑Acccgcgagg                                                                                                     >  2:111210/1‑144 (MQ=255)
                                                gcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAACGCTCCCACAAAAACGCCCCCGGC‑TT‑C‑ATC‑CGGAA‑GGGGCtt                                                                                              >  2:716856/1‑151 (MQ=255)
                                                 cTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCt                                                                                               <  2:777482/149‑1 (MQ=255)
                                                         tGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGATGGAAAAACCCTCCAACAAAACCCCCCGCGTC‑TTCC‑CCC‑GGAAG‑GGGCctttcttt                                                                                         >  1:626237/1‑148 (MQ=255)
                                                          gAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTggc                                                                                     <  2:1070687/150‑1 (MQ=255)
                                                           aaCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGTAAAAACTCTCCCACAATACCCCCCGCGCCTTT‑T‑CCC‑GGTAA‑GGGCCtttttttttggg                                                                                    >  1:837567/1‑142 (MQ=255)
                                                             cTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGGAAAAAATCTCCCACAAAAACCCCCGCGTC‑TT‑T‑CCC‑GGTAA‑GGGGCTTTCtttt                                                                                        >  1:77575/1‑144 (MQ=255)
                                                              tGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAATCCCCCACCAAAAACCCCGCGTC‑CT‑CTCCC‑GGAAA‑GGGCCTTT‑ATTTtgcg                                                                                    >  1:1159697/1‑147 (MQ=255)
                                                                tATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGATGAAAAAAATCTCCCACAAAACCGCCCGCGGC‑TT‑T‑ACC‑GGGAA‑GGGCCTTTTAttttgggggcg                                                                                >  2:87892/1‑141 (MQ=255)
                                                                             gcgcGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGGAAAAAATCTCCAACAAAAACACCCCCGGC‑TT‑C‑ACC‑CGGAA‑GGGGCTTTTATTTTGCGGGCGGCTTATT‑TGccc                                                                   >  1:194494/1‑148 (MQ=255)
                                                                             gcgcGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGGTGAAAAAACTCTCCAACAATAACCCCCGCGTC‑TT‑C‑ACC‑CGGAA‑GGGGCTTTCATTTTGGGGCAGGCTttttttccccg                                                                  >  2:827054/1‑143 (MQ=255)
                                                                                cgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGAGATGAAAAAGAGCTCCAACAACAACGCCCTCGTC‑TT‑C‑ACC‑CGGAA‑AGGCCTTTTAATTTGGGGCAGGCTtttttaccc                                                                    >  2:766219/1‑140 (MQ=255)
                                                                                 ggATCCCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATAATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAg                                                              <  2:82043/149‑1 (MQ=255)
                                                                                 ggATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCt                                                            <  2:1084500/151‑1 (MQ=255)
                                                                                         aCGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGGGAAAAAACACTCCCACAAAAACGCCCGCGTC‑TT‑T‑CCC‑GGGAA‑GAGCCTTTTCTTTTGCGGCGGTCTTTTT‑CTCCCCGGCGCTACCg                                                        >  2:748741/1‑147 (MQ=255)
                                                                                          cGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGTTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATCCCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTAc                                                   <  2:1031814/151‑1 (MQ=255)
                                                                                              ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAATCTCCCACAAAAACGCCCGCGGC‑TT‑T‑ACC‑GGTAA‑GGGGCTTTCTTTTTTGGGCGGGCTTTTT‑CGCCCCGAAACAAccaggcaagg                                                 >  1:288858/1‑141 (MQ=255)
                                                                                              ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGt                                                <  1:1224189/150‑1 (MQ=255)
                                                                                              ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATAAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGt                                                <  1:1163363/150‑1 (MQ=255)
                                                                                               cGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGGTGAAAAAAAACTCCCACAAAAAACCCTGCGTC‑TT‑C‑ACC‑CGCAA‑AGGGCTTTTTTTTTGGGGCAGGCTTATT‑TTTCCCGAAGCTACCGGCCAc                                                   >  2:746958/1‑146 (MQ=255)
                                                                                                ggCAGGTGATATCGGGGTGGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGt                                                <  1:494317/148‑1 (MQ=255)
                                                                                                     gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGtt                                          <  1:163800/149‑1 (MQ=255)
                                                                                                     gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTa                                        <  2:1129886/151‑1 (MQ=255)
                                                                                                       gATATGGGGGTTGTGGGTAGTAAAGTGGGCTTTTTGGGGGGAATAATAAAAAAGTTCTCCAAAAATACCCCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTTCCGGCTACTGTCTGGTTTAt                                       <  1:1340673/150‑1 (MQ=255)
                                                                                                       gATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTAt                                       <  2:1184190/150‑1 (MQ=255)
                                                                                                        atatCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTAttt                                     <  1:821841/151‑1 (MQ=255)
                                                                                                         tatCGGGGTTGGTGGTCGGGAAGCCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTAt                                       <  1:534931/148‑1 (MQ=255)
                                                                                                          atCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCACCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTca                                   <  2:1335190/151‑1 (MQ=255)
                                                                                                             ggggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGAAAAAAAACTCCCACAAAAACGCCCGCGTC‑TT‑C‑ACCAGGAAA‑GGGCC‑TTCATTTTGCGGCAGGCTTTTT‑TGCCCGGAAGCTACCCGCTAAGGGGTGGGttttttt                                    >  1:135092/1‑146 (MQ=255)
                                                                                                                      ggtCGTGAAGACGGCTTTATGGGGGGTATGATAAAAAGGCTATCCAGCATTGCCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGc                       <  2:499562/151‑1 (MQ=255)
                                                                                                                         cGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTa                     <  1:924725/150‑1 (MQ=255)
                                                                                                                            gAAGTCGGCTTTATGGCGGGGGGGATGAAAAAACTCCCCCAAAAAAACGCCCGCGTC‑TT‑C‑CCC‑GGGAAGGGGCCTTTC‑TTTTGCGGCGGGCTTATTCCCCCCGGAGGCAACC‑GCCACGGGCTGGTTTTTTTTTCCAAAGCAATACTaaaa                  >  1:728990/1‑150 (MQ=255)
                                                                                                                                  ggCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCAc            <  1:281084/150‑1 (MQ=255)
                                                                                                                                       tATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTAt       <  1:450176/150‑1 (MQ=255)
                                                                                                                                          ggCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATggg    <  1:796564/150‑1 (MQ=255)
                                                                                                                                          ggCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATggg    <  1:575840/150‑1 (MQ=255)
                                                                                                                                           gCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGtt  <  2:1052757/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTC‑TT‑C‑ACC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATT‑CGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTT  >  NC_000913/1842769‑1843058

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