Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,188	A→G	intergenic (+55/+502)	yjhV → / ← fecE	pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,188	0	A	G	100.0%	293.2 / NA	81	intergenic (+55/+502)	yjhV/fecE	pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (44/37);  total (44/37)

GCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGA  >  NC_000913/4510057‑4510321                                                                                                                                    |                                                                                                                                      gCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTa                                                                                                                                <  1:340256/139‑1 (MQ=255)               ggTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGt                                                                                                                  >  1:733642/1‑139 (MQ=255)                               aTGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGc                                                                                                  >  1:619012/1‑139 (MQ=255)                                    gATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTc                                                                                             <  2:733642/139‑1 (MQ=255)                                       atCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACTCACGCTTTTTAGAAAACTCGATAGTCAGAAACTGGATTAGCAATTCCGt                                                                                          >  1:1332829/1‑139 (MQ=255)                                       atCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGt                                                                                          >  2:63104/1‑139 (MQ=255)                                                aaCAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaaa                                                                                                            <  1:47014/112‑1 (MQ=255)                                                aaCAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaaa                                                                                                            >  2:47014/1‑112 (MQ=255)                                                   aaacgttcgaTCAACTGGCCAGTATAACGTGGCTGCTTGTGAATAAGAAAGCTGCACCTACAGGGGGAGAAGCAAACGCACGAGTTTTAGAACAATCGATGGACAGAATCTGGATTAGCAATTCCGTTCCATTTTTGCt                                                                              <  2:1332829/130‑1 (MQ=255)                                                   aaTCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                              <  2:660473/139‑1 (MQ=255)                                                     tCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        >  2:1205427/1‑111 (MQ=255)                                                     tCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        <  1:1205427/111‑1 (MQ=255)                                                        aTCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGa                                                                         >  2:812705/1‑139 (MQ=255)                                                                aaagacaaggattaCGTTGCCGCCGGTTATTAAAAAAGCTGCCCCCTCCTAAGTAGTAGCCAACGCACGCTTTTTAGAAAAAATGATGGTTAGAAACTGGATTAGCAATTTCGTTTCATGGTTGCTTTTGATTTACGtt                                                                 <  2:25632/126‑1 (MQ=255)                                                                      cTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAAc                                                                                                           <  1:633116/91‑1 (MQ=255)                                                                      cTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAAc                                                                                                           >  2:633116/1‑91 (MQ=255)                                                                       tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                          >  2:410617/1‑139 (MQ=255)                                                                       tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                          <  1:63104/139‑1 (MQ=255)                                                                         tataACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGa                                                        >  2:1030914/1‑139 (MQ=255)                                                                              cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        <  1:600963/86‑1 (MQ=255)                                                                              cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        >  2:600963/1‑86 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGa                                                                                                       >  1:262867/1‑82 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGa                                                                                                       <  2:262867/82‑1 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGTTTTTGATTTACTTTGGCGTCTGATcctgtcttt                                              >  1:25632/1‑131 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTc                                                           >  1:687362/1‑126 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTc                                                           <  2:687362/126‑1 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                          <  1:981141/127‑1 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                          >  2:981141/1‑127 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                              >  2:176390/1‑139 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                              >  2:1005544/1‑139 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                              >  1:1212227/1‑139 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                <  1:144221/137‑1 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                >  2:144221/1‑137 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                >  2:799670/1‑137 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                <  1:799670/137‑1 (MQ=255)                                                                                        ggTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                             >  1:682526/1‑71 (MQ=255)                                                                                        ggTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                             <  2:682526/71‑1 (MQ=255)                                                                                        ggTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCt                                         >  2:409248/1‑139 (MQ=255)                                                                                          tGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                       <  1:409248/139‑1 (MQ=255)                                                                                            aTTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCttt                                                                            <  1:5469/100‑1 (MQ=255)                                                                                            aTTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCttt                                                                            >  2:5469/1‑100 (MQ=255)                                                                                                 gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                                <  1:176390/139‑1 (MQ=255)                                                                                                 gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                                <  1:410617/139‑1 (MQ=255)                                                                                                    aGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                             <  1:1005544/139‑1 (MQ=255)                                                                                                      cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                             >  2:193017/1‑137 (MQ=255)                                                                                                      cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                             <  1:193017/137‑1 (MQ=255)                                                                                                              acctaAGTAGTAGCAAACGCACGCTTTTTAGAATAATCGATGGTCAGAAACTgg                                                                                                        <  2:792552/54‑1 (MQ=255)                                                                                                              acctaAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        >  1:792552/1‑54 (MQ=255)                                                                                                                   agtagtAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                       <  1:761060/114‑1 (MQ=255)                                                                                                                   agtagtAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                       >  2:761060/1‑114 (MQ=255)                                                                                                                   agtagtAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAtt                       >  2:1323542/1‑130 (MQ=255)                                                                                                                   agtagtAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAtt                       <  1:1323542/130‑1 (MQ=255)                                                                                                                     tagtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                             >  1:384890/1‑42 (MQ=255)                                                                                                                     tagtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                             <  2:384890/42‑1 (MQ=255)                                                                                                                      agtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                             >  2:587637/1‑41 (MQ=255)                                                                                                                      agtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                             <  1:587637/41‑1 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCaa                                                                                                <  2:769167/52‑1 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCaa                                                                                                >  1:769167/1‑52 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTaa                  >  1:106698/1‑130 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTaa                  <  2:106698/130‑1 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGc         >  2:1058133/1‑139 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        <  1:1116966/43‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                        >  2:1116966/1‑43 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                              <  1:225106/53‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                              >  2:225106/1‑53 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                          >  2:68208/1‑89 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                          <  1:68208/89‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                >  2:211139/1‑99 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                                <  1:211139/99‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCt        >  1:359849/1‑139 (MQ=255)                                                                                                                          gCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                                <  2:575579/114‑1 (MQ=255)                                                                                                                          gCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                                >  1:575579/1‑114 (MQ=255)                                                                                                                           cAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTcc      >  1:1126955/1‑139 (MQ=255)                                                                                                                             aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                                   >  2:1192457/1‑44 (MQ=255)                                                                                                                             aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                                   <  1:1192457/44‑1 (MQ=255)                                                                                                                             aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGt    >  2:576227/1‑139 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                                >  2:765224/1‑109 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                                <  1:765224/109‑1 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGa  >  2:209333/1‑139 (MQ=255)                                                                                                                                gCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGt                                                                                  >  1:367739/1‑58 (MQ=255)                                                                                                                                gCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGt                                                                                  <  2:367739/58‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      GCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGA  >  NC_000913/4510057‑4510321

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