Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,958,261	T→C	intergenic (‑830/‑31)	ppiC ← / → rep	peptidyl‑prolyl cis‑trans isomerase C/rep helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,958,261	0	T	C	100.0%	95.2 / NA	31	intergenic (‑830/‑31)	ppiC/rep	peptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/12);  total (19/12)

TTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑CCT‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGG  >  NC_000913/3958171‑3958349                                                                                                 |                                                                                             ttCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑t‑‑‑ccc‑‑‑cc                                                                                         <  2:1982896/93‑1 (MQ=255)  tCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑ccc‑‑‑ccc                                                                                        >  2:6583680/1‑93 (MQ=255)  tCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑ccc‑‑‑ccc                                                                                        <  2:4844878/93‑1 (MQ=255)   cGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑cccc                                                                                       <  1:1994115/93‑1 (MQ=255)   cGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑cccc                                                                                       <  1:6140179/93‑1 (MQ=255)     gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACCA‑‑‑T‑‑‑Ccc‑‑‑cccccc                                                                                     >  1:3809181/1‑93 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    >  1:5808894/1‑92 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    <  1:136865/93‑2 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    >  1:1444669/1‑92 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    >  2:3261684/1‑92 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    >  1:4164353/1‑92 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    >  2:2762691/1‑92 (MQ=255)      tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑CTC‑‑‑ccccccc                                                                                    >  2:2727152/1‑92 (MQ=255)      gcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccc                                                                                    <  2:5250140/92‑2 (MQ=255)       cttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑cccccccc                                                                                   >  1:5627417/1‑91 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTCCCCCCGGTAACATTCAAATGCGATTCTCCTCCAT‑‑‑T‑‑‑Ccc‑‑‑cccccccct                                                                                  >  2:422596/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  2:4208888/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  2:2329218/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  2:5511329/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  1:3209839/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  2:3037700/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  1:2737248/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  2:2092045/1‑90 (MQ=255)        ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccc                                                                                  >  1:1111818/1‑90 (MQ=255)         tctctctcTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccgt                                                                                 >  1:431849/1‑89 (MQ=255)         tctctctcTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑Ccc‑‑‑ccccccccct                                                                                 >  2:553426/1‑89 (MQ=255)                 tGATTTAAACATCCGCAGCCAACGGTTTAGCGGCTTACACGCGGTCCCATTAAAATGCGATTCTGCTCCAA‑‑‑TC‑‑CCC‑‑‑CCCCCCGTTCGAAGAtt                                                                          <  2:3336110/93‑1 (MQ=255)                 tGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGAt                                                                           >  2:992841/1‑93 (MQ=255)                  gATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TTC‑CCC‑‑‑CCCCCCCTTCGAGGAtt                                                                          >  2:4490928/1‑93 (MQ=255)                   aTTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTg                                                                         <  1:1882940/93‑1 (MQ=255)                    tttAAACATCCGCCGCCAACCGGTTAGCGGGTTTCCCGCGGTCACATTCAAATGCGATTCTGCTACAAAC‑C‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGa                                                                        <  2:5681889/93‑1 (MQ=255)                    tttAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGa                                                                        <  1:2077512/93‑1 (MQ=255)                    tttAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGGCACATTCAAATGCGATTTTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTg                                                                         <  1:3325343/93‑1 (MQ=255)                     ttAAACATCCCCAGCCAACCGGTTAGCGCCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGa                                                                        <  2:154431/93‑1 (MQ=255)                     ttAAACACCCCCAGCCAACCGGGTAAGGGGTTACCCGGGGTCACATTCAAATGCGATTCTGCTACAAACCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGa                                                                        <  2:4352522/93‑1 (MQ=255)                     ttAAAAATCCCCCGCCAACCGGGTAGCGGCTTTCCCGCGGGCCCATTCAAAAGCGATTCTTCTTCAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGa                                                                        <  2:1054358/93‑1 (MQ=255)                       taaCATCCGCCGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAAAGCGATTCTTCTACAAACCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGc                                                                      <  2:6550670/92‑1 (MQ=255)                       taaCATCCGCCGCCAACCCGTTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCa                                                                     <  2:496389/92‑1 (MQ=255)                        aaCATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCaa                                                                    <  1:2570396/93‑1 (MQ=255)                         aCATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TC‑‑CCC‑‑‑CCCCCCGTTCGAAGATTTAGCAATa                                                                  >  2:5964186/1‑93 (MQ=255)                              cGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAAAATG‑GCAATa‑ccacc                                                            >  1:4200921/1‑89 (MQ=255)                                cAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTa                                                           <  2:2153045/93‑1 (MQ=255)                                 cgCCAACCGGTTAGCGGCTTTCACGCGGGTCCATTCAAATGCGATTCTTCCACAAACCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTa                                                           <  2:1023570/92‑1 (MQ=255)                                  gCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAAGTTGACCAATA‑CACCTag                                                          >  2:3296991/1‑92 (MQ=255)                                  gCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCCTTCGAAGATTAGGCAATA‑CACCTATg                                                         >  2:704801/1‑93 (MQ=255)                                   ccAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGc                                                        <  1:3837574/93‑1 (MQ=255)                                    cAACCCGGTAGCGGCTTACACGCGGGCACATTCAAAAGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGc                                                        <  2:321940/93‑1 (MQ=255)                                     aaCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TC‑‑CCC‑‑‑CCCCCCCTTTCAAGATTTAGTAAAA‑CACCTATGCgct                                                     >  1:1317933/1‑91 (MQ=255)                                      aCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAAC‑‑C‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCt                                                    <  2:3841564/93‑1 (MQ=255)                                      aCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTCCGAAATTTGGGAAATA‑CACCTACGCGt                                                      >  2:3302460/1‑93 (MQ=255)                                      aCCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCCTTCGAAGATTGAGCAATA‑CACCTATGCGTc                                                     >  1:3386066/1‑93 (MQ=255)                                        cGGTTAGCGGCTTTCACGCGGTCACATTCAAATGCGGTTTTTCTTCAAACCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCt                                                    <  2:5768874/93‑1 (MQ=255)                                        cGGGTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCt                                                    <  2:4039878/93‑1 (MQ=255)                                         ggTTAGCGGCTTACACCCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTa                                                   <  2:1473519/93‑1 (MQ=255)                                           ttAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCCTTCGAAGATTGAGCAAAC‑CACCTATGGGTCTaaa                                                 >  2:4494312/1‑93 (MQ=255)                                             aGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TC‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAAcccc                                             <  1:331238/93‑1 (MQ=255)                                              ggggCTTACCCCCGGGCACATTCAAATGCGATTCTGCTACAACCCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAAccc                                              <  1:1939322/91‑1 (MQ=255)                                              gCGGCTTACACGCGGTCACATTCAAATGCCATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAAcccc                                             <  2:5858084/93‑1 (MQ=255)                                                 gCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGCGCAATA‑CACCTATACGTCTAAACCCCgc                                           >  2:1081865/1‑92 (MQ=255)                                                  cTTACCCGCGGTTCCATTCAAAAGCGAATCTGCCACAAACCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGc                                          <  2:624763/93‑1 (MQ=255)                                                   ttACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGcc                                         <  2:397568/93‑1 (MQ=255)                                                      cccgcgGTCCCATTCAAATGCGATTCTGCCTCAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCcaac                                      <  1:5915197/91‑1 (MQ=255)                                                        cggcgTCACATTCAAATGCGATTCTGCTGCAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCcaacaac                                   <  1:1689377/89‑1 (MQ=255)                                                         ccgGTCACATTCAAACGCCATTCTGCTACAACC‑C‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCcaacaaca                                  <  1:4217637/92‑1 (MQ=255)                                                          cgGTCACATTCAAATGCGATTCTGCTACAACC‑C‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCcaacaacaa                                 <  1:424677/93‑1 (MQ=255)                                                            gTCAAATTCAAATGCGATTCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGc                               <  2:3740706/93‑1 (MQ=255)                                                               ccaTTCAAATGCGATTCTCCTCCAACC‑C‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGt                            <  2:5324375/92‑1 (MQ=255)                                                                caTTCAAATGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCCTTCAAAGACTG‑GCAACAGAACCTTAGCGTGAAAACCCCGGTCAACAACAAGCTgc                            >  1:4403166/1‑92 (MQ=255)                                                                 aTTCAAATGCGA‑TCTGCCAACA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGa                         <  2:3431425/93‑1 (MQ=255)                                                                   tCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑CCCCCCCCCCCCCTTCTAAGATTTAACAAAA‑CACCCATCCGGCAAAACCCCCGCCACAAACAAGACGTCGa                         >  1:3215905/1‑93 (MQ=255)                                                                   tCAAATGCGATTCTGCTAAAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGa                         <  1:3182987/93‑1 (MQ=255)                                                                    cAAATGCGATTCTGCTACAA‑‑‑TC‑‑CCC‑‑‑CCCCCCCGTCGAAGGTTTAGCAAAA‑AACCCATTCCGCTTAAACCCGGCCCACAACCAACTGCCCAAtt                      >  2:1833882/1‑93 (MQ=255)                                                                     aaaTGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAAATTGAGCAATA‑CACCTATCCGTCTAAACCCCCGCCAACAACAAACTGTAGAAt                       >  2:2265810/1‑93 (MQ=255)                                                                      aaTGCGATTCTGCTACAA‑‑‑TCCCCCC‑‑‑CCCCCCGTTCGAAAATTGAGCAAAA‑CACATTTCCGCCTAAACCCCGGCCAAAAAAAAGCTTTCAAAtt                      >  2:1088564/1‑93 (MQ=255)                                                                       aTGCGATTCTGCTACAACCCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTc                     <  2:3882584/93‑1 (MQ=255)                                                                        agCGATTCTGCTACAACCCC‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCg                    <  2:6015265/92‑1 (MQ=255)                                                                           gATTCTGCCACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAGCCCCGGCCAACAACAAGCTGTCGAATTCGTTAc                <  2:2294406/93‑1 (MQ=255)                                                                             ttCTGCTACAA‑‑‑TCC‑CCC‑‑‑CCCCCCGTTCGAAAATTGGGCAATA‑ACCCTATGGGTTTAAACCCCGGCCAAAAAAAACCTGTCGAATTATTTACCg              >  2:4387677/1‑93 (MQ=255)                                                                             ttCTACTACAA‑‑‑T‑‑‑CCCCC‑CCCCCCGTTCGAAAATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAACCTGTCGAATTCGTTACCg              >  1:784278/1‑93 (MQ=255)                                                                                 gCTACAA‑‑‑T‑‑‑CCCC‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCAAACAACAAGCTGTCGAATTCGTTCCCGGcccc         >  2:2272228/1‑93 (MQ=255)                                                                                       a‑‑‑T‑‑‑CCC‑‑‑CCCCCCGTTCGAAGATTTAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAACTGTCGAATTCGTTACCGGCCCCTGCCTgg  >  1:13974/1‑93 (MQ=255)                                                                                                 |                                                                                             TTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑‑CCT‑‑‑CCCCCCGTTCGAAGATTGAGCAATA‑CACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGG  >  NC_000913/3958171‑3958349

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