Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,976,321	C→A	7.6%	noncoding (65/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,321	0	C	A	7.6%	229.2 / 7.7	66	noncoding (65/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 C (30/31);  new base A (2/3);  total (32/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.90e-01

CAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAGGTTTGC  >  NC_000913/2976173‑2976467                                                                                                                                                     |                                                                                                                                                   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                                                                                                                                       <  4:16809/149‑1 (MQ=255)              tGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGgt                                                                                                                                       <  1:714906/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                                                                                                                      >  3:1925/1‑149 (MQ=255)                               tCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCaa                                                                                                                      <  1:19639/149‑1 (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                                                                                                             <  2:105930/149‑1 (MQ=255)                                        cgacTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGtt                                                                                                             <  3:274444/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)                                                                                                                                                     |                                                                                                                                                   CAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTAAATTGTCGGTTACTGTTACAGATTGATGACCGGCAAAAAAAACCTGCGCATCTGCGCAGGCTGGTGTAATTCATGTGCTCAACCCGAAGTTGACTTCACCTATCAATACCTCTGGGATCACCACTTTAGCAACCTGAAGCCAAACGCCACCAGCGGTCAATCGCAACAGCGTTTCGCAAAGTGTAACCAAAGGTTTGC  >  NC_000913/2976173‑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: ‑