Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,407,399	C→A	10.6%	intergenic (‑758/‑162)	lrhA ← / → alaA	transcriptional repressor of flagellar, motility and chemotaxis genes/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,407,399	0	C	A	10.6%	226.9 / 14.0	67	intergenic (‑758/‑162)	lrhA/alaA	transcriptional repressor of flagellar, motility and chemotaxis genes/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base A (3/4);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GAAATATCATCATCGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAA  >  NC_000913/2407258‑2407546                                                                                                                                              |                                                                                                                                                    gAAATATCATCATCGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTAATTAttt                                                                                                                                              <  2:67882/149‑1 (MQ=255)    atatCATCATCGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGt                                                                                                                                           <  4:622537/149‑1 (MQ=255)      atcatcatcGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCa                                                                                                                                         >  3:309409/1‑149 (MQ=255)            atcGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCAttatta                                                                                                                                   >  3:78345/1‑149 (MQ=255)             tcGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTGATTATTTTGTCATTATTaa                                                                                                                                  >  4:545515/1‑149 (MQ=255)             tcGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTaa                                                                                                                                  <  4:561841/149‑1 (MQ=255)              cgcgTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAg                                                                                                                                 >  3:658775/1‑149 (MQ=255)                     cATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTaa                                                                                                                          <  2:705061/149‑1 (MQ=255)                     cATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTaa                                                                                                                          <  4:199109/149‑1 (MQ=255)                           tGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTAtttt                                                                                                                                             <  4:285307/124‑1 (MQ=255)                           tGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTAtttt                                                                                                                                             >  3:285307/1‑124 (MQ=255)                           tGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTAATTATTTTGTCATTATTAAGTTATTAACAGcac                                                                                                                    >  4:194809/1‑149 (MQ=255)                           tGAATTATTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAAGTCACTATATAACATTGCATGTAAAGCATATACACCTAATTATTTTGTCATTATTAAGTTATTAACAGcac                                                                                                                    >  4:783511/1‑149 (MQ=255)                              aTTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAAt                                                                                                                 >  4:224884/1‑149 (MQ=255)                                 cTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTCTATAACATTGCATGTAAAGCATATACACCTAATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGa                                                                                                              <  3:783511/149‑1 (MQ=255)                                 cTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTAATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGa                                                                                                              <  3:194809/149‑1 (MQ=255)                                  ttAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAg                                                                                                             >  4:170323/1‑149 (MQ=255)                                             aTTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCctct                                                                                                  >  4:550069/1‑149 (MQ=255)                                             aTTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCctct                                                                                                  <  3:550069/149‑1 (MQ=255)                                              ttttCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTg                                                                                                 >  1:619366/1‑149 (MQ=255)                                              ttttCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTg                                                                                                 >  3:78180/1‑149 (MQ=255)                                                   acacaATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCaaa                                                                                            >  3:258633/1‑149 (MQ=255)                                                   acacaATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCaaa                                                                                            >  1:47207/1‑149 (MQ=255)                                                     acaATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAAt                                                                                          >  1:596596/1‑149 (MQ=255)                                                        aTTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCtt                                                                                       <  2:604545/149‑1 (MQ=255)                                                         ttATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTa                                                                                      <  2:47207/149‑1 (MQ=255)                                                            tCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTc                                                                                   <  4:658775/149‑1 (MQ=255)                                                              ttACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTg                                                                                 >  4:367173/1‑149 (MQ=255)                                                               gtaaGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGc                                                                                >  2:738631/4‑149 (MQ=255)                                                                  aGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTAATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGa                                                                             >  2:74632/1‑149 (MQ=255)                                                                       aaaaaaCCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTc                                                                        <  2:266425/149‑1 (MQ=255)                                                                         aaaaCCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCaa                                                                      >  3:246505/1‑149 (MQ=255)                                                                          aaaCCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCaaa                                                                     <  1:738631/149‑1 (MQ=255)                                                                               aGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACAcc                                                                >  2:399256/1‑149 (MQ=255)                                                                               aGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACAcc                                                                >  3:470496/1‑149 (MQ=255)                                                                                  cTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTc                                                             <  2:619366/149‑1 (MQ=255)                                                                                       aaTTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTggg                                                        >  2:542100/1‑149 (MQ=255)                                                                                        aTTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTgggg                                                       <  1:423309/149‑1 (MQ=255)                                                                                          tAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGa                                                     <  4:32852/149‑1 (MQ=255)                                                                                          tAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGa                                                     <  3:85705/149‑1 (MQ=255)                                                                                             gAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTa                                                  >  1:87807/1‑149 (MQ=255)                                                                                                   gAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTTCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGt                                            >  1:99266/1‑149 (MQ=255)                                                                                                   gAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCa                                                                       <  4:277289/122‑1 (MQ=255)                                                                                                   gAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCa                                                                       >  3:277289/1‑122 (MQ=255)                                                                                                    aaCAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGtt                                           >  3:536163/1‑149 (MQ=255)                                                                                                        aTTCACTATATAACATTGCATGTAAAGCATATACACCTAATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTaa                                       <  2:724804/149‑1 (MQ=255)                                                                                                            aCTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTg                                   >  1:113232/1‑149 (MQ=255)                                                                                                            aCTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCATATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTg                                   >  3:402321/1‑149 (MQ=255)                                                                                                              tatataACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGAc                                 >  3:482674/1‑149 (MQ=255)                                                                                                                  taACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTcc                             <  4:470496/149‑1 (MQ=255)                                                                                                                    aCATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCAc                           <  2:113232/149‑1 (MQ=255)                                                                                                                      aTTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACgg                         >  1:365690/1‑149 (MQ=255)                                                                                                                       ttGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGc                        >  1:772214/1‑149 (MQ=255)                                                                                                                        tGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCa                       <  3:367173/149‑1 (MQ=255)                                                                                                                          cATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAgg                     <  3:224884/149‑1 (MQ=255)                                                                                                                          cATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAgg                     <  1:399256/149‑1 (MQ=255)                                                                                                                            tGTAAAGCATAGACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCAAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAATAGACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGa                   >  4:72032/1‑149 (MQ=255)                                                                                                                               aaaGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGACAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTg                <  2:610372/149‑1 (MQ=255)                                                                                                                                aaGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTgg               >  1:169169/1‑149 (MQ=255)                                                                                                                                aaGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGACCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTgg               <  3:545515/149‑1 (MQ=255)                                                                                                                                     tataCACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGc                                     <  2:721906/122‑1 (MQ=255)                                                                                                                                     tataCACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGc                                     >  1:721906/1‑122 (MQ=255)                                                                                                                                     tataCACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATa          >  2:115719/1‑149 (MQ=255)                                                                                                                                       taCACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAAc        <  2:169169/149‑1 (MQ=255)                                                                                                                                        acacCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACa       <  4:536163/149‑1 (MQ=255)                                                                                                                                        acacCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACa       <  1:542100/149‑1 (MQ=255)                                                                                                                                             tCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCaaa  >  3:250784/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                    GAAATATCATCATCGCGTTTCATTAGTGAATTCTTAAATGAGGCATTTTCACACAATTATCTTACAGATAAAAAAACCAGACTTACAATTAAGAATCAGAACAATTCACTATATAACATTGCATGTAAAGCATATACACCTCATTATTTTGTCATTATTAAGTTATTAACAGCACAATCGAGCCTTCCCCTCTGGCAAAATCTTATTCTGCAGACCTTCAAAACACCGTCCTGGGGGAGTACATTGTTCTAAGCTGACTTCCACGGCAGGGAGTGGCGATAACAGCAAA  >  NC_000913/2407258‑2407546

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