Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,465,733	A→G	Y150H (TAC→CAC)	treB ←	trehalose‑specific PTS enzyme IIBC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,465,733	0	A	G	100.0%	197.1 / NA	57	Y150H (TAC→CAC)	treB	trehalose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/13);  total (44/13)

CGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTC  >  NC_000913/4465598‑4465868                                                                                                                                         |                                                                                                                                        cGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCAtt                                                                                                                                     >  1:1910806/1‑140 (MQ=255)      aGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTgcgc                                                                                                                                <  1:133914/140‑1 (MQ=255)      aGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGGGCATTTgcgc                                                                                                                                >  2:19379/1‑140 (MQ=255)       gCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGcc                                                                                                                               <  2:1271929/140‑1 (MQ=255)            aTACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCGCAGAAAATCGTAGATCATTTGCAGGGAA‑GGAGGCATTTGCGCCAGCGt                                                                                                                          >  2:649651/1‑140 (MQ=255)             tACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTc                                                                                                                         >  1:1316065/1‑140 (MQ=255)              aCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCt                                                                                                                        <  2:1460024/140‑1 (MQ=255)                  aGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGAcc                                                                                                                    >  1:1505402/1‑140 (MQ=255)                    gATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGt                                                                                                                  <  2:164994/140‑1 (MQ=255)                     aTCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGtt                                                                                                                 >  1:1680596/1‑140 (MQ=255)                     aTCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTAGCAGGGAA‑GGGGGCCTTTGCGACAGCGTCTGACCGtt                                                                                                                 >  1:639838/1‑140 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGGAGGGAA‑GAGTGAATTTGCGCCGGCGTCTGACCGTTGc                                                                                                               >  1:664807/1‑140 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGa                                                                                                                      <  2:1100410/133‑1 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGa                                                                                                                      >  1:1100410/1‑133 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGc                                                                                                               >  2:1744662/1‑140 (MQ=255)                       cGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGCAGATCGTTTGCAGGGAG‑GGGGGCATGTGCGACAGCGTCTGACCGTTGc                                                                                                               >  2:753288/1‑140 (MQ=255)                              ccgccCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATggg                                                                                                        >  1:1773815/1‑140 (MQ=255)                                gccCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCa                                                                                                      <  1:410214/140‑1 (MQ=255)                                  ccATTTTTTTCACCGGTGACCATCAAATAACGACCGGAAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCaaa                                                                                                    <  1:649651/140‑1 (MQ=255)                                  ccATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCaaa                                                                                                    <  2:1773815/140‑1 (MQ=255)                                  ccATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCaaa                                                                                                    <  2:1253812/140‑1 (MQ=255)                                    aTTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGAAGAAGAAGATCGCTTCACCGGTCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGGGCATTTGCGCCAGCGTCTGACCGTTGCTCAAGGGCAAATc                                                                                                  >  1:380923/1‑140 (MQ=255)                                      ttttttCACCGCTGACCAGCAATTACCGATCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGc                                                                                                <  1:753288/140‑1 (MQ=255)                                          ttCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGAt                                                                                            <  2:1505402/140‑1 (MQ=255)                                              ccGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCATCGTCTGACCGTAGCCCATGGGTAAATCGCCGATcaca                                                                                        >  1:786734/1‑140 (MQ=255)                                                gCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGc                                                                                                >  2:1451369/1‑130 (MQ=255)                                                gCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGc                                                                                                <  1:1451369/130‑1 (MQ=255)                                                                cgaccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCACGGGCAAAGCGCGGAGCACATTACGAAAACCGAGGATc                                                                      >  1:574284/1‑140 (MQ=255)                                                                  accgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGAGCAAATCGCCGCTCCCATTGCTAAAACCGAGGATCaa                                                                    >  1:780266/1‑140 (MQ=255)                                                                   ccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCaaa                                                                   >  1:1521560/1‑140 (MQ=255)                                                                    cgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCCAAACAAATCGTATATCTTTTGCGGGGAG‑GGGTGCATTTGCGCCAGCGTCTGCCCGTTGCTAAGGGGAAAACCGCCGATCACATTGCAAAAACCGAGGATAAAAc                                                                  <  2:639838/140‑1 (MQ=255)                                                                    cgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAAc                                                                  >  1:1863468/1‑140 (MQ=255)                                                                       cAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAAccgc                                                               >  1:1428497/1‑140 (MQ=255)                                                                                   gaccGCTTCACCGATCACCAAATAAAAATCGTAGATCGTGTGGGGGGAG‑GGGTGCTTTTGCGCCAGCGTCTGACCGTTGCTAAGGGGCAAATCGCCGATCACATTGCAAAACCCGAGGACCAAACCGCCGCTAATCAAcg                                                   <  2:786734/137‑1 (MQ=255)                                                                                    aTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATGGCCGATCAGATTCCGACAACCGAGGGTCAAACCGCCGCTAATCAAcgc                                                  >  2:531319/1‑140 (MQ=255)                                                                                    aTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAAcgc                                                  >  1:1189523/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCTAAAACCGAGGACCAAACCGCCACTAATCAACGCGGgc                                              >  1:767507/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                              >  2:1519755/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                              >  2:1733148/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                              >  1:1327547/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGgc                                              >  1:1361903/1‑140 (MQ=255)                                                                                        cTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCAATCACATTCCGACAACGGAGGATCAGACCGTCGCTAATCAACGCGGgc                                              >  2:609252/1‑140 (MQ=255)                                                                                           cACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCACACAGCCGCTAATCAACGCGGgcagc                                           >  1:695170/1‑140 (MQ=255)                                                                                            aCCGATCAACCACAGCAAATCGTAGATCGCTTGCAGGCAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAGCCGCCGCTAATCAACGCGGgcagca                                          <  1:531319/140‑1 (MQ=255)                                                                                              cGATCAACCACAAAAAATCGTAGATCG‑TTGGAGGGAAGGGGTGCATTGGCGCCAGCGTCTGCCCGTTGCTCAGGGGAAAAGCGCCGATCACATGGCAAAACCCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAAc                                        <  2:767507/140‑1 (MQ=255)                                                                                              agaccAACCACAAAAAATCGTAGATCGTTTGAAGGGAA‑GGGTGCTTTTGCCCCAGCGTCTGCCCGTTGCTCAGGGGAAAACCGCCAACCACATTGCAAAAACCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAAc                                        <  2:1619413/136‑1 (MQ=255)                                                                                                aTCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACgg                                      >  1:1048044/1‑140 (MQ=255)                                                                                                aTCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACgg                                      >  1:23017/1‑140 (MQ=255)                                                                                                 tCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACggg                                     >  1:1813208/1‑140 (MQ=255)                                                                                                      cacaGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaagaa                                >  1:1480734/1‑140 (MQ=255)                                                                                                        caGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATGAACGCGGGCAGCAACGGGAagaaga                              >  1:660350/1‑140 (MQ=255)                                                                                                          gAAAATCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATc                            >  2:1510644/1‑140 (MQ=255)                                                                                                           aaaaTCGTAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATct                           >  1:1619413/1‑139 (MQ=255)                                                                                                                  tAGATCGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATAAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaa                    >  1:1435883/1‑140 (MQ=255)                                                                                                                       cGTTTGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag               >  1:1231996/1‑140 (MQ=255)                                                                                                                        gTTTGCAGGGAG‑GGGTGCATTTGCGCCAGCGTCTGCCCTTTGCTCATGGAAAAATCGCCGACCACATTGAAAAACCCGAGAAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgaga              <  2:1231996/140‑1 (MQ=255)                                                                                                                         tttGCAGGGAA‑GGGTGCTTTGGCGCCAGCGTCTAGCGTTTCCTCAGGGGCAAATCCCCAACCACATTACAAAACCCGAGCATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag             <  2:380923/140‑1 (MQ=255)                                                                                                                         tttGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGTGGGCAGCGACGGGAGGAAGATGACCGCGAAATgagag             >  1:593153/1‑140 (MQ=255)                                                                                                                         tttGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag             >  1:1042659/1‑140 (MQ=255)                                                                                                                         tttGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag             >  1:1107700/1‑140 (MQ=255)                                                                                                                         tttGCAGGGAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGAGGGCAGCAACGCGAAGAAGATGACCGCGAAATgagag             >  1:486437/1‑140 (MQ=255)                                                                                                                         ttggaagggag‑ggGTGCTTTGGCGCCGGCGTCTGACCGTTGCTAAGGGGAAAATCGCCGATCACATTGCAAAAACCGAGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag             <  2:1428497/135‑1 (MQ=255)                                                                                                                         ttggCAGGGAG‑GGGTGCATTGGCGCCAGCGTCTGCCGGTTGCCCAGGGGAAAACCGCCGATCACATGGCAAAACCCGAGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag             <  2:1480734/137‑1 (MQ=255)                                                                                                                              ggggAG‑GGGTGCTTTTGCGCCGTCGTCTGACCGTGTCTCTGGGGAAAACCGCCGATCACATTGCAAAACCCGAGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCaa        <  2:593153/139‑1 (MQ=255)                                                                                                                               gggAA‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAAc       >  1:188727/1‑140 (MQ=255)                                                                                                                                  aa‑GGGTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGc    >  2:1647476/1‑140 (MQ=255)                                                                                                                                     gggTGCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTc  >  1:1607086/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        CGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAA‑GGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTC  >  NC_000913/4465598‑4465868

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