Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	intergenic (‑1/+1)	gatC ← / ← gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,361	0	C	.	90.9%	159.9 / 4.1	44	pseudogene (2/442 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (25/15);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	2,173,362	0	C	.	90.9%	158.5 / 2.5	44	pseudogene (1/442 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base . (25/15);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACA  >  NC_000913/2173123‑2173544                                                                                                                                                                                                                                                    ||                                                                                                                                                                                       cgttgtgccgccAGTT‑GGTGTGCAAAACCATTGTGTGTCTACCGAACCCCAGGTATATCGCCCAGAATTAGGACCCCGGAG‑ATTAGGGTGGGGAACAGATTCCCGACATGCGCGGCGCCCGCAATCCCCAAAAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                                                                         <  4:706/233‑1 (MQ=255)            cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                                                                             >  3:257602/1‑238 (MQ=255)                           aCCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       >  4:78835/1‑214 (MQ=255)                           aCCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       >  4:104583/1‑214 (MQ=255)                           aCCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       <  3:78835/215‑2 (MQ=255)                           aCCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       <  3:104583/215‑2 (MQ=255)                                     tGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                   >  3:405779/1‑239 (MQ=255)                                       cGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                                                                      >  3:396818/1‑218 (MQ=255)                                                       gTGAT‑GCACATAATGCTGACA‑CCGAGAATTAGGGTGCGGACCAGTTTCCCACGATGCACGGCGACCCCCACCGCCACAAAGAAG‑CCAATGGTGCCAAGACCGCAAAACGGCAGCCCCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGAAAAAAACCAGCC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                     <  4:405779/235‑1 (MQ=255)                                                              cTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGa                                                                                                                            >  3:234083/1‑238 (MQ=255)                                                               tCATAATGATGACACCCGAG‑ATTAAGGTCCGGAACAGATTTCCACGATCCACGGCGACCGCCATCGCCCCGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                 <  1:111191/216‑1 (MQ=255)                                                                        tgaCACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGCATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                      >  3:389387/1‑234 (MQ=255)                                                                         gaCACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCAAGAAGAAAAGACGGACACAGGCCAATCAGGAACTc                                                                                                                                       >  2:111191/1‑216 (MQ=255)                                                                                             cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       >  2:136501/1‑150 (MQ=255)                                                                                             cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       <  1:136501/151‑2 (MQ=255)                                                                                             cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                         >  1:52865/1‑227 (MQ=255)                                                                                                     aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                     >  3:398526/1‑239 (MQ=255)                                                                                                             gttGCACGGCGACCGCCATCCCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAACCGGCGGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                             <  3:60568/237‑1 (MQ=255)                                                                                                                        aCCGCCATCGGCAATAAGAAGCCCGTTGGTGGAAAGACCGCAAAAGGCCAGCACTTAATTCCCCGGAACAC‑AAACAGCATTAAAAAGGGTGAGTGGAAAAAAAACCAGGC‑TTCCAGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                    <  4:398526/238‑1 (MQ=255)                                                                                                                         ccGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                 >  4:368756/1‑239 (MQ=255)                                                                                                                                      gaagaaG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                    <  3:368756/239‑1 (MQ=255)                                                                                                                                              ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       <  3:233665/103‑2 (MQ=255)                                                                                                                                              ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                                                                       >  4:233665/1‑102 (MQ=255)                                                                                                                                                        cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                   >  3:115291/1‑239 (MQ=255)                                                                                                                                                        cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                   >  3:283098/1‑239 (MQ=255)                                                                                                                                                               gCCAAAGGCGAGCCCCTGATT‑CCCGGGCCAC‑AAAACCGAAATAAAAAGGGGGGAGGGAAAAAAAACCGGCTTTGCCAAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                            <  4:115291/239‑1 (MQ=255)                                                                                                                                                                           caCCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                        <  4:27562/135‑1 (MQ=255)                                                                                                                                                                           caCCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                        >  3:27562/1‑135 (MQ=255)                                                                                                                                                                                gATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                               >  1:198390/1‑187 (MQ=255)                                                                                                                                                                                 aaaaagaaacacacaAAACAGAAATAAAAATGGTGAGTGGAATAAAAATCAGGC‑TTGCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg            <  1:103939/229‑1 (MQ=255)                                                                                                                                                                                          cacac‑aAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                                                   <  1:6213/125‑1 (MQ=255)                                                                                                                                                                                          cacac‑aAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                                                   >  2:6213/1‑125 (MQ=255)                                                                                                                                                                                           acac‑aAACACCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                                                               <  3:33410/160‑1 (MQ=255)                                                                                                                                                                                           acac‑aAACACCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                                                               >  4:33410/1‑160 (MQ=255)                                                                                                                                                                                                        ataaaaaTGGTGAGTGGAAAAAAAATCAGGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                     <  2:212028/91‑1 (MQ=255)                                                                                                                                                                                                        aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATAACCCGGGGCATTAGCAGCATCACTGCCGCCGTTTTTACAGCCAGCTGCAATaga  >  2:103939/1‑223 (MQ=255)                                                                                                                                                                                                                ggggAGTGGGAAAAAAATCAGCC‑TTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                <  4:215459/136‑1 (MQ=255)                                                                                                                                                                                                                  tGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                                                           >  1:212028/1‑94 (MQ=255)                                                                                                                                                                                                                   gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                                                                        >  2:214692/1‑48 (MQ=255)                                                                                                                                                                                                                   gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                                                                        <  1:214692/48‑1 (MQ=255)                                                                                                                                                                                                                   gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                     <  3:68520/67‑1 (MQ=255)                                                                                                                                                                                                                   gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                     >  4:68520/1‑67 (MQ=255)                                                                                                                                                                                                                   gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                             >  4:143408/1‑139 (MQ=255)                                                                                                                                                                                                                   gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                             <  3:143408/139‑1 (MQ=255)                                                                                                                                                                                                                    aGTGGGATAAAAATCAGGC‑TTGCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACACCTGACCGCCGAACTTCGCCTGTAAACGGCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                          <  3:376547/157‑1 (MQ=255)                                                                                                                                                                                                                          aTAAAAATCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                                                                      >  3:215459/1‑139 (MQ=255)                                                                                                                                                                                                                              aaaTCAGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAATCCATCCATGATGGGATTATTCACCCGTgg                                                >  4:376547/1‑157 (MQ=255)                                                                                                                                                                                                                                   aGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                    <  3:412720/52‑1 (MQ=255)                                                                                                                                                                                                                                   aGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                    >  4:412720/1‑52 (MQ=255)                                                                                                                                                                                                                                   aGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                <  3:66872/72‑1 (MQ=255)                                                                                                                                                                                                                                   aGGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                >  4:66872/1‑72 (MQ=255)                                                                                                                                                                                                                                    ggC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc           <  3:62599/188‑1 (MQ=255)                                                                                                                                                                                                                                    ggC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc           >  4:62599/1‑188 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                       GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGAT‑GCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAG‑CCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAC‑AAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACA  >  NC_000913/2173123‑2173544

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