Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,414,992	G→A	6.0%	D377N (GAT→AAT)	acrE →	cytoplasmic membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,414,992	0	G	A	6.0%	210.5 / 2.6	67	D377N (GAT→AAT)	acrE	cytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base A (2/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

AAAAGTCAGGTTGAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAATTGCCCGTCGCTCAGTATC  >  NC_000913/3414845‑3415138                                                                                                                                                    |                                                                                                                                                    aaaaGTCAGGTTGAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGa                                                                                                                                                    >  5:156014/1‑149 (MQ=255)  aaaGTCAGGTTGAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGAt                                                                                                                                                   <  6:226358/149‑1 (MQ=255)   aaGTCAGGTTGAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATa                                                                                                                                                  >  1:248211/1‑149 (MQ=255)   aaGTCAGGTTGAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATa                                                                                                                                                  >  2:116225/1‑149 (MQ=255)    aGTCAGGTTGAAGCGCGCCCTGTCGTTGCCGGTCAGGCGATTGGCGATAAATGGTTGAGTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGGGTCCGGGAGAGCAGGTTAAAGCCCCTACCGAGAc                                                                                                                                                 >  4:114292/1‑149 (MQ=255)             gAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAATGACTGAAATCTGGCGATCAAGTCATTGTAAGCGGCCTGCAAAAAGCGCGTCAGGGAGAGCAGGTTAAAGCCACTACCAATCCACCCg                                                                                                                                            >  8:138750/1‑145 (MQ=255)             gAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCAATACCCCCg                                                                                                                                            >  8:135245/1‑145 (MQ=255)                    gcCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCa                                                                                                                                 >  7:28961/1‑149 (MQ=255)                     cccTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCAt                                                                                                                                <  7:147924/149‑1 (MQ=255)                      ccTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTg                                                                                                                                   >  2:230188/1‑145 (MQ=255)                      ccTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTg                                                                                                                                   <  1:230188/145‑1 (MQ=255)                      ccTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATc                                                                                                                               <  4:117722/149‑1 (MQ=255)                      ccTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATc                                                                                                                               <  7:244059/149‑1 (MQ=255)                       cTGTCGTTGCCAGTCAGGCGATTGGCGATAATTGGTTGATTAGTGAAGGACTGAAATCTGGCTTTCTAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCg                                                                                                                              <  6:88043/149‑1 (MQ=255)                       cTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCg                                                                                                                              <  3:114292/149‑1 (MQ=255)                          tCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGaa                                                                                                                            <  3:33817/148‑1 (MQ=255)                          tCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGaa                                                                                                                            <  2:248211/148‑1 (MQ=255)                          tCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGaa                                                                                                                            <  2:47358/148‑1 (MQ=255)                            gTTGCAATTAAGGCGATTGGAGATAAATGGTTGATTAGTGAAGGACTGACATCTGGCGATCAAGTCATTGTCAGCTGCCTGCAAACAGCGGGTCCGGGAGAGCAGTTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTa                                                                                                                         <  3:10389/149‑1 (MQ=255)                               gCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAgg                                                                                                                      >  7:221864/1‑149 (MQ=255)                               gCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCAATACCCCCGCAGATACTGCATCGAAGTAAg                                                                                                                       <  7:138750/148‑1 (MQ=255)                               gCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCAATACCCCCGCAGATACTGCATCGAAGTAAg                                                                                                                       <  7:135245/148‑1 (MQ=255)                                 cAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGt                                                                                                                          >  5:144505/1‑143 (MQ=255)                                 cAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGt                                                                                                                          <  6:144505/143‑1 (MQ=255)                                   gTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCTTTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAAt                                                                                                                  <  8:26508/149‑1 (MQ=255)                                   gTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCg                                                                                                                              >  2:185263/1‑137 (MQ=255)                                   gTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCg                                                                                                                              <  1:185263/137‑1 (MQ=255)                                       ggCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGa                                                                                                              >  2:114898/1‑149 (MQ=255)                                         cGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACa                                                                                                            <  1:99874/149‑1 (MQ=255)                                         cGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACa                                                                                                            <  8:174866/149‑1 (MQ=255)                                               gCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCa                                                                                                       >  4:191777/1‑148 (MQ=255)                                                      aTGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAC‑Cttttt                                                                                              >  5:179699/1‑149 (MQ=255)                                                         gTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTa                                                                                            <  6:156014/148‑1 (MQ=255)                                                                 gTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCTCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCTCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAACCTTTTTTATTcgacga                                                                                    <  6:199317/149‑1 (MQ=255)                                                                    aaGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCg                                                                                 <  8:28961/148‑1 (MQ=255)                                                                            aaaTCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCa                                                                        >  7:221905/1‑149 (MQ=255)                                                                             aaTCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCAt                                                                       <  4:54494/149‑1 (MQ=255)                                                                             aaTCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCAt                                                                       <  2:30470/149‑1 (MQ=255)                                                                                   tcGATCAAGTCATTGTCAGCGGCCTTCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATgg                                                                     <  1:104556/144‑1 (MQ=255)                                                                                   gCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCTATATTTGCATGGGTg                                                                  <  8:221864/148‑1 (MQ=255)                                                                                   gCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGc                                                                 <  8:183115/149‑1 (MQ=255)                                                                                   gCGATCAAGTCATTGACAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGGAAGGTAAGCTGACATGGCAAC‑CTTTTTTATTCGACGCCCGATATTTGCATGGGTGc                                                                 >  4:155723/1‑149 (MQ=255)                                                                                    cGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCt                                                                >  2:82373/1‑149 (MQ=255)                                                                                              aTTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGAAATGGCAAA‑CTTTTTTATTCGACGCCCGATATTTGCATGGGTGCTGGCCattatt                                                      >  7:112031/1‑149 (MQ=255)                                                                                              aTTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCAACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGAGCGATATTTGCATGGGTGCTGGCCattatt                                                      >  5:241508/1‑149 (MQ=255)                                                                                               ttGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTc                                                     <  1:116225/149‑1 (MQ=255)                                                                                                tGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCt                                                    >  4:225567/1‑149 (MQ=255)                                                                                                   cAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCtgat                                                 <  5:202811/149‑1 (MQ=255)                                                                                                   cAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCtgat                                                 <  1:114898/149‑1 (MQ=255)                                                                                                     gCGGCCTGCAAAAAGAGCGTCCGGGAGAGCAGGTTAAAGCCACTAACGATACCCCCGCAGATACTGCATCGAAGTAATGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGAGCTGGCCATTATTCtgatga                                               >  4:150167/1‑149 (MQ=255)                                                                                                           tGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTc                                                     >  6:5670/1‑137 (MQ=255)                                                                                                           tGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTc                                                     <  5:5670/137‑1 (MQ=255)                                                                                                            gCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGAcc                                                                                  <  6:223978/107‑1 (MQ=255)                                                                                                            gCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGAcc                                                                                  >  5:223978/1‑107 (MQ=255)                                                                                                             cAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATggcgggc                                       >  1:199457/1‑149 (MQ=255)                                                                                                              aaaaaGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATggcgggcg                                      >  5:131412/1‑149 (MQ=255)                                                                                                                    cccgTCCGGGAGATCAGGGTAAAGCCACTACCGATCCCCCCGCAGACACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTgg                                <  6:145079/147‑1 (MQ=255)                                                                                                                     gcgTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGc                               >  3:64676/1‑149 (MQ=255)                                                                                                                     gcgTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGc                               <  3:225567/149‑1 (MQ=255)                                                                                                                          cGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGGGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATcc                          <  3:155723/149‑1 (MQ=255)                                                                                                                           gggAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATcc                          >  5:109446/1‑148 (MQ=255)                                                                                                                           gggAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATcc                          >  5:154974/1‑148 (MQ=255)                                                                                                                                      gtAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAATTAACGTAATCTGACATTGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAATTGCCCg              <  4:121520/148‑1 (MQ=255)                                                                                                                                             cACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAATTGCCCGTCGCTCa       <  1:9389/149‑1 (MQ=255)                                                                                                                                                  ccGCTACCCCCGCAGATACTGCATCGAAGTACGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAATTGCCCGTCGCTCAGTATc  <  2:199457/149‑1 (MQ=255)                                                                                                                                                   cGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAAtt                   >  4:37210/1‑131 (MQ=255)                                                                                                                                                   cGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAAtt                   <  3:37210/131‑1 (MQ=255)                                                                                                                                                    gATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATc                           >  2:105346/1‑122 (MQ=255)                                                                                                                                                    gATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATc                           <  1:105346/122‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAAAGTCAGGTTGAAGCGCGCCCTGTCGTTGCCAGTCAGGCGATTGGCGATAAATGGTTGATTAGTGAAGGACTGAAATCTGGCGATCAAGTCATTGTCAGCGGCCTGCAAAAAGCGCGTCCGGGAGAGCAGGTTAAAGCCACTACCGATACCCCCGCAGATACTGCATCGAAGTAAGGTAATCTGACATGGCAAA‑CTTTTTTATTCGACGACCGATATTTGCATGGGTGCTGGCCATTATTCTGATGATGGCGGGCGCACTGGCGATCCTACAATTGCCCGTCGCTCAGTATC  >  NC_000913/3414845‑3415138

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