Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,765,794	(C)7→8	intergenic (‑18/‑124)	yfjM ← / → rnlA	CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,765,787	1	.	C	100.0%	239.6 / NA	69	intergenic (‑11/‑131)	yfjM/rnlA	CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (28/41);  total (28/41)

CTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAA  >  NC_000913/2765665‑2765918                                                                                                                            |                                                                                                                                     cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACg                                                                                                                      >  2:1091017/1‑139 (MQ=255)   gTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGAt                                                                                                                    >  2:103094/1‑139 (MQ=255)    tttCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATc                                                                                                                   <  2:138188/139‑1 (MQ=255)       cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTg                                                                                                                <  2:281766/139‑1 (MQ=255)        gACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGc                                                                                                               <  1:1091017/139‑1 (MQ=255)        gACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGc                                                                                                               <  1:311071/139‑1 (MQ=255)          cGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGccc                                                                                                             >  2:1327550/1‑139 (MQ=255)                        ttGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                               <  2:575972/139‑1 (MQ=255)                         tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                              >  2:988638/1‑139 (MQ=255)                                 ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCa                                                                                     <  2:1448556/139‑1 (MQ=255)                                                aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGccc                                                                                                             <  2:1126624/101‑1 (MQ=255)                                                aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGccc                                                                                                             >  1:1126624/1‑101 (MQ=255)                                                aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAGTTTAg                                                                      >  2:928011/1‑139 (MQ=255)                                                  tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGCCGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGc                                                                     >  1:1195178/1‑139 (MQ=255)                                                  tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCa                                                                    >  1:1156692/1‑139 (MQ=255)                                                  tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCa                                                                    <  1:988638/139‑1 (MQ=255)                                                  tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCa                                                                    >  1:397993/1‑139 (MQ=255)                                                   aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaa                                                                   <  1:98295/139‑1 (MQ=255)                                                    tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaaa                                                                  >  1:1136483/1‑139 (MQ=255)                                                     cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                 >  1:1364556/1‑139 (MQ=255)                                                     cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                 <  2:397993/139‑1 (MQ=255)                                                     cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                 >  2:1212471/1‑139 (MQ=255)                                                     cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACAGCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaaag                                                                 <  1:918341/139‑2 (MQ=255)                                                       aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGccc                                                                                                             <  1:1460752/94‑1 (MQ=255)                                                       aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGcca                                                                                                             >  2:1460752/1‑93 (MQ=255)                                                         ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTGTAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                             <  2:128057/139‑1 (MQ=255)                                                          tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAGTTTAGCAAATCAATa                                                            <  1:928011/139‑1 (MQ=255)                                                          tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                            <  2:1100066/139‑1 (MQ=255)                                                          tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                            <  2:1156692/139‑1 (MQ=255)                                                          tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                            <  1:195968/139‑1 (MQ=255)                                                          tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                            <  1:1220832/139‑1 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGa                                                                                                                        >  1:1035502/1‑78 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGa                                                                                                                        >  2:1054961/1‑78 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGa                                                                                                                        <  2:1035502/78‑1 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                          >  1:742289/1‑139 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                          >  1:1461198/1‑139 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                          >  1:928229/1‑139 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                          >  2:971021/1‑139 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                                                                   <  1:1098443/66‑1 (MQ=255)                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cc                                                                                                                                   >  2:1098443/1‑66 (MQ=255)                                                            cggctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGa                                                                                                                        <  1:1054961/76‑1 (MQ=255)                                                                  ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg                                                                      <  1:1151182/121‑1 (MQ=255)                                                                  ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                    <  1:1399679/139‑1 (MQ=255)                                                                  ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGAAGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg                                                                      >  2:1151182/1‑121 (MQ=255)                                                                  ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccc                                                                                                                              <  2:1109602/65‑1 (MQ=255)                                                                  ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccccc                                                                                                                              >  1:1109602/1‑65 (MQ=255)                                                                       gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                               <  1:1212471/139‑1 (MQ=255)                                                                       gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                               <  2:1136483/139‑1 (MQ=255)                                                                       gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCAGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                               <  2:928229/139‑1 (MQ=255)                                                                          cAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                             <  2:1195178/139‑1 (MQ=255)                                                                                  ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTag                                    >  2:1011462/1‑139 (MQ=255)                                                                                   gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                  <  2:774216/124‑1 (MQ=255)                                                                                   gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                  >  1:774216/1‑124 (MQ=255)                                                                                   gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                                                             >  2:232433/1‑48 (MQ=255)                                                                                   gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑cccccccc                                                                                                                             <  1:232433/49‑2 (MQ=255)                                                                                     tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                 <  1:971021/139‑1 (MQ=255)                                                                                      aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGtt                                >  2:242039/1‑139 (MQ=255)                                                                                      aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGtt                                <  2:1461198/139‑1 (MQ=255)                                                                                       ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccc                                                                                                                                >  2:819560/1‑42 (MQ=255)                                                                                       ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑ccccc                                                                                                                                <  1:819560/42‑1 (MQ=255)                                                                                         gAAACGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                             <  2:1364556/139‑1 (MQ=255)                                                                                             cGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATc                                                                                                                   <  2:612086/50‑1 (MQ=255)                                                                                             cGATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATc                                                                                                                   >  1:612086/1‑50 (MQ=255)                                                                                              gATGGACATCATCCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCata                        <  1:103094/139‑1 (MQ=255)                                                                                                       catcCCCCATTTTGTTGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCGAt               <  1:312325/139‑1 (MQ=255)                                                                                                       catcCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                  <  2:722475/104‑1 (MQ=255)                                                                                                       catcCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                  >  1:722475/1‑104 (MQ=255)                                                                                                        atcCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                 >  1:1120222/1‑120 (MQ=255)                                                                                                        atcCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                 <  2:1120222/120‑1 (MQ=255)                                                                                                        atcCCCCATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg              <  1:1327550/139‑1 (MQ=255)                                                                                                          cccccATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                    <  2:273023/99‑1 (MQ=255)                                                                                                          cccccATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                    >  1:273023/1‑99 (MQ=255)                                                                                                             ccATTTTGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc         <  2:742289/139‑1 (MQ=255)                                                                                                                 tttGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtat                                            <  1:1059420/100‑1 (MQ=255)                                                                                                                 tttGTGGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGtaa                                            >  2:1059420/1‑99 (MQ=255)                                                                                                                   tgtgGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                            >  2:1457353/1‑82 (MQ=255)                                                                                                                   tgtgGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                            <  1:1457353/82‑1 (MQ=255)                                                                                                                    gtgGACTGC‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaaa  <  2:1323892/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                     CTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑‑CCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAA  >  NC_000913/2765665‑2765918

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