Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,976,320	T→G	5.4%	noncoding (66/82 nt)	omrB ←	sRNA antisense regulator downregulates OM proteins and curli; positively regulated by OmpR/EnvZ, Hfq‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,976,320	0	T	G	5.4%	229.8 / 4.1	69	noncoding (66/82 nt)	omrB	sRNA antisense regulator downregulates OM proteins and curli; positively regulated by OmpR/EnvZ, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base G (2/2);  total (35/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 = 4.11e-01

CCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAGGTTTGC  >  NC_000913/2976172‑2976467                                                                                                                                                     |                                                                                                                                                    ccAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCAt                                                                                                                                                     >  2:215231/1‑149 (MQ=255)  caatcaatACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATc                                                                                                                                                    >  2:714906/1‑149 (MQ=255)  caatcaatACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATc                                                                                                                                                    >  4:220976/1‑149 (MQ=255)  caatcaatACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATc                                                                                                                                                    >  4:321858/1‑149 (MQ=255)        atACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCa                                                                                                                                              >  1:768728/1‑149 (MQ=255)         tACCTCTGGGATCTTGATTGTTGTGTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCAGATGCGCAg                                                                                                                                             <  2:164768/149‑1 (MQ=255)         tACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAg                                                                                                                                             <  3:220976/149‑1 (MQ=255)          aCCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAgg                                                                                                                                            <  2:534736/149‑1 (MQ=255)           ccTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCAGATGCGCAggt                                                                                                                                           <  2:659706/149‑2 (MQ=255)            ctctGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCt                                                                                                                                          >  2:17242/1‑149 (MQ=255)             tctGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTg                                                                                                                                         <  2:362435/149‑1 (MQ=255)               tGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGgt                                                                                                                                       >  3:703752/1‑149 (MQ=255)               tGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGgt                                                                                                                                       <  1:714906/149‑1 (MQ=255)               tGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGgt                                                                                                                                       <  4:16809/149‑1 (MQ=255)               tGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGgt                                                                                                                                       <  1:17242/149‑1 (MQ=255)                   atattGATTGTGGTCTTCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCACAAAAAACCTGCGCATCTGCGCAGGCTGGTGTaa                                                                                                                                   <  3:480866/146‑1 (MQ=255)                   aTCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTaa                                                                                                                                   >  4:243989/1‑149 (MQ=255)                    tCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAAt                                                                                                                                  >  1:105930/1‑149 (MQ=255)                        gATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCAt                                                                                                                              <  3:449934/149‑1 (MQ=255)                          ttGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCAtgt                                                                                                                            <  4:257790/149‑1 (MQ=255)                           tgtgGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCAtgtg                                                                                                                           <  1:702147/149‑1 (MQ=255)                               gTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCa                                                                                                                       >  4:275507/1‑149 (MQ=255)                                tCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCaa                                                                                                                      <  1:19639/149‑1 (MQ=255)                                tCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCaa                                                                                                                      >  3:1925/1‑149 (MQ=255)                                  tGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAAcc                                                                                                                    <  1:414367/149‑1 (MQ=255)                                      cgacgaCTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGaa                                                                                                                <  4:1925/149‑1 (MQ=255)                                       gacgacTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAg                                                                                                               <  4:504465/149‑1 (MQ=255)                                       gacgacTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAg                                                                                                               <  2:230077/149‑1 (MQ=255)                                         cgacTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGtt                                                                                                             <  3:274444/149‑1 (MQ=255)                                         cgacTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGtt                                                                                                             <  2:105930/149‑1 (MQ=255)                                             tctTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACt                                                                                                         >  1:185076/1‑149 (MQ=255)                                             tctTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACt                                                                                                         >  3:425126/1‑149 (MQ=255)                                                       cGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCa                                                                                               >  2:147429/1‑149 (MQ=255)                                                        gagaAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCaa                                                                                              <  1:147429/149‑1 (MQ=255)                                                          gaAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTATGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATa                                                                                            >  2:677037/1‑149 (MQ=255)                                                           aaaaCGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATAc                                                                                           >  2:263281/1‑149 (MQ=255)                                                           aaaaCGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATAc                                                                                           <  4:703752/149‑1 (MQ=255)                                                                  aaaGGAATGAAGGGAAATGCAACGAGGTATGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTggg                                                                                    <  1:677037/149‑1 (MQ=255)                                                                     ggAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATc                                                                                 <  3:108645/149‑1 (MQ=255)                                                                           aaGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACt                                                                           >  3:402989/1‑149 (MQ=255)                                                                                                 tAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGc                                                     >  3:214731/1‑149 (MQ=255)                                                                                                 tAAATTGTAGGTTACTGTTACAGATTGATGAGCGGCAAAAAGAACCTGCGCATCTGCGCAGGTTGGTGTACTCCATGTGCTCAACCGGAAGTTGACTTCACCTATTAATTCCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGc                                                     <  3:468337/149‑1 (MQ=255)                                                                                                   aaTTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCAGATGCGCAggttt                                                                                                                                         >  4:57240/1‑60 (MQ=38)                                                                                                   aaTTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCAAATGCGCAggttt                                                                                                                                         <  3:57240/63‑4 (MQ=255)                                                                                                    aTTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGccac                                                  <  2:97317/149‑1 (MQ=255)                                                                                                           gTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTCCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCCCCACTTTAGCAACCTGAAGCCAAACGCAACCAGCGGt                                           >  3:27247/1‑149 (MQ=255)                                                                                                                  ttACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGc                                    >  4:675887/1‑149 (MQ=255)                                                                                                                  ttACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGc                                    >  4:263093/1‑149 (MQ=255)                                                                                                                  ttACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGc                                    >  2:99202/1‑149 (MQ=255)                                                                                                                  ttACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGc                                    >  1:101754/1‑149 (MQ=255)                                                                                                                    aCAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCaa                                  <  4:299647/149‑1 (MQ=255)                                                                                                                        aTTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGc                              <  4:425126/149‑1 (MQ=255)                                                                                                                         ttGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCg                             >  1:729141/1‑149 (MQ=255)                                                                                                                           gatgaCCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGtt                           <  1:219799/149‑1 (MQ=255)                                                                                                                             tgaCCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTc                         <  2:768728/149‑1 (MQ=255)                                                                                                                             tgaCCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTc                         <  2:729141/149‑1 (MQ=255)                                                                                                                                ccGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCa                      >  3:306840/1‑149 (MQ=255)                                                                                                                                    cAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAgt                  <  4:214731/149‑1 (MQ=255)                                                                                                                                     aaaaaaaaCCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAgtg                 >  4:82899/1‑149 (MQ=255)                                                                                                                                       aaaaaaCCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTa               >  2:695966/1‑149 (MQ=255)                                                                                                                                       aaaaaaCCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTa               >  3:359212/1‑149 (MQ=255)                                                                                                                                        aaaaaCCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTaa              >  1:641404/1‑149 (MQ=255)                                                                                                                                         aaaaCCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAAc             >  2:384623/1‑149 (MQ=255)                                                                                                                                          cagccTGCGCAGATGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAAcc            >  4:468777/4‑149 (MQ=255)                                                                                                                                             ccTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCaaa         <  3:89151/149‑1 (MQ=255)                                                                                                                                               tGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGttt                          <  2:343889/130‑1 (MQ=255)                                                                                                                                               tGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGttt                          >  1:343889/1‑130 (MQ=255)                                                                                                                                               tGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAgg       >  1:299964/1‑149 (MQ=255)                                                                                                                                                    aTCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAGGTTTGc  >  3:546089/1‑149 (MQ=255)                                                                                                                                                    aTCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAGGTTTGc  >  3:125238/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAGGTTTGC  >  NC_000913/2976172‑2976467

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