Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	716,572	T→C	5.9%	pseudogene (26/651 nt)	ybfG ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	716,572	0	T	C	5.9%	216.6 / 3.6	68	pseudogene (26/651 nt)	ybfG	pseudogene
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACATTGCCACTTTTT  >  NC_000913/716425‑716720                                                                                                                                                    |                                                                                                                                                     aCAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTAtt                                                                                                                                                     >  4:472927/1‑149 (MQ=255) aCAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTAtt                                                                                                                                                     >  4:5022/1‑149 (MQ=255)   aaGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTa                                                                                                                                                   >  2:584002/1‑149 (MQ=255)       gcACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTaa                                                                                                                                               <  4:187845/149‑1 (MQ=255)        cACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAAc                                                                                                                                              <  1:41181/149‑1 (MQ=255)         aCCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACt                                                                                                                                             <  1:525632/149‑1 (MQ=255)            cATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGAt                                                                                                                                          >  2:141749/1‑149 (MQ=255)                ggTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCAt                                                                                                                                      >  2:130819/1‑149 (MQ=255)                 gTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCAtt                                                                                                                                     >  1:137696/1‑149 (MQ=255)                 gTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCAtt                                                                                                                                     <  2:54705/149‑1 (MQ=255)                   gATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTc                                                                                                                                   <  2:555360/149‑1 (MQ=255)                    aTTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGa                                                                                                                                           <  2:291786/140‑1 (MQ=255)                    aTTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGa                                                                                                                                           >  1:291786/1‑140 (MQ=255)                    aTTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCa                                                                                                                                  <  4:480682/149‑1 (MQ=255)                     ttttGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCAc                                                                                                                                 >  2:412309/1‑149 (MQ=255)                     ttttGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCAc                                                                                                                                 >  4:394075/1‑149 (MQ=255)                            gtCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGc                                                                                                                          <  1:252431/149‑1 (MQ=255)                            gtCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTACTTACTAACTGATTCATTTCACGCATCGc                                                                                                                          >  2:102111/1‑149 (MQ=255)                            gtCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTACTTACTAACTGATTCATTTCACGCATCGc                                                                                                                          >  4:223377/1‑149 (MQ=255)                             tCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCt                                                                                                                         >  4:324873/1‑149 (MQ=255)                               aCCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGTCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTAc                                                                                                                       >  2:264025/1‑149 (MQ=255)                                ccTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACt                                                                                                                      <  4:517220/149‑1 (MQ=255)                                  tGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTg                                                                                                                    >  1:365560/1‑149 (MQ=255)                                      aTTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTcc                                                                                                                <  3:377588/149‑1 (MQ=255)                                          aTAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTtgt                                                                                                            >  2:104214/1‑149 (MQ=255)                                           tAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTtgtg                                                                                                           >  4:135236/1‑149 (MQ=255)                                                tGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTa                                                                                                      <  1:584002/149‑1 (MQ=255)                                                          gCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCAtt                                                                                            <  4:246021/149‑1 (MQ=255)                                                            tAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTa                                                                                          >  1:74038/1‑149 (MQ=255)                                                            tAAATCTTTATTAACCGGTAAGAGTCCGGTATGGTTTGTGGCGTCTTTTAAACAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTa                                                                                          <  1:259209/149‑1 (MQ=255)                                                               aTCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTatga                                                                                       <  1:505577/149‑1 (MQ=255)                                                                tCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTatgat                                                                                      <  2:137696/149‑1 (MQ=255)                                                                   ttattaACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAg                                                                                   >  1:554616/1‑149 (MQ=255)                                                                   ttattaACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACCGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAg                                                                                   <  1:500851/149‑1 (MQ=255)                                                                   ttattaACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTACTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAg                                                                                   <  3:223377/149‑1 (MQ=255)                                                                   ttattaACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTACTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAg                                                                                   <  1:102111/149‑1 (MQ=255)                                                                    tattaACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAgg                                                                                  >  1:152612/1‑149 (MQ=255)                                                                        aaCCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAAtt                                                                              <  4:3625/149‑1 (MQ=255)                                                                          ccGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGt                                                                            >  2:325772/1‑149 (MQ=255)                                                                          ccGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGt                                                                            >  2:559743/1‑149 (MQ=255)                                                                             gTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATg                                                                         >  1:343949/1‑149 (MQ=255)                                                                             gTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATg                                                                         <  3:481318/149‑1 (MQ=255)                                                                                 gagTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAAc                                                                     >  2:454868/1‑149 (MQ=255)                                                                                  agTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACa                                                                    >  3:86539/1‑149 (MQ=255)                                                                                    ttCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAg                                                                                   <  4:146375/132‑1 (MQ=255)                                                                                    ttCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAg                                                                                   >  3:146375/1‑132 (MQ=255)                                                                                         tATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTa                                                             >  3:481175/1‑149 (MQ=255)                                                                                           tGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTaaa                                                           <  2:74038/149‑1 (MQ=255)                                                                                                    gCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGtt                                                  >  4:37545/1‑149 (MQ=255)                                                                                                    gCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGtt                                                  >  3:344486/1‑149 (MQ=255)                                                                                                           ttAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCCACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTtata                                           <  1:130819/149‑1 (MQ=255)                                                                                                              aCCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATAttt                                        <  1:104214/149‑1 (MQ=255)                                                                                                              aCCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATAttt                                        >  3:247245/1‑149 (MQ=255)                                                                                                              aCCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATAttt                                        >  1:407918/1‑149 (MQ=255)                                                                                                                   ttCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATa                                   >  2:359498/1‑149 (MQ=255)                                                                                                                      tGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATaaaa                                >  1:132175/1‑149 (MQ=255)                                                                                                                      tGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTGATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATaaaa                                >  1:187531/1‑149 (MQ=255)                                                                                                                        gATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGCTTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTAGGGTGTGATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTTAATaaaaaa                              <  3:324873/149‑1 (MQ=255)                                                                                                                         attattGGCCAGTCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAAt                             <  1:462345/149‑1 (MQ=255)                                                                                                                           tattGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATTGCTACTTGCTCCTTGTGTGCTACTGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGa                           <  4:178888/149‑1 (MQ=255)                                                                                                                           tattGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGa                           >  2:265246/1‑149 (MQ=255)                                                                                                                             ttGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACt                         <  2:343949/149‑1 (MQ=255)                                                                                                                                   aGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATcc                   <  2:159442/149‑1 (MQ=255)                                                                                                                                        gAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACa              >  2:299198/1‑149 (MQ=255)                                                                                                                                          aTCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATa                >  4:292797/1‑145 (MQ=255)                                                                                                                                          aTCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATa                <  3:292797/145‑1 (MQ=255)                                                                                                                                          aTCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACAtt            <  1:141749/149‑1 (MQ=255)                                                                                                                                             ttCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACATTGcc         >  1:471441/1‑149 (MQ=255)                                                                                                                                             ttCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACATTGcc         >  4:275816/1‑149 (MQ=255)                                                                                                                                              tCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGTTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAGATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACATTGCCa        <  4:481175/149‑1 (MQ=255)                                                                                                                                                    tttACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACATTGCCACttttt  <  1:265246/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTGGCGTCTTTTAACCAGTTCTGGATTATTGGCCAGGCTGAATCTTCAGTATTTACTAACTGATTCATTTCACGCATCGCTACTTGCTCCTTGTGTGCTACGGTGTCATTTATGATGAGGAATTGTATGGAACATTTTTTAAATATTAAGTTCCTTATATTTGAATAAAAAATGACTTTATCCATACATTGCCACTTTTT  >  NC_000913/716425‑716720

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