Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,140,400	G→A	5.8%	R405W (CGG→TGG)	asmA ←	suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,140,400	0	G	A	5.8%	225.8 / 2.4	68	R405W (CGG→TGG)	asmA	suppressor of OmpF assembly mutants; putative outer membrane protein assembly factor; inner membrane‑anchored periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTC  >  NC_000913/2140258‑2140548                                                                                                                                               |                                                                                                                                                     gCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTg                                                                                                                                                <  1:160125/149‑1 (MQ=255)  cttGTCCTTGCCAGTTGTGGCGGAATGCGGCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTgg                                                                                                                                               <  2:475971/147‑1 (MQ=255)  ccTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTgg                                                                                                                                               >  4:555674/1‑149 (MQ=255)   cTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGa                                                                                                                                              >  4:577984/1‑149 (MQ=255)     gTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAg                                                                                                                                            >  3:591806/1‑149 (MQ=255)     gTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAg                                                                                                                                            >  4:101656/1‑149 (MQ=255)        cTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGttt                                                                                                                                         >  2:234010/1‑149 (MQ=255)           gCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATc                                                                                                                                      <  1:121759/149‑1 (MQ=255)           gCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATc                                                                                                                                      <  4:344907/149‑1 (MQ=255)            ccAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATcc                                                                                                                                     >  3:289205/1‑149 (MQ=255)                   tgGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGAtt                                                                                                                              >  2:332849/1‑149 (MQ=255)                      tggAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAt                                                                                                                           <  4:458774/148‑1 (MQ=255)                      cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAt                                                                                                                           <  3:555674/149‑1 (MQ=255)                      cGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTAt                                                                                                                           <  3:101656/149‑1 (MQ=255)                        gAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTg                                                                                                                         >  4:413475/1‑149 (MQ=255)                         aaTGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGa                                                                                                                        <  1:536168/149‑1 (MQ=255)                             cGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTc                                                                                                                    <  4:496089/149‑1 (MQ=255)                              gTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcg                                                                                                                   >  1:472118/1‑149 (MQ=255)                               tCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTcgc                                                                                                                  <  2:476512/149‑1 (MQ=255)                                  gCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTc                                                                                                               >  3:583544/1‑149 (MQ=255)                                        aTGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGt                                                                                                         >  1:449549/1‑149 (MQ=255)                                                  cGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGt                                                                                               >  1:487418/1‑149 (MQ=255)                                                   ggAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTg                                                                                              >  3:332088/1‑149 (MQ=255)                                                    gagaAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCGGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGa                                                                                             >  2:200347/1‑149 (MQ=255)                                                       aaaTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAAc                                                                                          <  1:332849/149‑1 (MQ=255)                                                         aTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCt                                                                                        >  1:73282/1‑149 (MQ=255)                                                         aTCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCt                                                                                        >  4:18310/1‑149 (MQ=255)                                                           cACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGt                                                                                      <  2:346298/149‑1 (MQ=255)                                                           cACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGt                                                                                      <  1:24630/149‑1 (MQ=255)                                                           cACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGt                                                                                      >  4:371969/1‑149 (MQ=255)                                                            aCCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGGTGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTc                                                                                     <  3:144264/149‑1 (MQ=255)                                                               agccagTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCAc                                                                                  >  4:327822/1‑149 (MQ=255)                                                                     tGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTaa                                                                            >  3:152122/1‑149 (MQ=255)                                                                     tGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGATTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTaa                                                                            <  2:376389/149‑1 (MQ=255)                                                                      gACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCTGATTTATTGATGTCGCGACCAGCTTGCCCGGCAGTGAAACCTGTCCACCGTTAAg                                                                           <  2:472118/149‑1 (MQ=255)                                                                       aCATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGt                                                                          >  1:536877/1‑149 (MQ=255)                                                                        cATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGtt                                                                         >  1:524035/1‑149 (MQ=255)                                                                                    cAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGtt                                                             <  1:7225/149‑1 (MQ=255)                                                                                      aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCGAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                           <  2:533386/149‑1 (MQ=255)                                                                                      aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                           <  1:249834/149‑1 (MQ=255)                                                                                      aCGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                           <  1:557501/149‑1 (MQ=255)                                                                                        gAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGt                                                         >  2:101739/1‑149 (MQ=255)                                                                                           aTCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAAt                                                      <  4:152122/149‑1 (MQ=255)                                                                                                aTAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCa                                                 <  3:327822/149‑1 (MQ=255)                                                                                                    ttAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCaa                                             >  3:316646/1‑149 (MQ=255)                                                                                                      aaaCGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAc                                           >  2:299160/1‑149 (MQ=255)                                                                                                        aCGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAAccc                                         <  4:289205/149‑1 (MQ=255)                                                                                                          gCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTa                                                                                                                            >  1:526605/1‑64 (MQ=255)                                                                                                          gCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTa                                                                                                                            <  2:526605/64‑1 (MQ=255)                                                                                                           ccTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGAc                                      <  1:110633/149‑1 (MQ=255)                                                                                                             ttCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                           <  2:455101/126‑1 (MQ=255)                                                                                                             ttCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGa                                                           >  1:455101/1‑126 (MQ=255)                                                                                                                aGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 <  1:87100/149‑1 (MQ=255)                                                                                                                aGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 <  4:125405/149‑1 (MQ=255)                                                                                                                aGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 >  1:534063/1‑149 (MQ=255)                                                                                                                 gAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgttg                                <  2:187336/149‑1 (MQ=255)                                                                                                                 gAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgttg                                <  4:583544/149‑1 (MQ=255)                                                                                                                                cAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCa                 <  2:524035/149‑1 (MQ=255)                                                                                                                                 aaCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTaa                                                                            >  1:523407/1‑89 (MQ=255)                                                                                                                                 aaCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTaa                                                                            <  2:523407/89‑1 (MQ=255)                                                                                                                                  aCGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAAc               >  3:414685/1‑149 (MQ=255)                                                                                                                                    gTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 >  3:518407/1‑129 (MQ=255)                                                                                                                                    gTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACttgtt                                 <  4:518407/129‑1 (MQ=255)                                                                                                                                     acttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCa                                                 >  2:396804/3‑112 (MQ=255)                                                                                                                                     acttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGGCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCa                                                 <  1:396804/110‑1 (MQ=255)                                                                                                                                       ttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGt                                                                          <  3:68811/85‑1 (MQ=255)                                                                                                                                       ttCCAGCCACGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGt                                                                          >  4:68811/1‑85 (MQ=255)                                                                                                                                            gCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTaaaa     <  2:449549/149‑1 (MQ=255)                                                                                                                                               gTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAattc  <  4:332088/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     GCCTGTCCTTGCCAGTTGTGGCGGAATGCGTCGGCATCTATGTCAGCACCGGAGAAATCACCAGCCAGTGACATTTTTCCGGTCAACGAAATCGGATAGTTAAACGCCTTCAGAATGGTACCAATCTCAACGTTTTCCAGCCGTGGCTGGAAGTTTATCCGCGGATTTATTGATGTCGCGTCCAGCGTGCCCGGCAGTGAAACCTGTCCACCGTTAAGTTTGCCCTGCAGTTGAGTAATTTCCAGCAAACCCGACTTGTTGGTCATTTGCGTGGCAACATCTGTAAAATTC  >  NC_000913/2140258‑2140548

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