Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,554,471	T→C	intergenic (+110/‑105)	uxuB → / → uxuR	D‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,554,471	0	T	C	98.3%	209.2 / ‑3.9	60	intergenic (+110/‑105)	uxuB/uxuR	D‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (34/25);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

GCGAACGCGGTGCGCGCGAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAA  >  NC_000913/4554324‑4554614                                                                                                                                                    |                                                                                                                                                gCGAACGCGGTGCGCGCGAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGctt                                                                                                                                               <  2:759129/150‑1 (MQ=255) gCGAACGCGGTGCGCGCGAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTg                                                                                                                                              <  2:192865/151‑1 (MQ=255)  cGAACGCGGTGCGCGCGAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTga                                                                                                                                             <  2:785364/151‑1 (MQ=255)  cGAACGCGGTGCGCGCGAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTga                                                                                                                                             <  2:77256/151‑1 (MQ=255)             cgcgcgAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACgggg                                                                                                                                  <  2:480671/151‑1 (MQ=255)                 cgAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCg                                                                                                                              >  1:58122/1‑151 (MQ=255)                     tgtgtgGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTaa                                                                                                                           >  1:154117/1‑150 (MQ=255)                     tgtgtgGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTaa                                                                                                                           >  2:87687/1‑150 (MQ=255)                      gtgtgGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTaaa                                                                                                                          <  2:475837/150‑1 (MQ=255)                        gtgGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAAt                                                                                                                        <  2:667457/150‑1 (MQ=255)                          ggCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTg                                                                                                                      <  2:321493/150‑1 (MQ=255)                          ggCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTg                                                                                                                      <  1:109336/150‑1 (MQ=255)                             tGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCg                                                                                                                  >  2:401288/1‑151 (MQ=255)                               cgcgGTAACTAACTGATTACCCTGCAGACTTAGGGGGCGATGAAACTGATATTTGGTTGGCCAGTTTTCGTTTTTTTGCCCACCTGAACGGGCGAGATTCCAGTGCTAATGGTGTAGTTTGAGATTAATGGGGCCCGTAAGATTGccgggt                                                                                                                >  1:326832/1‑147 (MQ=255)                                    taactaacTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAgg                                                                                                           >  1:225394/1‑151 (MQ=255)                                     aactaactGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGc                                                                                                          >  2:447524/1‑151 (MQ=255)                                         aactGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCg                                                                                                        <  2:272160/149‑1 (MQ=255)                                                tACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCg                                                                                               >  2:777415/1‑151 (MQ=255)                                                  cccTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGtt                                                                                             >  2:306599/1‑151 (MQ=255)                                                    cTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTc                                                                                            >  2:271838/1‑150 (MQ=255)                                                        aGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                       >  1:79071/1‑151 (MQ=255)                                                                tGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTT‑AGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                              >  2:575472/1‑151 (MQ=255)                                                                 ggTCAATCAAACTGTAATTTGGTTGACCAGTTTTCGTTTTTTAGCCCACCTGCACGTGCCAACTTCCAGGGCTAATGGTATAGCTTGGGATTAACGGTGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCGGTATCCGGc                                                                              <  1:380238/151‑1 (MQ=255)                                                                 ggTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                <  2:556096/149‑1 (MQ=255)                                                                  gTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                              >  2:592190/1‑150 (MQ=255)                                                                     aatcaaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaaaa                                                                          <  1:317095/151‑1 (MQ=255)                                                                        caaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAtt                                                                       >  1:858265/1‑151 (MQ=255)                                                                        caaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAt                                                                        <  1:31933/150‑1 (MQ=255)                                                                         aaaCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                      <  1:686524/151‑1 (MQ=255)                                                                             tGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAtt                                                                  >  2:480731/1‑151 (MQ=255)                                                                              gATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAtt                                                                  <  2:642223/150‑1 (MQ=255)                                                                              gATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGAtt                                                                  >  2:26744/1‑150 (MQ=255)                                                                                      gTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCAcc                                                         >  2:20364/1‑151 (MQ=255)                                                                                       ttGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCt                                                        >  2:389951/1‑151 (MQ=255)                                                                                        tGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCt                                                        >  1:752629/1‑150 (MQ=255)                                                                                          aCCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                     >  1:490699/1‑151 (MQ=255)                                                                                          aCCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                     >  2:322118/1‑151 (MQ=255)                                                                                          aCCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAAc                                                     >  1:347233/1‑151 (MQ=255)                                                                                                    gTTTTTTTGCCCACCTGTACGTGCCACCTGCCAGTGCTAAGGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGCCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCCCCTAACGGACacaac                                            <  1:583409/150‑1 (MQ=39)                                                                                                    gTTTTTTTGCCCACCTGTACGTGCCAACTTTCAGTGCTAATGGTCTGGCTTGAGATTAACGGCGGCCGTAAAATTGCCCGATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaa                                             >  1:707286/1‑149 (MQ=38)                                                                                                            gCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGaa                                    <  2:47003/150‑1 (MQ=255)                                                                                                            gCCCACCTGAACGTGCCAACTTCCAGCGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCTCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGaa                                    >  1:614064/1‑150 (MQ=255)                                                                                                             cccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCt                                                        <  2:566611/129‑1 (MQ=255)                                                                                                             cccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCt                                                        >  1:566573/1‑129 (MQ=255)                                                                                                             cccACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAAt                                  >  2:269110/1‑151 (MQ=255)                                                                                                               cACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCt                                <  1:811244/151‑1 (MQ=255)                                                                                                                 ccTGTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGc                              >  2:210440/1‑151 (MQ=255)                                                                                                                    gTACGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGcc                             <  2:496829/149‑1 (MQ=255)                                                                                                                      aCGTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACct                         <  2:749525/151‑1 (MQ=255)                                                                                                                        gTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctct                       >  2:34388/1‑151 (MQ=255)                                                                                                                        gTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctct                       >  2:341179/1‑151 (MQ=255)                                                                                                                        gTGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACctct                       <  1:414703/151‑1 (MQ=255)                                                                                                                         tGCCAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTg                      >  2:604548/1‑151 (MQ=255)                                                                                                                           ccAACTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTgcg                    <  1:365486/151‑1 (MQ=255)                                                                                                                              aCTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCa                  >  2:434238/1‑150 (MQ=255)                                                                                                                              aCTTCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCa                  >  1:760230/1‑150 (MQ=255)                                                                                                                                ttCCAGTGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCaacg               >  1:422693/1‑149 (MQ=255)                                                                                                                                ttCCAGCGCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCaa                 >  1:256522/1‑149 (MQ=255)                                                                                                                                       gCTAATGGTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTa         <  2:392139/150‑1 (MQ=255)                                                                                                                                             ggTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGaa  <  1:625309/151‑1 (MQ=255)                                                                                                                                             ggTATAGCTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGa   >  1:729557/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                GCGAACGCGGTGCGCGCGAGTGTGTGGCTGCGCTGTAACTAACTGATTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAA  >  NC_000913/4554324‑4554614

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