Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	582,152	C→T	pseudogene (653/653 nt)	ybcY ←	pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	582,152	0	C	T	100.0%	327.5 / NA	99	pseudogene (653/653 nt)	ybcY	pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (52/47);  total (52/47)

GTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCT  >  NC_000913/582028‑582282                                                                                                                             |                                                                                                                                   gtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGaa                                                                                                                     >  1:769252/1‑139 (MQ=255)   tGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGGAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGc                                                                                                                   >  1:1097072/1‑139 (MQ=255)        aaCCGTGTTGATACATCAACTGCACCTGATATTGAGTTGCCTACGAACCCTGTCTGGGAGTAATCAGTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGATT‑GATGTTATTTTTTTCCTGAAGCGGaaa                                                                                                              >  1:1568633/1‑139 (MQ=255)          ccGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAAc                                                                                                            <  1:292314/139‑1 (MQ=255)                  aTACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               <  1:1291981/128‑1 (MQ=255)                  aTACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               >  2:1291981/1‑128 (MQ=255)                    aCATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                  <  2:731948/139‑1 (MQ=255)                       tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAc                                                                                               >  2:549488/1‑139 (MQ=255)                              cACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAc                                                                                        >  2:1046930/1‑139 (MQ=255)                               aCCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                       <  1:205359/139‑1 (MQ=255)                                ccTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                      <  2:769252/139‑1 (MQ=255)                                   gATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTgg                                                                                   >  1:652948/1‑139 (MQ=255)                                   gATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTgg                                                                                   >  2:215829/1‑139 (MQ=255)                                      atTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGttt                                                                                >  2:388384/1‑139 (MQ=255)                                         gAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                             <  1:1046930/139‑1 (MQ=255)                                           gTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                             <  2:47474/105‑1 (MQ=255)                                           gTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                             >  1:47474/1‑105 (MQ=255)                                           gTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACa                                                                           >  2:921712/1‑139 (MQ=255)                                             ggCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACAtt                                                                         >  1:85061/1‑139 (MQ=255)                                              gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                        <  1:1457934/107‑1 (MQ=255)                                              gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                        >  2:1457934/1‑107 (MQ=255)                                               ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct                                                                       >  2:8537/1‑139 (MQ=255)                                                   cGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               >  2:115040/1‑95 (MQ=255)                                                   cGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               <  1:115040/95‑1 (MQ=255)                                                   cGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaa                                                                   >  2:1367823/1‑139 (MQ=255)                                                      aCCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATg                                                                <  1:17694/139‑1 (MQ=255)                                                        ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                              <  1:1367823/139‑1 (MQ=255)                                                        ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                              <  1:173691/139‑1 (MQ=255)                                                         cTGTCAGGGAGTAATCATTGGGATTATGCCGCAACACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctc                                                             >  1:53629/1‑139 (MQ=255)                                                          tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                            <  2:1097072/139‑1 (MQ=255)                                                            tCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                          <  2:85061/139‑1 (MQ=255)                                                             cAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTa                                                                              >  2:155854/1‑118 (MQ=255)                                                             cAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTa                                                                              <  1:155854/118‑1 (MQ=255)                                                                ggAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               <  1:1068463/82‑1 (MQ=255)                                                                ggAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               >  2:1068463/1‑82 (MQ=255)                                                                  aGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                            >  1:1611006/1‑131 (MQ=255)                                                                  aGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                            <  2:1611006/131‑1 (MQ=255)                                                                   gTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                   <  1:107780/139‑1 (MQ=255)                                                                      aTCATTGGGATTATGTCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                          >  2:1280985/1‑129 (MQ=255)                                                                      aTCATTGGGATTATGTCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                          <  1:1280985/129‑1 (MQ=255)                                                                      aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                          <  2:1408261/81‑1 (MQ=255)                                                                      aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                          >  1:1408261/1‑81 (MQ=255)                                                                      aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                   >  1:284215/1‑88 (MQ=255)                                                                      aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                   <  2:284215/88‑1 (MQ=255)                                                                      aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATAtgt                                                >  2:540767/1‑139 (MQ=255)                                                                                  aTGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCcttct                                    >  2:925001/1‑139 (MQ=255)                                                                                      cGCAGCACGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCgg                                <  1:921712/139‑1 (MQ=255)                                                                                       gcagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAg                                                                                                                    >  2:564351/1‑54 (MQ=38)                                                                                       gcagcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAg                                                                                                                    <  1:564351/54‑1 (MQ=38)                                                                                         agcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGTGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAAt                             >  1:286912/1‑139 (MQ=255)                                                                                         agcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAAt                             >  1:1132978/1‑139 (MQ=255)                                                                                         agcaCGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAAt                             >  1:189444/1‑139 (MQ=255)                                                                                             cGTCTTAAGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCttt                         >  1:624909/1‑139 (MQ=255)                                                                                                   aaGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               <  1:884103/47‑1 (MQ=38)                                                                                                   aaGCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGaa                                                                                                               >  2:884103/1‑47 (MQ=37)                                                                                                     gCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                      >  2:233498/1‑70 (MQ=255)                                                                                                     gCAAGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                      <  1:233498/70‑1 (MQ=255)                                                                                                        aGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaa              <  1:215829/139‑1 (MQ=255)                                                                                                        aGAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaa              >  2:733546/1‑139 (MQ=255)                                                                                                         gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                   >  2:420895/1‑53 (MQ=38)                                                                                                         gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                   <  1:420895/53‑1 (MQ=38)                                                                                                         gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaa                                                                   >  1:1490715/1‑85 (MQ=255)                                                                                                         gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGaa                                                                   <  2:1490715/85‑1 (MQ=255)                                                                                                         gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCg                                 >  1:407006/1‑119 (MQ=255)                                                                                                         gAACATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCg                                 <  2:407006/119‑1 (MQ=255)                                                                                                          agcATGCTGCGGTT‑GATGTTAGTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                      >  1:968163/3‑65 (MQ=255)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                        <  2:835953/47‑1 (MQ=37)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                        >  1:835953/1‑47 (MQ=38)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                  <  2:215337/53‑1 (MQ=38)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                  >  1:215337/1‑53 (MQ=38)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                       <  2:875709/64‑1 (MQ=255)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                       >  1:875709/1‑64 (MQ=255)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                      <  2:968163/65‑1 (MQ=255)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtt                          <  2:925504/125‑1 (MQ=255)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtt                          >  1:925504/1‑125 (MQ=255)                                                                                                          aaCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAg            >  1:679589/1‑139 (MQ=255)                                                                                                           aCATGCTGCGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGa           <  1:8537/139‑1 (MQ=255)                                                                                                              tgctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                          >  2:956805/1‑89 (MQ=255)                                                                                                              tgctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                          <  1:956805/89‑1 (MQ=255)                                                                                                              tgctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGc        >  2:1070694/1‑139 (MQ=255)                                                                                                               gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                   <  2:797565/47‑1 (MQ=38)                                                                                                               gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACa                                                                                                   >  1:797565/1‑47 (MQ=38)                                                                                                               gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAcc                                                                                              <  2:653874/52‑1 (MQ=38)                                                                                                               gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAcc                                                                                              >  1:653874/1‑52 (MQ=38)                                                                                                               gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                       >  2:130966/1‑59 (MQ=255)                                                                                                               gctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                       <  1:130966/59‑1 (MQ=255)                                                                                                                ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt      >  2:683503/1‑139 (MQ=255)                                                                                                                ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt      <  2:1132978/139‑1 (MQ=255)                                                                                                                ctgcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCt      <  2:189444/139‑1 (MQ=255)                                                                                                                  gcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt    <  2:1416949/139‑1 (MQ=255)                                                                                                                  gcGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt    <  1:388384/139‑1 (MQ=255)                                                                                                                   cGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt    <  1:901278/138‑1 (MQ=255)                                                                                                                   cGGTT‑GATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCGTTTGTGTTTGAAAAGATGCCttt    >  2:901278/1‑138 (MQ=255)                                                                                                                     ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCt  <  1:733546/135‑1 (MQ=255)                                                                                                                     ggttgATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCt  >  1:251392/5‑139 (MQ=255)                                                                                                                         gATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                  <  1:976027/39‑1 (MQ=38)                                                                                                                         gATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAg                                                                                                  >  2:976027/1‑39 (MQ=38)                                                                                                                         gATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGa           <  2:929771/126‑1 (MQ=255)                                                                                                                         gATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTCAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTGTGTGTTTGAAAAGa           >  1:929771/1‑126 (MQ=255)                                                                                                                             |                                                                                                                                   GTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTTCT  >  NC_000913/582028‑582282

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