Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091332,6530TC43.4% 99.9 / 30.5 77S613P (TCC→CCC) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (38/5);  new base C (0/33);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 2.79e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTT  >  NC_000913/32517‑32799
                                                                                                                                              |                                                                                                                                                  
cggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTAc                                                                                                                                      >  2:425861/1‑151 (MQ=255)
 ggcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACt                                                                                                                                     >  1:1126928/1‑151 (MQ=255)
      gAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCg                                                                                                                                  >  2:982657/1‑149 (MQ=255)
            tATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGt                                                                                                                           >  2:675911/1‑150 (MQ=255)
               cGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAAc                                                                                                                        >  2:502667/1‑150 (MQ=255)
                ggTCGGGGTATCGATTTCGACTACTGTTGCGTACCCGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAACCCTTTATGGTTAA‑CTGTAC‑CCCGG‑AAC‑CC‑GTCCCCC‑CCGACTACGACCCTTCCGACCGCCTCTACTTCGAGCCGGTAAc                                                                                                                        <  1:617647/149‑1 (MQ=255)
                ggTCGGGGTATCGAATTCGACTCCTTTTGCGTACACCCCTCGCTGGCGCTGCGGGAAGACGGTTACGAACCCATTATGTTTAA‑CTGTAC‑CCCGG‑AAA‑CC‑GTCCCCC‑CCGACTACGCCACTCCCGACCGCCTCTACTTCGAGCCGGTaa                                                                                                                         <  1:884859/148‑1 (MQ=255)
                     ggttATCGATTTCGACTACTGTTGCGTCCCCGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAACCCATTATGGTTAA‑CTGTAC‑CCCGG‑AAA‑CC‑GTCTCCC‑CCGACTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGa                                                                                                                  <  1:308480/147‑1 (MQ=255)
                      gttaccaaaTTCGCCTCCTGTGCCGTCCCCGCCTCGCTGGCGCTGCGCGAAGACGTTTACGAAACCTTTATGTTAAA‑CT‑TAA‑CCCGGAAAA‑CC‑GTCCCCC‑CCACCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAg                                                                                                                <  2:1063196/144‑1 (MQ=255)
                               tccgccttcttttGCGTCCACTCCTCGTGGGCGCTGCGCGAAGGCTGTTCCGAACCCATTTTGGTTAA‑CT‑TAA‑CCCCG‑AAA‑CCGGTCTCCA‑CCGACTACGACTCCCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTg                                                                                                        <  2:309220/139‑1 (MQ=255)
                                 cGACTCCTGTTGCTACCCCCCCTCCCTGGCGCGGCGCGAAGACGGTTACAAAACCATTATGTTTAA‑CTGTAC‑CCCGA‑AAA‑CG‑GTCCCCC‑CCGACGACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGa                                                                                                      <  2:1237851/150‑1 (MQ=255)
                                    ctactGTTGCGTACACGCCTCGCTGGCGCTGCGCGAATACCGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGCGCTGGAAAt                                                                                                   >  2:785195/1‑150 (MQ=255)
                                    ctactGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATc                                                                                                  >  2:407706/1‑151 (MQ=255)
                                    ctactGTGGCGTCCACCCCTCGCGGGCGCTGCGGGAAGACGGTTACAAACCCTTTATGGTTAA‑CTGTAC‑CCCGG‑AAC‑CC‑GTCTCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGaaa                                                                                                    <  1:853923/149‑1 (MQ=255)
                                     tcctGTTGCGTCCCCCCCTCGCGGGCCCGGGGCGAAGCCGGTTACAAACCTATTATGGTAAA‑CTGT‑A‑CCCGG‑AAACCC‑GTCTCCC‑CCGACTACGCCACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGaaa                                                                                                    <  1:444941/146‑1 (MQ=255)
                                     tactGTGGCGTACCCGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAACCCATTATGTTAAC‑CTTTAC‑CCCGA‑AAA‑CC‑TTCCCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGaaa                                                                                                    <  2:812796/148‑1 (MQ=255)
                                      actGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGt                                                                                                >  2:704300/1‑151 (MQ=255)
                                        tttggcttcccccccTCGTGGGCGCTGCGCGAAGACGTTTACAAAACCATTATGTTTAA‑TTGT‑A‑CCCGG‑AAACCC‑GTCCCCA‑CCGACTACGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGt                                                                                                <  2:1169906/136‑1 (MQ=255)
                                         gTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCg                                                                                             >  1:224201/1‑151 (MQ=255)
                                          ttGCGTCCCCCCCCCGCTGCCGCTGCGCGAGGGGGATTCCAAAACCTTTATGTTTAC‑CTGT‑A‑CCCGG‑AAACCC‑GTCCCCC‑CCCCCTCCGCAACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGc                                                                                              <  1:73609/149‑1 (MQ=255)
                                          ttGCGTCCACGCCCCGCGGGCGCTGCGCAAAGCGGGTTCGAAAACCATTATGGTTAA‑CTGTAC‑CCCGG‑AAC‑CC‑GTCCCCC‑CCGATTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTgg                                                                                                       <  1:125564/140‑1 (MQ=255)
                                          ttGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGt                                                                                            >  2:484649/1‑151 (MQ=255)
                                          ttGCGAACCCCCCTCGCGGGCGCTGCGCGAAGACGTTTACAAACCCATTATGTTAAC‑CTGT‑A‑CCCGG‑AAACCC‑GTCTCCC‑CCGACTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTg                                                                                               <  2:510413/148‑1 (MQ=255)
                                             cGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATc                                                                                         >  2:48567/1‑151 (MQ=255)
                                              gTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                        >  2:453570/1‑151 (MQ=255)
                                                acacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                        >  2:1084152/1‑149 (MQ=255)
                                                acacGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                        >  1:332536/1‑149 (MQ=255)
                                                 ccccccTCGCGGGCGTGGCGCGAAGCGGTTTACAAACCCATTTTGTTTAA‑CTGAAC‑CCCGA‑AAC‑CC‑CTCCCCA‑CCGATTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCgaga                                                                                     <  2:1081729/147‑1 (MQ=255)
                                                   cgccgcGCTGGCGCTGCGCGAAGAGGGTTACGAAACCATTATGGTTAA‑CTGAAC‑CCGGG‑AAC‑CC‑GTCCCCC‑CCGACTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGaa                                                                                    <  1:323688/150‑1 (MQ=255)
                                                     ccTCGCTGGCGCTGCGCAAAGACGTTTACAAACCCATTAT‑GTTAA‑CCGGAACCCCGG‑AAC‑CC‑GTCTCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc                                                                                  <  2:184934/150‑1 (MQ=255)
                                                     ccTCGCGGGCGCTGCGCGAAGCGGGTTACAAACCCTTTTTGTTTAA‑CTGT‑A‑CCCGG‑AAACCC‑GTCCCCC‑CCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTAt                                                                                          <  2:367090/142‑1 (MQ=255)
                                                     ccTCGCGGGCGCTGCGCAAAGACGGTAAC‑AAACCCTTTTTGTTAA‑CTTGAA‑CCCCG‑AAACCC‑TTCCCCC‑CCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGcc                                                                                 <  2:827667/151‑1 (MQ=255)
                                                         gCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGaa                                                                              >  1:1150115/1‑150 (MQ=255)
                                                          cTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGaa                                                                              >  1:781170/1‑149 (MQ=255)
                                                           tGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAgg                                                                           >  1:339309/1‑151 (MQ=255)
                                                           tGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAgg                                                                           >  2:924259/1‑151 (MQ=255)
                                                            ggcgcggcgcG‑AGGCGGTTCCAAAAACCTTTAGGTTAACCTGTAC‑CCCGA‑AAC‑CC‑GTCTCCC‑CCAGTTACAACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGc                                                                          <  2:882483/151‑1 (MQ=255)
                                                            ggCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGc                                                                          >  2:492072/1‑151 (MQ=255)
                                                                 tGCGCAAAGACGATAACAAAACCATTATGGTTAA‑CT‑CAC‑ACCGG‑TAA‑CC‑CGCTCCACCCGACCACTCCACTGCCCTCCGCCTCTCCTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAgg                                                                           <  2:8751/145‑1 (MQ=255)
                                                                   cgggAAGCCGTTTACAAACCCATTATGGTTAC‑CTGTAC‑CCGGG‑AAC‑CC‑GTCCCCC‑CGGCCTAGACCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTAt                                                                     <  2:39989/146‑1 (MQ=255)
                                                                    gcgAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGt                                                                  >  1:948343/1‑151 (MQ=255)
                                                                    gcgAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGt                                                                  >  1:413236/1‑151 (MQ=255)
                                                                         ggcggtgccGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAC‑CC‑GTCTCCA‑CCGACTACGCCACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGt                                                             <  1:72816/143‑1 (MQ=255)
                                                                           cGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAc                                                           >  2:1046512/1‑151 (MQ=255)
                                                                            gtttaaaaacccATTATGGTAAA‑CTGTAC‑CCCGA‑AAC‑CC‑GCCCCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTa                                                            <  2:101031/139‑1 (MQ=255)
                                                                              ttACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggc                                                        >  2:1238374/1‑151 (MQ=255)
                                                                               tACAAAACAATTATGGTTAC‑GTG‑CA‑CCCGG‑AAACCC‑GTCCCCC‑CCGACTACGCCCCTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcgg                                                         <  2:436333/149‑1 (MQ=255)
                                                                                  gAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTc                                                    >  2:1106391/1‑151 (MQ=255)
                                                                                   aaccctttttgttaaa‑CTGAAC‑CCCGA‑AAC‑CC‑GTCTCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcg                                                       <  1:1014374/133‑1 (MQ=255)
                                                                                   aacccATTATGGTTAA‑CTGTAC‑CCCGA‑AAC‑CC‑GTCCCCC‑CCGACTCCGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTc                                                    <  1:433082/147‑1 (MQ=255)
                                                                                    aaCAATTATGG‑TAA‑CTTGTA‑CCCGG‑AAACCC‑GTCCCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAg                                                  <  1:360121/151‑1 (MQ=255)
                                                                                           aTGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGc                                           >  2:1189644/1‑151 (MQ=255)
                                                                                              gTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTg                                         >  2:176799/1‑150 (MQ=255)
                                                                                               ttAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGAcccc                                             >  2:616698/1‑145 (MQ=255)
                                                                                               gttaacctgaa‑CCCGG‑AAACCC‑GCCCCCC‑CCGCCTCCGACCTTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGAcccc                                             <  2:975688/138‑1 (MQ=255)
                                                                                                tAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa                                       >  2:185647/1‑150 (MQ=255)
                                                                                                 aa‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaaa                                      >  1:81258/1‑150 (MQ=255)
                                                                                                    cTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTg                                   >  2:1208771/1‑151 (MQ=255)
                                                                                                    atctgaaacgagaaaa‑CC‑TTTTCCC‑CAGACTACGTCACTCCGGACCGCCTCTTCTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa                                       <  1:1095535/136‑1 (MQ=255)
                                                                                                     tGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTgg                                  >  1:135392/1‑151 (MQ=255)
                                                                                                      tttacccggg‑AAC‑CC‑GCCCCCC‑CCGACTCCGACCTTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa                                       <  2:92927/138‑1 (MQ=255)
                                                                                                         accccGA‑AAC‑CC‑GTCCCCC‑CCGACTCCGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgc                                 <  1:1160497/147‑1 (MQ=255)
                                                                                                         accccGA‑AAC‑CC‑GTCCCCC‑CCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgc                               <  2:1226889/149‑1 (MQ=255)
                                                                                                         a‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcg                              >  1:537986/1‑151 (MQ=255)
                                                                                                             ccagaaaccc‑GCCTCCC‑CCGACTACCACCTTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTg                                   <  1:387900/137‑1 (MQ=255)
                                                                                                               g‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcgcgcg                          >  2:695635/1‑150 (MQ=255)
                                                                                                                 aaa‑CC‑GTCCCCC‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGgcgcgcgc                           <  1:14219/148‑1 (MQ=255)
                                                                                                                   accc‑GTCCCCC‑CCGACTCCGCCACTCCCGCCCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTg                       <  2:75764/149‑1 (MQ=255)
                                                                                                                           tcCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctg                >  1:800125/1‑151 (MQ=255)
                                                                                                                            ccA‑CCGCCTCCGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAg                   <  2:494429/147‑1 (MQ=255)
                                                                                                                                ccGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCCCACCCCGCTAAAACTGGCGCGCACGCTCGAAGCCGCTgg            >  1:575391/1‑151 (MQ=255)
                                                                                                                                    cTACAACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGctgctg             <  2:24738/146‑1 (MQ=255)
                                                                                                                                      aCGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTAcc      >  2:948100/1‑151 (MQ=255)
                                                                                                                                       cGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCCCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTAcc      >  2:138896/1‑150 (MQ=255)
                                                                                                                                          cacTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGtt  >  2:689190/1‑151 (MQ=255)
                                                                                                                                          cacTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTa        <  2:239168/145‑1 (MQ=255)
                                                                                                                                          cacTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGtt  <  1:441347/151‑1 (MQ=255)
                                                                                                                                             ttCCGACCGCCTCTACTTCGAGCCGGTAACTCCGGAAGATGTGCTGGAAATCGGGCGTATCGAGAACCCGAAAGGCGTAACCGTCCGGTACGGCGGGCCGACCCCGCTGAAACGGGGCGGCGCCCGGGGAGCCGCTGGCGTACCGGt   >  2:673961/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAA‑CTGTAA‑CCCGG‑AAA‑CC‑GTCTCCA‑CCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTT  >  NC_000913/32517‑32799

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: ‑