Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,556,769	(T)9→8	coding (116/144 nt)	ytiD →	protein YtiD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,556,761	0	T	.	97.1%	281.7 / 2.4	70	coding (108/144 nt)	ytiD	protein YtiD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (33/35);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAT  >  NC_000913/4556627‑4556894                                                                                                                                       |                                                                                                                                      ttACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGtttttt                                                                                                                                  >  2:744190/1‑140 (MQ=255)  tACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGttttttt                                                                                                                                 >  1:817228/1‑140 (MQ=255)           aTTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtga                                                                                                                       <  1:495448/140‑1 (MQ=255)           aTTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTtga                                                                                                                       <  1:359575/140‑1 (MQ=255)                 tCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGAttt                                                                                                                 <  2:596500/140‑1 (MQ=255)                    aataatACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTa                                                                                                              >  2:594723/1‑140 (MQ=255)                       aataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCa                                                                                                           >  1:335763/1‑140 (MQ=255)                        ataCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCaa                                                                                                          <  1:257888/140‑1 (MQ=255)                         taCATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCaaa                                                                                                         >  2:290263/1‑140 (MQ=255)                               gCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGa                                                                                                   >  1:34298/1‑140 (MQ=255)                                cTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGAc                                                                                                  >  1:828736/1‑140 (MQ=255)                                cTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGAc                                                                                                  >  1:926948/1‑140 (MQ=255)                                   aTTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                               >  1:593747/1‑140 (MQ=255)                                   aTTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                               >  1:80984/1‑140 (MQ=255)                                   aTTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                               >  2:741971/1‑140 (MQ=255)                                   aTTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                               >  1:987751/1‑140 (MQ=255)                                   aTTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACttt                                                                                               >  2:10310/1‑140 (MQ=255)                                         cGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCt                                                                                         >  2:240038/1‑140 (MQ=255)                                            aaataaatAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaa                                                                                      <  1:259622/140‑1 (MQ=255)                                             aataaataATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaa                                                                                     >  2:888739/1‑140 (MQ=255)                                              ataaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaaa                                                                                    >  1:157598/1‑140 (MQ=255)                                              ataaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGaaaa                                                                                    >  2:545399/1‑140 (MQ=255)                                                 aataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                 <  2:837164/140‑1 (MQ=255)                                                 aataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                 <  1:888739/140‑1 (MQ=255)                                                 aataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTc                                                                                 <  1:156670/140‑1 (MQ=255)                                                   taatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAg                                                                               <  2:157598/140‑1 (MQ=255)                                                   taatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAg                                                                               <  1:585779/140‑1 (MQ=255)                                                   taatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAg                                                                               <  2:34298/140‑1 (MQ=255)                                                     atTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTa                                                                             <  2:593747/140‑1 (MQ=255)                                                       ttCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGa                                                                           <  1:744190/140‑1 (MQ=255)                                                       ttCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGa                                                                           <  1:129740/140‑1 (MQ=255)                                                            aCCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTc                                                                      >  1:371428/1‑140 (MQ=255)                                                             ccGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCa                                                                     <  1:694740/140‑1 (MQ=255)                                                                gAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAg                                                                  <  1:594723/140‑1 (MQ=255)                                                                       gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTa                                                           <  2:828736/140‑1 (MQ=255)                                                                       gTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTa                                                           <  2:918880/140‑1 (MQ=255)                                                                           cagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGt                                                       >  2:490746/1‑140 (MQ=255)                                                                            agcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGtt                                                      <  1:972701/140‑1 (MQ=255)                                                                                gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                  >  2:598279/1‑140 (MQ=255)                                                                                gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                  <  2:277774/140‑1 (MQ=255)                                                                                gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTtat                                                  <  2:926948/140‑1 (MQ=255)                                                                                     aaGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACttt                                             <  2:817228/140‑1 (MQ=255)                                                                                      aGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTg                                            <  1:490746/140‑1 (MQ=255)                                                                                        tATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCt                                          >  1:124662/1‑140 (MQ=255)                                                                                           tttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGAc                                       <  1:262899/140‑1 (MQ=255)                                                                                            ttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACa                                      <  1:257477/140‑1 (MQ=255)                                                                                               aTCTACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGa                                   >  2:273440/1‑140 (MQ=255)                                                                                                  tACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTc                                >  1:520739/1‑140 (MQ=255)                                                                                                  tACGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTc                                >  2:728656/1‑140 (MQ=255)                                                                                                    cGTAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAg                              >  2:1016887/1‑140 (MQ=255)                                                                                                      tAACTATTTGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGc                            <  2:1055809/140‑1 (MQ=255)                                                                                                      tAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGc                            <  1:240038/140‑1 (MQ=255)                                                                                                      tAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGc                            <  2:80984/140‑1 (MQ=255)                                                                                                      tAACTATTCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGc                            <  1:290263/140‑1 (MQ=255)                                                                                                            ttCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtct                      >  1:897443/1‑140 (MQ=255)                                                                                                            ttCGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtct                      >  1:276238/1‑140 (MQ=255)                                                                                                              cGGAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCa                    >  1:169556/1‑140 (MQ=255)                                                                                                                gAATATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGa                  <  2:152571/140‑1 (MQ=255)                                                                                                                 aaTATTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGAc                 >  2:623820/1‑140 (MQ=255)                                                                                                                   tatTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTg               >  2:656395/1‑140 (MQ=255)                                                                                                                    atTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGa              <  2:897443/140‑1 (MQ=255)                                                                                                                    atTCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCGGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGa              >  2:1008462/1‑140 (MQ=255)                                                                                                                     ttCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGAc             <  1:1008462/140‑1 (MQ=255)                                                                                                                     ttCAGAATATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGAc             <  1:46409/140‑1 (MQ=255)                                                                                                                          aaTATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtc                       <  2:53935/126‑1 (MQ=255)                                                                                                                          aaTATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGtc                       >  1:53935/1‑126 (MQ=255)                                                                                                                          aaTATACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTAAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGa        >  2:543337/1‑140 (MQ=255)                                                                                                                           atatACTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGAt       >  2:55256/1‑140 (MQ=255)                                                                                                                               aCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTaa   >  1:426732/1‑140 (MQ=255)                                                                                                                               aCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTaa   >  2:849300/1‑140 (MQ=255)                                                                                                                               aCTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTaa   <  2:253624/140‑1 (MQ=255)                                                                                                                                cTAGCGG‑TTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAt  <  1:560544/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TTACAAACTAATTGTTTCAAATAATACATGGCTGATTATGCGGAAATAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAAT  >  NC_000913/4556627‑4556894

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