Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,421,120	Δ1 bp	noncoding (60/121 nt)	rrfF ←	5S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,421,119	0	G	.	95.2%	200.6 / 6.8	83	noncoding (61/121 nt)	rrfF	5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base . (47/32);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA  >  NC_000913/3421027‑3421211                                                                                             |                                                                                             cAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcg                                                                                              >  2:4225268/1‑93 (MQ=255)  aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgt                                                                                             <  1:1082308/93‑2 (MQ=39)   gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                            >  2:935234/1‑91 (MQ=255)   gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                            >  1:3249847/1‑91 (MQ=255)   gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                            >  1:2925361/1‑91 (MQ=255)   gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAcggcgtt                                                                                            >  2:4635681/1‑91 (MQ=255)   gTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGt                                                                                            >  1:2076162/1‑93 (MQ=255)     cTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttc                                                                                          >  2:4364973/1‑92 (MQ=255)      ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtttcc                                                                                         >  1:521871/1‑91 (MQ=255)       tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCAc                                                                                       <  2:1233560/93‑1 (MQ=16)       tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCAc                                                                                       >  1:2338624/1‑93 (MQ=255)       tCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCAc                                                                                       <  1:3118767/93‑1 (MQ=16)        ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACt                                                                                      >  2:1862972/1‑93 (MQ=12)         cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                     >  1:440154/1‑93 (MQ=16)         cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                     >  2:4015559/1‑93 (MQ=12)         cGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACtt                                                                                     <  1:2028329/93‑1 (MQ=9)          gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTTCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTACACTTc                                                                                    >  1:1640203/1‑93 (MQ=11)          gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                    <  1:1833238/93‑1 (MQ=9)          gACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTc                                                                                    <  1:3625411/93‑1 (MQ=9)           aCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCt                                                                                   <  1:2969132/93‑1 (MQ=9)            cTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTg                                                                                  >  1:3892180/1‑93 (MQ=16)              gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                                >  1:688226/1‑93 (MQ=12)              gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                                <  2:1344658/93‑1 (MQ=9)              gAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAg                                                                                >  2:1699323/1‑93 (MQ=16)               aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                               >  1:4592875/1‑93 (MQ=16)               aGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGt                                                                               >  1:3669084/1‑93 (MQ=12)                gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                              >  1:3315097/1‑93 (MQ=16)                gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                              >  1:130024/1‑93 (MQ=16)                gCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGAGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                              >  2:680243/1‑93 (MQ=255)                gCCTTTTGGTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGtt                                                                              >  2:4580009/1‑93 (MQ=9)                 ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGttg                                                                             >  1:4089600/1‑92 (MQ=9)                 ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                             >  2:305529/1‑93 (MQ=12)                 ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                             >  2:2363063/1‑93 (MQ=9)                 ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                             <  2:4306236/93‑1 (MQ=9)                 ccTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTc                                                                             >  2:1909026/1‑93 (MQ=16)                  cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCg                                                                            >  2:36382/1‑93 (MQ=9)                  cTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGGTCg                                                                            >  2:1923662/1‑93 (MQ=9)                   ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                           >  1:602791/1‑93 (MQ=9)                   ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                           >  2:1011801/1‑93 (MQ=12)                   ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                           >  2:3058136/1‑93 (MQ=12)                   ttttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCgg                                                                           >  2:4310842/1‑93 (MQ=12)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          <  1:1941508/93‑1 (MQ=9)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          >  2:3224342/1‑93 (MQ=12)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          >  2:2028524/1‑93 (MQ=12)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          >  1:4452682/1‑93 (MQ=9)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          >  1:945575/1‑93 (MQ=12)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          >  1:3970180/1‑93 (MQ=12)                    tttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGc                                                                          >  1:2645437/1‑93 (MQ=12)                     ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                         >  1:2727236/1‑93 (MQ=0)                     ttGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCa                                                                         >  1:3239440/1‑93 (MQ=0)                         tttATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGG‑TTCACTTCTGAGCTCGGCATggg                                                                     >  1:4043236/1‑93 (MQ=9)                                                                    cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                          >  1:82723/1‑93 (MQ=0)                                                                    cccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCttt                          <  2:1349198/93‑1 (MQ=0)                                                                     ccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                         <  2:740363/93‑1 (MQ=9)                                                                     ccACACTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCtttt                         >  1:2846698/1‑93 (MQ=9)                                                                      cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                        >  1:665023/1‑93 (MQ=12)                                                                      cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                        <  2:2152655/93‑1 (MQ=12)                                                                      cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                        <  1:3489282/93‑1 (MQ=16)                                                                      cacaCTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTc                        <  1:4061369/93‑1 (MQ=12)                                                                       acacTACCATCGGCGCTACGGC‑TTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                       <  2:2999117/93‑1 (MQ=255)                                                                       acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                       >  2:581612/1‑93 (MQ=16)                                                                       acacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCt                       >  2:2979513/1‑93 (MQ=16)                                                                        cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                      >  2:3060501/1‑93 (MQ=19)                                                                        cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                      <  2:3041979/93‑1 (MQ=255)                                                                        cacTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTa                      <  2:2779361/93‑1 (MQ=255)                                                                         ccTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                     <  1:4027870/92‑1 (MQ=255)                                                                         acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                     >  2:4536902/1‑93 (MQ=255)                                                                         acTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTaa                     >  2:160770/1‑93 (MQ=255)                                                                          cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                    <  2:1034482/93‑1 (MQ=255)                                                                          cTACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                    <  1:4509163/93‑1 (MQ=255)                                                                           tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc                   <  1:2309170/93‑1 (MQ=255)                                                                           tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc                   >  2:2242987/1‑93 (MQ=255)                                                                           tACCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATc                   <  1:926982/93‑1 (MQ=255)                                                                            aCCATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCt                  <  2:1047528/93‑1 (MQ=255)                                                                             ccATCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTg                 >  1:1109740/1‑93 (MQ=255)                                                                                tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg              <  2:276372/93‑1 (MQ=255)                                                                                tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg              <  1:4130613/93‑1 (MQ=255)                                                                                tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg              <  2:2921791/93‑1 (MQ=255)                                                                                tCGGCGCTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCg              <  2:2190133/93‑1 (MQ=255)                                                                                   gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAAc           <  1:4029513/93‑1 (MQ=255)                                                                                   gcgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGCGGGACCACCGCGCTACTGCCGCCAGACAAATTCTGTTCTAATCTGCCGatc           >  2:109002/1‑91 (MQ=14)                                                                                    cgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt          <  1:2901007/93‑1 (MQ=255)                                                                                    cgcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACt          <  1:1416360/93‑1 (MQ=255)                                                                                     gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGCACtt         >  2:2514304/1‑93 (MQ=255)                                                                                     gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt         >  2:1143449/1‑93 (MQ=255)                                                                                     gcTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACtt         >  2:761585/1‑93 (MQ=255)                                                                                      cTACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACttt        <  1:3752465/93‑1 (MQ=255)                                                                                       tACGGC‑GTTTCGCTTCTGAGTTCGGCATGAGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa       <  2:256442/93‑1 (MQ=17)                                                                                       tACGGC‑GTTTCACTTCTGAGTTCGGCCTGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTACTCCGCCGATCtttc       >  2:4458294/1‑92 (MQ=14)                                                                                       tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa       >  2:3213113/1‑93 (MQ=255)                                                                                       tACGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTa       >  2:3568756/1‑93 (MQ=255)                                                                                        aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa      <  1:488505/93‑1 (MQ=255)                                                                                        aCGGC‑GTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTaa      <  2:540180/93‑1 (MQ=255)                                                                                             ggTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTa  <  1:1579843/93‑1 (MQ=255)                                                                                             |                                                                                             CAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTA  >  NC_000913/3421027‑3421211

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