Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,412,371	T→G	K289Q (AAA→CAA)	intR ←	Rac prophage; integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,412,371	0	T	G	100.0%	229.5 / NA	60	K289Q (AAA→CAA)	intR	Rac prophage; integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (33/27);  total (33/27)

CCCAAACTGCTTACCGAGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTTGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCC  >  NC_000913/1412226‑1412517                                                                                                                                                  |                                                                                                                                                   cccAAACTGCTTACCGAGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTAtt                                                                                                                                               <  2:338346/151‑1 (MQ=255)         gCTTACCGAGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCt                                                                                                                                         <  2:346559/149‑1 (MQ=255)             aCCGAGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTg                                                                                                                                   <  1:118425/151‑1 (MQ=255)              ccGAGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGc                                                                                                                                  <  1:801571/151‑1 (MQ=255)                 aGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTc                                                                                                                               >  2:642053/1‑151 (MQ=255)                      taaGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAag                                                                                                                            >  1:643241/1‑149 (MQ=255)                       aaGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAagag                                                                                                                          <  1:75452/150‑1 (MQ=255)                          gTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAgg                                                                                                                      >  1:811776/1‑151 (MQ=255)                            ttttGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGc                                                                                                                     <  2:195608/150‑1 (MQ=255)                               tGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGc                                                                                                                  >  2:662296/1‑150 (MQ=255)                                  tttCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCagt                                                                                                               >  1:608505/1‑150 (MQ=255)                                          acacacCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACt                                                                                                      <  2:918144/151‑1 (MQ=255)                                           cacacCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACt                                                                                                      <  2:681742/150‑1 (MQ=255)                                              acCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTc                                                                                                  <  2:649329/151‑1 (MQ=255)                                               cccTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCg                                                                                                 >  1:396800/1‑151 (MQ=255)                                                ccTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTc                                                                                                  <  1:631012/149‑1 (MQ=255)                                                ccTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCgg                                                                                                >  2:743899/1‑151 (MQ=255)                                                     gCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATccc                                                                                           <  1:374410/151‑1 (MQ=255)                                                       aTGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCgg                                                                                         >  2:867187/1‑151 (MQ=255)                                                          aaaaCAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCAt                                                                                      >  1:787151/1‑151 (MQ=255)                                                                  aTGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGttttt                                                                              >  1:444539/1‑151 (MQ=255)                                                                    gCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTgg                                                                            <  2:808789/151‑1 (MQ=255)                                                                      aGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTgggg                                                                          >  1:368110/1‑151 (MQ=255)                                                                      aGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTgggg                                                                          >  2:300773/1‑151 (MQ=255)                                                                       gTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTggggg                                                                         <  1:645433/151‑1 (MQ=255)                                                                         ttttGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGAc                                                                       >  1:691436/1‑151 (MQ=255)                                                                             gCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGa                                                                    <  1:416881/150‑1 (MQ=255)                                                                             gCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGa                                                                    >  2:152395/1‑150 (MQ=255)                                                                             gCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGa                                                                    <  1:884274/150‑1 (MQ=255)                                                                              cTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGaa                                                                   >  2:847821/1‑150 (MQ=255)                                                                               ttttCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAAc                                                                  >  2:930227/1‑150 (MQ=255)                                                                               ttttCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAAc                                                                  >  2:485305/1‑150 (MQ=255)                                                                               ttttCGGTTCTGCCGTACTCCCGATGATAAAAGGCGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAAc                                                                  >  1:294766/1‑150 (MQ=255)                                                                                tttCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACat                                                                >  1:832800/1‑151 (MQ=255)                                                                                  tCGGTTCTGCCGTACTCCCGATGATAAAAGGGGATTTCTCTTTTGCGATGATGCCCTGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGTGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACat                                                                <  1:4383/149‑1 (MQ=255)                                                                                   cGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATc                                                             >  2:454677/1‑151 (MQ=255)                                                                                    ggTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATc                                                             >  2:698063/1‑150 (MQ=255)                                                                                     gTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCaa                                                           <  2:417885/151‑1 (MQ=255)                                                                                        cTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGt                                                         >  1:27037/1‑150 (MQ=255)                                                                                          gCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATc                                                      <  2:800014/151‑1 (MQ=255)                                                                                                cTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAgg                                                >  2:670921/1‑151 (MQ=255)                                                                                                  cccGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAgg                                                <  1:807516/149‑1 (MQ=255)                                                                                                  cccGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGtt                                              <  2:166629/151‑1 (MQ=255)                                                                                                  cccGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGt                                               <  2:622213/150‑1 (MQ=255)                                                                                                    cGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTttc                                            >  1:259022/1‑151 (MQ=255)                                                                                                     gatgatAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTttct                                           >  2:247217/1‑151 (MQ=255)                                                                                                            aaaGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGAc                                     >  1:496303/1‑150 (MQ=255)                                                                                                                gTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCa                                <  2:515385/151‑1 (MQ=255)                                                                                                                      tCGCTTTTGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTAtt                           <  2:129802/150‑1 (MQ=255)                                                                                                                        gCTATTGCGATGATGCGCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGTGCAGGCTGCAGTAGGGTTACTGTTCGGATCGCGGCAGTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTccc                        >  1:238176/1‑151 (MQ=255)                                                                                                                          ttttGCGATGATGCCCGGTGAGTTGGTATTGCTCCTTCAGCGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCCTCAGGTTTCTTCTGCCATTCACTATTCCCtt                      >  1:850017/1‑151 (MQ=255)                                                                                                                          ttttGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCtt                      >  1:25157/1‑151 (MQ=255)                                                                                                                           tttGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCt                       <  2:560953/149‑1 (MQ=255)                                                                                                                            ttGCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCttt                     <  2:923594/150‑1 (MQ=255)                                                                                                                              gCGATGATGCCCGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTc                   <  2:507583/150‑1 (MQ=255)                                                                                                                                        ccGGTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACAt        >  1:242896/1‑151 (MQ=255)                                                                                                                                           gTCAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACAtcct     >  2:357057/1‑151 (MQ=255)                                                                                                                                             cAGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACAtcctcc   <  2:811796/151‑1 (MQ=255)                                                                                                                                              aGTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTccc  >  2:568262/1‑151 (MQ=255)                                                                                                                                               gTTGGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTccc  >  2:182719/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCCAAACTGCTTACCGAGTAATAAGGTTTTTGTTTTCCGTTACACACCCTGGGCATGAAAACAAAATGCAGTTTTTGCTTTTCGGTTCTGCCGTACTCCCGATGATAAAAGGTGATTTCGCTTTTGCGATGATGCCCGGTCAGTTTGTATTGCTCCTTCAGTGCTTCAAGAGCAGGCTGCAGTAGTGTTACTGTTCGGATCCCGGCATTTGTTTTTGGGGGACCGAACATATCAAGTATCGTCAGGTTTCTTCTGACATTCACTATTCCCTTTTCGAGATCCACATCCTCCC  >  NC_000913/1412226‑1412517

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