Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,258,987	T→C	T104A (ACC→GCC)	psuK ←	putative pseudouridine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,258,987	0	T	C	100.0%	251.0 / NA	68	T104A (ACC→GCC)	psuK	putative pseudouridine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/28);  total (40/28)

GTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCA  >  NC_000913/2258852‑2259118                                                                                                                                        |                                                                                                                                    gTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGcatta                                                                                                                                 <  2:974694/140‑1 (MQ=255)  tGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGcattat                                                                                                                                >  1:894410/1‑140 (MQ=255)  tGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGcattat                                                                                                                                >  2:1152779/1‑140 (MQ=255)   gCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATc                                                                                                                               <  2:907232/140‑1 (MQ=255)     ctctTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGa                                                                                                                             <  2:204555/140‑1 (MQ=255)      tctTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGTTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAg                                                                                                                            <  2:505602/140‑1 (MQ=255)      tctTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAg                                                                                                                            <  2:898805/140‑1 (MQ=255)      tctTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAg                                                                                                                            <  1:304384/140‑1 (MQ=255)        ttCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTa                                                                                                                          <  1:3822/140‑1 (MQ=255)         tCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTaa                                                                                                                         >  1:492142/1‑140 (MQ=255)          cACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAAt                                                                                                                        <  2:489051/140‑1 (MQ=255)               atatTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGataa                                                                                                                   >  2:89151/1‑140 (MQ=255)                 atTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGataaat                                                                                                                 >  2:781409/1‑140 (MQ=255)                  ttACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAg                                                                                                                            >  1:1168959/1‑128 (MQ=255)                  ttACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAg                                                                                                                            <  2:1168959/128‑1 (MQ=255)                    aCAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAAc                                                                                                              >  1:699534/1‑140 (MQ=255)                      aGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTc                                                                                                            <  2:440497/140‑1 (MQ=255)                        tCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGa                                                                                                          >  2:304177/1‑140 (MQ=255)                         ccGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATACCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGAc                                                                                                         >  2:158535/1‑140 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTaa                                                                                                                         <  2:1061698/123‑1 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTaa                                                                                                                         <  2:543261/123‑1 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTaa                                                                                                                         >  1:543261/1‑123 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTaa                                                                                                                         >  1:1061698/1‑123 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGAc                                                                                                         <  2:1058114/139‑1 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGAc                                                                                                         >  1:1058114/1‑139 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACg                                                                                                        >  1:1008117/1‑140 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACg                                                                                                        >  2:961613/1‑140 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACg                                                                                                        >  1:728844/1‑140 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACg                                                                                                        >  2:407147/1‑140 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACg                                                                                                        >  2:203627/1‑140 (MQ=255)                          cgcgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACg                                                                                                        >  1:397025/1‑140 (MQ=255)                            cgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTa                                                                                                      >  1:660257/1‑140 (MQ=255)                            cgACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTa                                                                                                      >  2:240321/1‑140 (MQ=255)                             gACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTAt                                                                                                     >  2:31833/1‑140 (MQ=255)                               cAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTAttt                                                                                                   <  1:203627/140‑1 (MQ=255)                               cAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTAttt                                                                                                   <  1:407147/140‑1 (MQ=255)                                  tGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTtct                                                                                                >  2:456076/1‑140 (MQ=255)                                  tGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTtct                                                                                                >  1:197871/1‑140 (MQ=255)                                   gACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTtctc                                                                                               <  2:660257/140‑1 (MQ=255)                                        ttGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGc                                                                                          <  2:699534/140‑1 (MQ=255)                                         tGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGca                                                                                         <  2:599633/140‑1 (MQ=255)                                             ctctGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAAt                                                                                     <  1:781409/140‑1 (MQ=255)                                                tGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAg                                                                                  >  1:843615/1‑140 (MQ=255)                                                        ttCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTAt                                                                          <  1:456076/140‑1 (MQ=255)                                                         tCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATc                                                                         <  2:365736/140‑1 (MQ=255)                                                          cACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAgg                                                                                 >  2:774654/1‑131 (MQ=255)                                                          cACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAgg                                                                                 <  1:774654/131‑1 (MQ=255)                                                          cACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCg                                                                        >  2:742792/1‑140 (MQ=255)                                                            ccGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGAc                                                                      >  1:1144205/1‑140 (MQ=255)                                                                 cTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATaaa                                                                 >  1:23082/1‑140 (MQ=255)                                                                     gCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGcc                                                             <  1:220219/140‑1 (MQ=255)                                                                     gCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAAACAGGCATTTATCGACATAAACGcc                                                             <  2:1008117/140‑1 (MQ=255)                                                                      cGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCa                                                            <  1:620021/140‑1 (MQ=255)                                                                             ttCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAgattgat                                                     <  2:306200/140‑1 (MQ=255)                                                                                  cTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGtt                                                <  1:240321/140‑1 (MQ=255)                                                                                            gTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTagca                                      >  2:681495/1‑140 (MQ=255)                                                                                                cTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGa                                  >  1:544213/1‑140 (MQ=255)                                                                                                  aaTATTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGAtt                                <  1:304177/140‑1 (MQ=255)                                                                                                    tatTCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGa                              <  2:728844/140‑1 (MQ=255)                                                                                                      ttCATGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGAcc                            <  1:31833/140‑1 (MQ=255)                                                                                                          tGTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATaa                        >  1:883985/1‑140 (MQ=255)                                                                                                           gTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATaaa                       >  1:50704/1‑140 (MQ=255)                                                                                                           gTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATaaa                       >  1:1125277/1‑140 (MQ=255)                                                                                                           gTCATTTATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATaaa                       >  2:823612/1‑140 (MQ=255)                                                                                                                ttATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCa                  >  1:797270/1‑140 (MQ=255)                                                                                                                ttATAGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCa                  >  1:726970/1‑140 (MQ=255)                                                                                                                    aGCAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGc              >  2:1128890/1‑140 (MQ=255)                                                                                                                      cAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCt            >  1:967459/1‑140 (MQ=255)                                                                                                                      cAACCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCt            >  2:813577/1‑140 (MQ=255)                                                                                                                        aCCAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTAc          >  1:267183/1‑140 (MQ=255)                                                                                                                         ccAGCATTTCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACg         >  1:420032/1‑140 (MQ=255)                                                                                                                                ttCACCGGCATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCa  >  1:715612/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    GTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCA  >  NC_000913/2258852‑2259118

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