Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	525,183	C→T	6.2%	Y641Y (TAC→TAT)	rhsD →	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	525,183	0	C	T	6.2%	200.4 / 2.8	64	Y641Y (TAC→TAT)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/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 = 5.87e-01

AACGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACA  >  NC_000913/525045‑525330                                                                                                                                           |                                                                                                                                                    aaCGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCg                                                                                                                                           >  6:312641/1‑149 (MQ=255)      ggTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACgt                                                                                                                                      <  5:312641/149‑1 (MQ=255)       gTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGAtgctg                                                                                                                                              >  7:100215/1‑140 (MQ=255)       gTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGAtgctg                                                                                                                                              <  8:100215/140‑1 (MQ=255)       gTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACgtg                                                                                                                                     >  4:329809/1‑149 (MQ=255)       gTATGAATACAACGCCGCAGGCGACCAGAATGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGTGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGAcct                                                                                                                                      >  8:88529/1‑146 (MQ=255)          tGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACgt                                                                                                                                      >  3:346176/1‑145 (MQ=255)          tGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACgt                                                                                                                                      <  4:346176/145‑1 (MQ=255)            aaTACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGtcat                                                                                                                                 <  2:340100/148‑1 (MQ=255)             aTACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCACTATGGAGAACGATGCTCCCGGACGTGtcat                                                                                                                                 >  2:358844/1‑147 (MQ=255)                 aaCGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCTCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCt                                                                                                                           >  7:221538/1‑149 (MQ=255)                      cgcAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGAtgctg                                                                                                                                              <  2:50715/125‑1 (MQ=255)                      cgcAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGAtgctg                                                                                                                                              >  1:50715/1‑125 (MQ=255)                       gcAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCaa                                                                                                                     >  4:280905/1‑149 (MQ=255)                        cAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAAc                                                                                                                    <  1:358844/149‑1 (MQ=255)                        cAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAAc                                                                                                                    <  4:112611/149‑1 (MQ=255)                        cAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAAc                                                                                                                    <  8:371530/149‑1 (MQ=255)                        cAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAAc                                                                                                                    <  1:73058/149‑1 (MQ=255)                          ggCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACga                                                                                                                  <  8:221538/149‑1 (MQ=255)                            cGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACgaga                                                                                                                <  8:107761/149‑1 (MQ=255)                                cTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACgg                                                                                                            >  4:185205/1‑149 (MQ=255)                                   aCTGCGGTTATCACCCCGGACGGCAACCGGAGCGAGCCACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCTGGGCGGGCTGACGCGCAGTATTGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCa                                                                                                          <  3:288053/148‑1 (MQ=255)                                        cGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCcaca                                                                                                    <  7:88279/149‑1 (MQ=255)                                          ttATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGc                                                                                                  >  1:306203/1‑149 (MQ=255)                                                            aCCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATgcg                                                                                >  5:211741/1‑149 (MQ=255)                                                             ccGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATgcgc                                                                               >  7:124730/1‑149 (MQ=255)                                                                   cGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGAcc                                                                         >  8:179241/1‑149 (MQ=255)                                                                        cacaGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTg                                                                    <  4:251400/149‑1 (MQ=255)                                                                        cacaGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTg                                                                    <  1:252703/149‑1 (MQ=255)                                                                         acaGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTgg                                                                   >  3:72160/1‑149 (MQ=255)                                                                         acaGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTgg                                                                   >  8:196131/1‑149 (MQ=255)                                                                           aGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTa                                                                 >  5:251719/1‑149 (MQ=255)                                                                                gATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCtt                                                                                             <  4:242806/116‑1 (MQ=255)                                                                                gATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCtt                                                                                             >  3:242806/1‑116 (MQ=255)                                                                                  tGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAgg                                                          <  5:89617/149‑1 (MQ=255)                                                                                     gTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTATGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCgg                                                                       <  6:238750/133‑1 (MQ=255)                                                                                     gTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTATGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCgg                                                                       >  5:211103/1‑133 (MQ=255)                                                                                     gTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTATGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCgg                                                                       <  6:211103/133‑1 (MQ=255)                                                                                     gTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTATGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCgg                                                                       >  5:238750/1‑133 (MQ=255)                                                                                      tGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGggcgg                                                      >  5:156772/1‑149 (MQ=255)                                                                                         ggAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCtt                                                   <  3:280905/149‑1 (MQ=255)                                                                                          gAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCttt                                                  >  4:54225/1‑149 (MQ=255)                                                                                             aGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGAc                                               <  3:185205/149‑1 (MQ=255)                                                                                               gCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACgg                                             <  3:160710/149‑1 (MQ=255)                                                                                                   acAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAAAGAGAACGGCCGCCCCAGCGTCTTCAGTTACGATGCGCTGGACCGGCGGGTACAGCCGGGCGGCTGTGACcggcgg                                         >  8:91050/2‑149 (MQ=21)                                                                                                     aGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGAc                                       <  7:196131/149‑1 (MQ=255)                                                                                                        accacGCAGGGCGGGATGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGc                                     >  7:16449/1‑148 (MQ=255)                                                                                                          cacGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGc                                           >  8:81645/1‑140 (MQ=255)                                                                                                          cacGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGc                                           <  7:81645/140‑1 (MQ=255)                                                                                                          cacGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAAc                                  >  5:114858/1‑149 (MQ=255)                                                                                                              cAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACggg                                            >  3:99210/1‑135 (MQ=255)                                                                                                              cAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACggg                                            <  4:99210/135‑1 (MQ=255)                                                                                                               aGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTAt                             <  3:329809/149‑1 (MQ=255)                                                                                                               aGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTAt                             <  6:156772/149‑1 (MQ=255)                                                                                                               aGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGGCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTAt                             <  6:119937/149‑1 (MQ=255)                                                                                                                     ggcTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTAt                       >  4:288367/1‑149 (MQ=255)                                                                                                                        tGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTAtgac                    <  8:16449/149‑1 (MQ=255)                                                                                                                                  tATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGaa          <  2:306203/149‑1 (MQ=255)                                                                                                                                  tATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGa           >  2:51549/1‑148 (MQ=255)                                                                                                                                     ggAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTg                                                                    >  6:219131/1‑88 (MQ=255)                                                                                                                                     ggAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTg                                                                    <  5:219131/88‑1 (MQ=255)                                                                                                                                      gAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACt      >  4:81069/1‑149 (MQ=255)                                                                                                                                         tACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcac   <  3:288367/149‑1 (MQ=255)                                                                                                                                          aCGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGCCGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTcaca  >  8:322350/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                    AACGCGGTATGAATACAACGCCGCAGGCGACCTGACTGCCGTTATCACCCCGGACGGCAACCGGAGCGAGACACAGTACGATGCGTGGGGAAAGGCGGTCAGCACCACGCAGGGCGGGCTGACGCGCAGTATGGAGTACGATGCTGCCGGACGTGTCATCAGCCTGACCAACGAGAACGGCAGCCACAGCGTCTTCAGTTACGATGCGCTGGACCGGCTGGTACAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACA  >  NC_000913/525045‑525330

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