Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	40.8%	127.2 / 26.8	76	L79V (TTA→GTA)	maeB	malic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (40/5);  new base C (0/31);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.97e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCT  >  NC_000913/2577994‑2578283                                                                                                                                                             |                                                                                                                                               tggttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTaa                                                                                                                                                <  2:1436251/150‑1 (MQ=255)       gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAg                                                                                                                                         >  1:893848/1‑150 (MQ=255)       gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAg                                                                                                                                         >  1:1006862/1‑150 (MQ=255)       gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGc                                                                                                                                        >  2:77902/1‑151 (MQ=255)        aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAg                                                                                                                                         >  1:1252821/1‑149 (MQ=255)        aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGc                                                                                                                                        >  2:314881/1‑150 (MQ=255)        aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGc                                                                                                                                        >  2:567032/1‑150 (MQ=255)         gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCa                                                                                                                                       <  1:1451968/150‑1 (MQ=255)         gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑c                                                                                                                                     <  1:666550/151‑1 (MQ=255)              gcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑ccgccg                                                                                                                                >  1:1096001/1‑151 (MQ=255)                 aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTAc                                                                                                                             >  2:231742/1‑151 (MQ=255)                         aataatTTTGTCCGGGTCGAGTTCGTCAACTTAAATGTCAAATACATCAACCCCGGCGATTTTCTT‑AAACAGAACGCCCTTGCCTC‑CCA‑CCCCGGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAC‑CCCCAGCA‑CCGCCGTACCGTTagag                                                                                                                     <  2:1236046/151‑1 (MQ=255)                           tAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTag                                                                                                                       >  2:639113/1‑147 (MQ=255)                              aTTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCAc                                                                                                                >  1:1569451/1‑151 (MQ=255)                               tttGTCCGGGTCGTTTCCCTCACCTTCAATGTAAAATCCATCACCCCCGGCGATTTTCTT‑AAACAGACCCCCCTTCCCTC‑CCA‑CCCCCGT‑TT‑T‑GCCGG‑CCGGCGCGCCATTTTCGCCAAC‑CCCCGGCT‑CCGCCGTACCGTTAGAGATCAc                                                                                                                <  1:1312503/150‑1 (MQ=255)                               tttGCCCGGTTCGATTCCGTCAACTTCAATGTAAAATACATCACCCCCGGCGATTTTCTA‑AAACAGACCCCCCTTCCCTC‑CCA‑CCCCGGT‑TT‑T‑GCCGG‑CGGGCGCGCAATTGTTCCCTAC‑CCCCAGCC‑CCGCCGTACCGTTAGAGATCa                                                                                                                 <  1:92568/149‑1 (MQ=255)                                    cggggTCGAGTTCGTCAACTTCAATGTCAAATACATCAACCCCGGCGATTTTCTT‑AAACAGACCGCCCTTGCCTT‑CCA‑CCCCCGG‑TT‑T‑GCCTG‑CCGGCCCGCCAATGTTCCCAAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCg                                                                                                              <  1:93656/145‑1 (MQ=255)                                           gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                    >  1:965014/1‑150 (MQ=255)                                                    aaCTTAAATTTAAAATACATCAACCCCGGCAATTTTCTA‑AAACAAACCCCCTTTCCCTT‑CCA‑‑CCCCGG‑TT‑TGGCCTG‑CCAGCGCGCAAATGTTCCCAAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                               <  2:1237502/146‑1 (MQ=255)                                                          actgaaaaaaccATCAACCCGGGGTATTTTCTA‑AAACACACCCCCCTTCCTTC‑CCA‑CCACCGT‑TT‑G‑CCCCG‑CCAGCCCCCCAATGGTGCTTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtata                                                                                     <  1:56450/139‑1 (MQ=255)                                                          aaTGTAAAATCCATCACCCCCGGC‑AATTTTTT‑AAAACAAACCCCCTTTCCTT‑CCC‑ACCCCGGTTT‑T‑GCCGC‑CCGGCCCGCAATTGTGCCCAAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                    <  2:927242/151‑1 (MQ=255)                                                            tgccAAATACAAAAACCCCGCCGATTTCCTA‑AAACAGACCCCCCTTCCCTC‑CCC‑CCCCCGT‑TT‑T‑GCCCC‑CCAGCGC‑CCCATTTTGCCCAACCCCCGGCC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                           <  1:1351051/139‑1 (MQ=255)                                                            tGTCAAATACATCAATCCCGGCGAATTTCTT‑‑AACAGAACCCCCTTGCCTT‑CCC‑ATACCGG‑TT‑T‑GCCCG‑CCCACCCGCCAATGTTGCCCAA‑CCCCAACACCCGCCG‑AACGGTAAAGATCCCCGGCACCAAGTTAACTCcggcgggtttttt                                                                                  >  1:1204090/1‑143 (MQ=255)                                                                aaacaccacAACCCCGCCGATTTTCTA‑AAACAGAAC‑CCCCTGCCTTCCCA‑CCCCCGT‑TT‑T‑GCCTG‑CCAGCGCGCCAATTTCGCCCAC‑CCCCAGCG‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAtt                                                                                   <  2:1513834/138‑1 (MQ=255)                                                                   tACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt                                                                           >  2:412002/1‑151 (MQ=255)                                                                     cATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTCGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                          >  1:660459/1‑150 (MQ=255)                                                                         aaccccGGAGAATTTCAA‑AAACAGACCGCCTTTCCCTT‑CCA‑CCCCGGT‑TT‑T‑GCCGG‑CCGGCGCGCAATTGTTGCCTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                        <  2:1080772/145‑1 (MQ=255)                                                                           tCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                   >  1:970262/1‑151 (MQ=255)                                                                            ccgggCAAATTTCTT‑AAACAAACCCCCTTTCCCTC‑CCA‑CCCCCGT‑TT‑T‑CCCTG‑CCAGCGCGCCATGGTTCCCAAC‑CCCCACCC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                    <  2:1361984/146‑1 (MQ=255)                                                                            cccGGCAATTTTCTA‑AAACAAACCCCCCTTCCCTC‑CCA‑CCCCGGT‑TT‑T‑CCCCC‑CCGGCGCGCCATTTTTGCCTAC‑CCCCAGCT‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                    <  1:898645/149‑1 (MQ=255)                                                                                ggcaaTTTTTTAAAAAAGACCGCCTTTCCCTC‑CCC‑CCCCCGT‑TT‑T‑CCCTG‑CCGGCCCCCCAATGTTCCCTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                               <  1:103068/148‑1 (MQ=255)                                                                                gCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtagaatcggaag                                                                >  2:855890/1‑138 (MQ=255)                                                                                  gattttttt‑AAAGAGAAC‑CCCTTTGCTT‑CCAGACACGTT‑TT‑G‑GCCGG‑CTAGCGCGCAATTTTTGCCAAC‑CCCCTGCT‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                              <  2:338014/146‑1 (MQ=255)                                                                                   aaTTTTTA‑AAACAG‑ACCCCCTTTCCCT‑CCC‑ACCCCGG‑TTGT‑CCCTA‑CCGCCGCGCCAGTTATCCCTAC‑CCCCAGCT‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                            <  2:1510178/150‑1 (MQ=255)                                                                                   aaTTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                           >  1:83219/1‑151 (MQ=255)                                                                                     tttCTA‑AAACAGACCTCCCTTCCCTC‑CCA‑TCCCGGG‑TT‑T‑GCCTC‑CCGGCGCACAAATGTCGCCTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                           <  1:802618/149‑1 (MQ=255)                                                                                      ttCTA‑AAACAGACCCCCTTTCCCTC‑CCA‑CCCCCGT‑TT‑T‑CCCTG‑CCGGCGCGCAAATGTTGCCTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                         <  2:1335574/150‑1 (MQ=255)                                                                                       tCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑TCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                       >  2:1301177/1‑151 (MQ=255)                                                                                         tt‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                     >  2:1389477/1‑151 (MQ=255)                                                                                           aaaacaaacCCCCCTTCCCTC‑CCA‑‑CCCCGG‑TT‑T‑TCCCGCCCGGCGCCCCAGTTTTGCCTAC‑CCCCGGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                      <  2:1472336/148‑1 (MQ=255)                                                                                           aaaaCAGACCCCCCTGCCCTC‑CCC‑CCCCCGT‑TT‑C‑CCCTC‑CCGGCGCCCAATTTTTCCCTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                    <  1:933398/150‑1 (MQ=255)                                                                                            aaaCAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                     >  2:501530/1‑149 (MQ=255)                                                                                            aaaCAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                    >  1:1261249/1‑150 (MQ=255)                                                                                             aaCAGAACTCCTTTGCCTT‑CAA‑CCACTGG‑TT‑T‑GCTTG‑CCCGCGCCCTATTGTTGCCCAC‑CCCCCAAA‑CCGCCGTACCGTTAGAGATCACCGCCCCCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                     <  2:167832/148‑1 (MQ=255)                                                                                                aGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                               >  1:169625/1‑151 (MQ=255)                                                                                                  accGCCTTTCCCTC‑CCA‑CCCCCGG‑TT‑T‑GCCTG‑CCAGCGCGCAATTGTGCCCTAC‑CCCAAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                              <  2:352543/148‑1 (MQ=255)                                                                                                  aaCGCCCTTGCCTT‑CCA‑TCACCGG‑AT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTCTATTTGTCGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                              >  1:669071/1‑150 (MQ=255)                                                                                                    cccccTTGCCTT‑CCA‑TCCCCGG‑TT‑T‑GCCTG‑CCAGCCCGCCATTGTTCCCAAC‑CCCCAGCC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACaa                                             <  2:1017106/147‑1 (MQ=255)                                                                                                      cccTTGCCTC‑CCA‑CCCCCGT‑TT‑T‑GCCGG‑CCAGCGCCCAAATGTGCCCAAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                           <  2:1070963/149‑1 (MQ=255)                                                                                                      cccTTCCCTC‑CCA‑CCCCGGG‑TT‑T‑GCCGG‑CCGGCGCGCAAATGTTCCCAAC‑CCCAAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTg                                         <  2:169646/151‑1 (MQ=255)                                                                                                        cTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGAGCg                                       >  1:1536424/1‑151 (MQ=255)                                                                                                         ttGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                      >  2:590139/1‑151 (MQ=255)                                                                                                         ttGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                      >  1:354038/1‑151 (MQ=255)                                                                                                          tGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCGGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGc                                     >  2:315176/1‑151 (MQ=255)                                                                                                             cTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                  >  2:646590/1‑151 (MQ=255)                                                                                                                 cccaccccggg‑TT‑T‑CCCTG‑CCAGCGC‑CCAATTGTGCCCAACCCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                 <  1:380219/141‑1 (MQ=255)                                                                                                                   cacccccGG‑TT‑T‑CCCGG‑CAAGCGCCCCAGGGTCCCCCAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                  <  1:253037/142‑1 (MQ=255)                                                                                                                     tCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTgg                            >  1:463968/1‑151 (MQ=255)                                                                                                                      cacttg‑TT‑G‑GCCTG‑CCGGCCCGCCATGGTCGCCTAC‑CCCCGGCT‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCt                              <  2:704485/143‑1 (MQ=255)                                                                                                                          gg‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAg                        >  1:1558625/1‑150 (MQ=255)                                                                                                                          gg‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                       >  1:815759/1‑151 (MQ=255)                                                                                                                          gg‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                       >  2:599204/1‑151 (MQ=255)                                                                                                                           g‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                       >  2:1159022/1‑150 (MQ=255)                                                                                                                                ttggccgcccAGCGCCCAATTGTTGCCTAC‑CCCCAGCC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAg                        <  2:1551221/140‑1 (MQ=255)                                                                                                                                t‑gcctg‑ccAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGc                   >  2:413241/1‑151 (MQ=255)                                                                                                                                    ctg‑ccAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                  >  2:1271718/1‑149 (MQ=255)                                                                                                                                    ctg‑ccAACGCGCCAATGCTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                  >  1:12141/1‑149 (MQ=255)                                                                                                                                        ccAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc             >  1:411740/1‑151 (MQ=255)                                                                                                                                              cgcCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc            >  1:1323033/1‑146 (MQ=255)                                                                                                                                               gcCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc            >  2:613070/1‑145 (MQ=255)                                                                                                                                               cccAATTTTTCCCAACCCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAg          <  2:45623/147‑1 (MQ=255)                                                                                                                                                ccATTGTTCCCTCC‑CCCAAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcg      <  1:1012233/150‑1 (MQ=255)                                                                                                                                                caaaTGTCGCCCAC‑CCCCACCG‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcg      <  2:1086054/148‑1 (MQ=255)                                                                                                                                                  aaTGTTGCCTAC‑CCCCACCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAg          <  2:454142/144‑1 (MQ=255)                                                                                                                                                   aTGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt  >  1:1217815/1‑151 (MQ=255)                                                                                                                                                    ttttGCCAAC‑CCCCAGCC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt  <  2:841600/148‑1 (MQ=255)                                                                                                                                                      tccccTAC‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAg          <  1:610212/137‑1 (MQ=255)                                                                                                                                                             |                                                                                                                                               TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAACAGAACGCCCTTGCCTT‑CCA‑TCACCGG‑TT‑T‑GCCTG‑CCAGCGCGCCAATGTTGCCTAA‑CCCCAGCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCT  >  NC_000913/2577994‑2578283

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