Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,072,813	G→T	5.9%	L425L (CTG→CTT)	flu →	CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,072,813	0	G	T	5.9%	237.4 / 2.7	68	L425L (CTG→CTT)	flu	CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTG  >  NC_000913/2072666‑2072954                                                                                                                                                    |                                                                                                                                              caccacCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTg                                                                                                                                               <  5:39273/148‑1 (MQ=255)  accaccGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGa                                                                                                                                              <  6:53193/148‑1 (MQ=255)    caccGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAAc                                                                                                                                            >  6:54659/1‑148 (MQ=255)     accGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACg                                                                                                                                           >  5:322700/1‑148 (MQ=255)         tATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCgg                                                                                                                                       >  3:357893/1‑148 (MQ=255)           tCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTg                                                                                                                                     <  3:302338/148‑1 (MQ=255)            ccATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTgc                                                                                                                                    >  1:383420/1‑147 (MQ=255)             cATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTg                                                                                                                                     >  4:69075/1‑146 (MQ=255)             cATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTg                                                                                                                                     <  3:69075/146‑1 (MQ=255)                gggAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGt                                                                                                                                      <  8:379547/142‑1 (MQ=255)                gggAAATGGCGGTGTACTGCGGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGt                                                                                                                                      >  7:379547/1‑142 (MQ=255)                   aaaTGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGcc                                                                                                                             >  6:103963/1‑148 (MQ=255)                     aTGGCTGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCAc                                                                                                                           >  4:135567/1‑148 (MQ=255)                        gCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGCCGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGa                                                                                                                        >  1:29961/1‑148 (MQ=255)                              tACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCAc                                                                                                                  <  2:56794/148‑1 (MQ=255)                                     ggCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTg                                                                                                                                               >  5:50741/1‑112 (MQ=255)                                     ggCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTg                                                                                                                                               <  6:50741/112‑1 (MQ=255)                                       ccGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATgg                                                                                                         <  5:66023/148‑1 (MQ=255)                                          aTTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCgg                                                                                                      <  4:357893/148‑1 (MQ=255)                                            tCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGAc                                                                                                    >  1:251047/1‑148 (MQ=255)                                                         tCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGcc                                                                                         <  2:29961/146‑1 (MQ=255)                                                         tCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGcc                                                                                         <  6:322700/146‑1 (MQ=255)                                                                      gAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTg                                                                          >  8:262987/1‑148 (MQ=255)                                                                      gAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTg                                                                          >  8:269714/1‑148 (MQ=255)                                                                       aGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTgg                                                                         <  7:21969/148‑1 (MQ=255)                                                                       aGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTgg                                                                         <  2:251047/148‑1 (MQ=255)                                                                           aCGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCggg                                                                     >  6:214858/1‑148 (MQ=255)                                                                            cGGAAAGGCATTCAGTATCGGAGGCGGTCAGTCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGc                                                                    <  5:54659/148‑1 (MQ=255)                                                                            cGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCg                                                                       <  7:131566/145‑1 (MQ=255)                                                                            cGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCg                                                                       >  8:131566/1‑145 (MQ=255)                                                                             ggAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACgg                                                                                                                         >  7:114305/1‑94 (MQ=255)                                                                             ggAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACgg                                                                                                                         <  8:114305/94‑1 (MQ=255)                                                                             ggAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGca                                                                   <  3:135567/148‑1 (MQ=255)                                                                              gAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTaa                                                                                                             <  8:158542/105‑1 (MQ=255)                                                                              gAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTaa                                                                                                             >  7:158542/1‑105 (MQ=255)                                                                                  ggCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcaccacc                                                              <  5:103963/148‑1 (MQ=255)                                                                                   gCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcaccacca                                                             >  5:299453/1‑148 (MQ=255)                                                                                      ttCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGc                                                          <  5:15460/148‑1 (MQ=255)                                                                                       tCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcacc                                                                 <  1:166262/140‑1 (MQ=255)                                                                                       tCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGcacc                                                                 >  2:166262/1‑140 (MQ=255)                                                                                          gTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGaa                                                      >  4:409190/1‑148 (MQ=255)                                                                                           tATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAAt                                                     <  3:409190/148‑1 (MQ=255)                                                                                             tCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATa                                                    <  5:362880/147‑1 (MQ=255)                                                                                              cGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGt                                                                                                 <  2:406148/101‑1 (MQ=255)                                                                                              cGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGt                                                                                                 >  1:406148/1‑101 (MQ=255)                                                                                                   gCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCAc                                                                                                                           >  2:369281/1‑70 (MQ=255)                                                                                                   gCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCAc                                                                                                                           <  1:369281/70‑1 (MQ=255)                                                                                                        cAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATAc                                        >  5:283618/1‑148 (MQ=255)                                                                                                                      gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  1:166444/1‑147 (MQ=255)                                                                                                                      gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  6:28510/1‑147 (MQ=255)                                                                                                                      gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  6:164162/1‑147 (MQ=255)                                                                                                                      gATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  5:236859/1‑147 (MQ=255)                                                                                                                       aTGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           >  6:93337/1‑146 (MQ=255)                                                                                                                       aTGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCgg                           <  5:93337/146‑1 (MQ=255)                                                                                                                            ggAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACg                    >  7:183300/1‑148 (MQ=255)                                                                                                                             gAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACgg                   <  7:262987/148‑1 (MQ=255)                                                                                                                             gAAAAAGGCAGTTCATTCACGCTTAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACgg                   <  7:269714/148‑1 (MQ=255)                                                                                                                             gAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACgg                   <  6:236859/148‑1 (MQ=255)                                                                                                                                 aaGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAcc                                   >  2:124574/1‑128 (MQ=255)                                                                                                                                 aaGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAcc                                   <  1:124574/128‑1 (MQ=255)                                                                                                                                 aaGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaa               <  5:214858/148‑1 (MQ=255)                                                                                                                                     cAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaacaac           >  3:145047/1‑148 (MQ=255)                                                                                                                                      aGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACAGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaacaac           >  5:62444/1‑147 (MQ=255)                                                                                                                                        ttcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGAt        <  2:166444/148‑1 (MQ=255)                                                                                                                                        ttcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGAt        <  5:164162/148‑1 (MQ=255)                                                                                                                                        ttcattcaCGCTGAACGCCGGTGACACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGa         >  7:408153/1‑147 (MQ=255)                                                                                                                                         tcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaac              <  5:215463/141‑1 (MQ=255)                                                                                                                                         tcattcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGaac              >  6:215463/1‑141 (MQ=255)                                                                                                                                             tcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAc                                    <  7:64023/115‑1 (MQ=255)                                                                                                                                             tcaCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTAc                                    >  8:64023/1‑115 (MQ=255)                                                                                                                                               aCGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTg  <  8:183300/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              CACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCACCGCCAGGGGCGGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTG  >  NC_000913/2072666‑2072954

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: ‑