Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,639,617	T→C	G147G (GGA→GGG)	rrrQ ←	Qin prophage; putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,639,617	0	T	C	100.0%	247.3 / NA	68	G147G (GGA→GGG)	rrrQ	Qin prophage; putative lysozyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/43);  total (25/43)

CCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTG  >  NC_000913/1639475‑1639759                                                                                                                                               |                                                                                                                                               ccGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATc                                                                                                                                        >  2:839084/1‑151 (MQ=255) ccGACGATAAATCCTGACAACACGAGAAAAATCACCATAAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATACCCTGACCGTAACATTTTTTTGAACGAATGCGGCAATCGCGCCCGCCAt                                                                                                                                         <  1:935298/150‑1 (MQ=255)   gacgaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCt                                                                                                                                      <  2:100544/151‑1 (MQ=255)      gaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTaa                                                                                                                                   <  2:897557/151‑1 (MQ=255)      gaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTa                                                                                                                                    >  1:1097951/1‑150 (MQ=255)                gACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCg                                                                                                                         <  1:319714/151‑1 (MQ=255)                 acaacaCGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGa                                                                                                                        <  2:602142/151‑1 (MQ=255)                    acacGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATc                                                                                                                     <  1:98176/151‑1 (MQ=255)                     cacGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCg                                                                                                                    <  1:1074071/151‑1 (MQ=255)                     cacGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCg                                                                                                                    <  1:427324/151‑1 (MQ=255)                       cGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCt                                                                                                                  <  2:184373/151‑1 (MQ=255)                         agaAAAATCACCATGAATAGCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCtt                                                                                                                 <  1:796869/150‑1 (MQ=255)                          gaAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCg                                                                                                               <  2:313794/151‑1 (MQ=255)                          gaAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCg                                                                                                               <  1:789964/151‑1 (MQ=255)                              aaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATg                                                                                                           <  1:627939/151‑1 (MQ=255)                                tCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGc                                                                                                          >  1:440664/1‑150 (MQ=255)                                 cACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCAc                                                                                                        >  2:1127690/1‑151 (MQ=255)                                   ccATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCt                                                                                                      >  1:698711/1‑151 (MQ=255)                                               gATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCa                                                                                             <  2:653830/148‑1 (MQ=255)                                               gATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCa                                                                                             <  1:772272/148‑1 (MQ=255)                                                aTTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAcc                                                                                           >  2:200599/1‑149 (MQ=255)                                                aTTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAg                                                                                         <  1:757635/151‑1 (MQ=255)                                                aTTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAg                                                                                         <  1:273256/151‑1 (MQ=255)                                                  tCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGc                                                                                        <  2:711951/150‑1 (MQ=255)                                                   cACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCAt                                                                                      >  1:276076/1‑151 (MQ=255)                                                       gTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAg                                                                                  <  2:745087/151‑1 (MQ=255)                                                       gTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTAACGACGGATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTc                                                                                    >  2:260456/1‑149 (MQ=255)                                                            aTCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCt                                                                             >  2:80482/1‑151 (MQ=255)                                                              cccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCtt                                                                            <  1:828244/150‑1 (MQ=255)                                                              cccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTta                                                                           <  2:776971/151‑1 (MQ=255)                                                               ccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTtat                                                                          >  1:1120074/1‑151 (MQ=255)                                                               ccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTtat                                                                          >  2:908105/1‑151 (MQ=255)                                                               ccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTtat                                                                          >  2:517926/1‑151 (MQ=255)                                                               ccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTtat                                                                          <  2:231112/151‑1 (MQ=255)                                                                   agcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATaaaa                                                                      <  1:437699/151‑1 (MQ=255)                                                                   agcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATaaa                                                                       >  2:1106467/1‑150 (MQ=255)                                                                   agcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATaaa                                                                       >  2:882069/1‑150 (MQ=255)                                                                       ggTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTc                                                                  <  2:980604/151‑1 (MQ=255)                                                                         ttAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCg                                                                 >  2:380692/1‑150 (MQ=255)                                                                          tAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGa                                                                <  2:748393/150‑1 (MQ=255)                                                                                        ggTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTa                                                   <  2:856602/149‑1 (MQ=255)                                                                                           cACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCggg                                              <  1:280079/151‑1 (MQ=255)                                                                                           cACGACGAATAACCTCACCTTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCTCCGTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGTGAAACACTTACCggg                                              <  1:824682/151‑1 (MQ=255)                                                                                              gacgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGcc                                           >  1:525835/1‑151 (MQ=255)                                                                                              gacgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGcc                                           >  1:568362/1‑151 (MQ=255)                                                                                               acgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGcc                                           <  2:750049/150‑1 (MQ=255)                                                                                               acgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCa                                          <  1:734083/151‑1 (MQ=255)                                                                                                 gaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAAt                                        >  2:820076/1‑151 (MQ=255)                                                                                                       ccTGACCGTAACAGTTATTTGAACGAATGCGGCAATTGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTtata                                  <  1:202526/151‑1 (MQ=255)                                                                                                       ccTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTtata                                  <  1:639829/151‑1 (MQ=255)                                                                                                         tGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAgg                                >  2:1046008/1‑151 (MQ=255)                                                                                                          gACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAgg                                <  2:304693/150‑1 (MQ=255)                                                                                                            ccGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGAc                             <  1:94180/151‑1 (MQ=255)                                                                                                            ccGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGAc                             <  1:980778/151‑1 (MQ=255)                                                                                                             cGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGAc                             >  2:582069/1‑150 (MQ=255)                                                                                                             cGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGAc                             >  1:621988/1‑150 (MQ=255)                                                                                                              gTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACaa                           <  1:834707/151‑1 (MQ=255)                                                                                                                   aGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGa                      <  2:1116030/151‑1 (MQ=255)                                                                                                                    gTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAc                     <  2:707546/151‑1 (MQ=255)                                                                                                                       aTTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATCAACCGTCGACGGGAAACCCTTACCGGGGCCAATGTTATAGGGACAAAATGAcgcg                  <  2:146667/151‑1 (MQ=255)                                                                                                                         ttGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGAt                >  2:902405/1‑151 (MQ=255)                                                                                                                         ttGAACGAAGGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGa                 <  2:753261/150‑1 (MQ=255)                                                                                                                               gAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGc          >  1:601018/1‑151 (MQ=255)                                                                                                                                 aTGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCtt        >  2:1177592/1‑151 (MQ=255)                                                                                                                                   gCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCttt       <  2:985405/150‑1 (MQ=255)                                                                                                                                      gCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgt   >  1:1010472/1‑151 (MQ=255)                                                                                                                                      gCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgt   <  1:1127438/151‑1 (MQ=255)                                                                                                                                       cAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgtg  <  1:92001/151‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               CCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTG  >  NC_000913/1639475‑1639759

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