Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	540,515	C→T	intergenic (+7/‑50)	allB → / → ybbY	allantoinase/putative uracil/xanthine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	540,515	0	C	T	100.0%	273.8 / NA	83	intergenic (+7/‑50)	allB/ybbY	allantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/43);  total (40/43)

TGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT  >  NC_000913/540388‑540647                                                                                                                                 |                                                                                                                                     tGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                           <  1:536666/139‑1 (MQ=255) tGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                           <  2:756418/139‑1 (MQ=255)     cGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTc                                                                                                                       >  1:574006/1‑139 (MQ=255)       tACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCt                                                                                                                     >  1:1685920/1‑139 (MQ=255)        aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCtt                                                                                                                    >  1:1007952/1‑139 (MQ=255)        aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCtt                                                                                                                    <  2:1007952/139‑1 (MQ=255)        aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCtt                                                                                                                    >  2:610174/1‑139 (MQ=255)        aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCtt                                                                                                                    >  1:786170/1‑139 (MQ=255)            ttGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgcgg                                                                                                                <  2:1334092/139‑1 (MQ=255)             tGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgcggc                                                                                                               <  1:914260/139‑1 (MQ=255)                     gTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtc                                                                                                       <  2:1178138/139‑1 (MQ=255)                           cGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGt                                                                                                 <  2:777055/139‑1 (MQ=255)                            gAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                <  1:610174/139‑1 (MQ=255)                             aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                               <  2:574006/139‑1 (MQ=255)                             aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                               <  2:503768/139‑1 (MQ=255)                             aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                               <  2:786170/139‑1 (MQ=255)                                     cTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtgt                                                                                       <  2:1674721/139‑1 (MQ=255)                                            ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                >  1:772889/1‑123 (MQ=255)                                            ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                <  1:1020395/123‑1 (MQ=255)                                            ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                >  2:1020395/1‑123 (MQ=255)                                            ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                <  2:772889/123‑1 (MQ=255)                                            ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGtt                                                                                >  1:538530/1‑139 (MQ=255)                                                  gtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTcc                                                                                                                      <  1:1047649/95‑1 (MQ=255)                                                  gtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTcc                                                                                                                      >  2:1047649/1‑95 (MQ=255)                                                    gatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTg                                                                        >  1:435260/1‑139 (MQ=255)                                                       ttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAg                                                                     >  2:1491848/1‑139 (MQ=255)                                                                tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                            >  1:879275/1‑139 (MQ=255)                                                                       ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                           >  2:1567781/1‑102 (MQ=255)                                                                       ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                           <  1:1567781/102‑1 (MQ=255)                                                                       ggCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                                 >  1:263708/1‑127 (MQ=255)                                                                       ggCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                                 <  2:263708/127‑1 (MQ=255)                                                                         cTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgcgg                                                                                                                <  2:1734738/78‑1 (MQ=255)                                                                         cTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgcgg                                                                                                                >  1:1734738/1‑78 (MQ=255)                                                                            cccTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAAt                                                                                                                              <  2:898080/61‑1 (MQ=255)                                                                            cccTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAAt                                                                                                                              >  1:898080/1‑61 (MQ=255)                                                                            cccTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                           <  2:683129/64‑1 (MQ=255)                                                                            cccTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                           >  1:683129/1‑64 (MQ=255)                                                                            cccTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                >  2:305873/1‑139 (MQ=255)                                                                            cccTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                <  2:498069/139‑1 (MQ=255)                                                                                 ttGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa                                           <  2:1685920/139‑1 (MQ=255)                                                                                      ccGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCa                                                                                                                                >  1:1182442/1‑49 (MQ=255)                                                                                      ccGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCa                                                                                                                                <  2:1182442/49‑1 (MQ=255)                                                                                      ccGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                      >  2:215108/1‑139 (MQ=255)                                                                                      ccGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                      >  1:787090/1‑139 (MQ=255)                                                                                      ccGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                      >  1:1496994/1‑139 (MQ=255)                                                                                      ccGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                                      >  1:360429/1‑139 (MQ=255)                                                                                         aaaGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                           <  2:1421382/67‑1 (MQ=255)                                                                                         aaaGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                           >  1:1421382/1‑67 (MQ=255)                                                                                         aaaGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                                 <  2:904896/109‑1 (MQ=255)                                                                                         aaaGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                                 >  1:904896/1‑109 (MQ=255)                                                                                                        ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                <  1:704872/64‑1 (MQ=255)                                                                                                        ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                                >  2:704872/1‑64 (MQ=255)                                                                                                            aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                                     <  1:201/55‑1 (MQ=255)                                                                                                            aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                                     >  2:201/1‑55 (MQ=255)                                                                                                             aaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                               <  1:174725/60‑1 (MQ=255)                                                                                                             aaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                               >  2:174725/1‑60 (MQ=255)                                                                                                                aTCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAac             <  2:1496994/139‑1 (MQ=255)                                                                                                                 tCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca            >  2:1453893/1‑139 (MQ=255)                                                                                                                 tCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAaca            <  1:215108/139‑1 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                                     >  2:10228/1‑49 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                                     <  1:10228/49‑1 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTaa                                                                                              >  2:357225/1‑56 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTaa                                                                                              <  1:357225/56‑1 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                           >  2:772514/1‑59 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                           <  1:772514/59‑1 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                <  2:1317896/102‑1 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                >  1:1317896/1‑102 (MQ=255)                                                                                                                  cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac           <  2:435260/139‑1 (MQ=255)                                                                                                                        tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGaa                                                          >  1:347102/1‑86 (MQ=255)                                                                                                                        tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGaa                                                          <  2:347102/86‑1 (MQ=255)                                                                                                                        tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa              >  1:1475991/1‑130 (MQ=255)                                                                                                                        tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac           <  2:312332/133‑1 (MQ=255)                                                                                                                        tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac           >  1:312332/1‑133 (MQ=255)                                                                                                                        tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGACTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCaa              <  2:1475991/130‑1 (MQ=255)                                                                                                                          aTCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                            <  1:603570/82‑1 (MQ=255)                                                                                                                          aTCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                            >  2:603570/1‑82 (MQ=255)                                                                                                                          aTCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                       >  2:1466112/1‑87 (MQ=255)                                                                                                                          aTCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                       <  1:1466112/87‑1 (MQ=255)                                                                                                                           tCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt  <  2:879275/139‑1 (MQ=255)                                                                                                                            cTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                >  1:939512/1‑92 (MQ=255)                                                                                                                            cTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                <  2:939512/92‑1 (MQ=255)                                                                                                                             tGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAgg                                                                                            <  1:921174/47‑1 (MQ=255)                                                                                                                             tGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAgg                                                                                            >  2:921174/1‑47 (MQ=255)                                                                                                                                 |                                                                                                                                     TGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTAT‑CCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT  >  NC_000913/540388‑540647

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