Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,104,933	A→G	10.3%	intergenic (+31/‑28)	csgA → / → csgC	curlin, major subunit/inhibitor of CsgA amyloid formation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,104,933	0	A	G	10.3%	187.7 / 10.7	68	intergenic (+31/‑28)	csgA/csgC	curlin, major subunit/inhibitor of CsgA amyloid formation
Reads supporting (aligned to +/- strand):  ref base A (28/33);  new base G (2/5);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104789‑1105077                                                                                                                                                 |                                                                                                                                                  gTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGccc                                                                                                                                              <  2:1259624/150‑1 (MQ=255)    aaaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑t                                                                                                                                          >  2:27115/1‑150 (MQ=255)     aaCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tt                                                                                                                                         >  2:145418/1‑150 (MQ=255)      aCAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑ttt                                                                                                                                        >  1:192836/1‑150 (MQ=255)        aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttct                                                                                                                                      <  2:192836/150‑3 (MQ=255)        aGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttct                                                                                                                                      <  1:39435/150‑3 (MQ=255)         gTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCGGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttctg                                                                                                                                     <  2:233658/150‑4 (MQ=255)         gTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tttttt                                                                                                                                     <  1:758167/150‑1 (MQ=255)           tCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑tttttttt                                                                                                                                   <  2:757397/150‑1 (MQ=255)            cGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttctgtaa                                                                                                                                  <  2:186076/150‑7 (MQ=255)            cGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttctgtaa                                                                                                                                  <  2:219955/150‑7 (MQ=255)             ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                 >  2:418808/1‑150 (MQ=255)             ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                 >  2:123841/1‑150 (MQ=255)              gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCgg                                                                                                                                >  2:1220256/1‑150 (MQ=255)               tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                               <  2:1076529/150‑1 (MQ=255)               tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                               <  2:202143/150‑1 (MQ=255)               tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                               <  2:891713/150‑1 (MQ=255)                   ggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaa                                                                                                                           >  2:447103/1‑150 (MQ=255)                    gCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaag                                                                                                                          >  2:126083/1‑150 (MQ=255)                    gCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaag                                                                                                                          <  2:561131/150‑1 (MQ=255)                           tgctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGa                                                                                                                   <  1:126083/150‑1 (MQ=255)                            gctgcAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGaa                                                                                                                  <  2:1108904/150‑1 (MQ=255)                                        cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTGTTTTTTTTTCGGGAGAAGAATATGAATACGTTATTAc                                                                                                       <  1:1001665/150‑1 (MQ=255)                                        cAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACt                                                                                                      >  1:1280494/1‑150 (MQ=255)                                            cTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCtt                                                                                                  >  2:20603/1‑150 (MQ=255)                                            cTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCtt                                                                                                  >  1:127460/1‑150 (MQ=255)                                              gCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGc                                                                                                <  2:1191271/150‑1 (MQ=255)                                               cATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACATGGCGCAAGCCCTG‑TTTTTTTTCGTGAGAAGAATATGAATACGTTATTACTCTTTTCg                                                                                               <  2:1031004/150‑1 (MQ=255)                                                aTCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCgg                                                                                              <  1:145418/150‑1 (MQ=255)                                                  cTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCa                                                                                            >  2:865665/1‑150 (MQ=255)                                                         ctccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCa                                                                                     >  2:1398066/1‑150 (MQ=255)                                                          tccGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAg                                                                                    <  1:418808/150‑1 (MQ=255)                                                            cGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTc                                                                                  >  1:539011/1‑150 (MQ=255)                                                                aaCGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGAt                                                                              >  2:46658/1‑150 (MQ=255)                                                                   gTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAAcc                                                                          >  2:159238/1‑150 (MQ=255)                                                                    tGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAAcc                                                                          <  2:1152994/150‑1 (MQ=255)                                                                      aCTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCtt                                                                        <  2:864472/150‑1 (MQ=255)                                                                       cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                       <  2:981886/150‑1 (MQ=255)                                                                       cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                       >  1:1037154/1‑150 (MQ=255)                                                                          aGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAAt                                                                    <  1:265319/150‑1 (MQ=255)                                                                               ggCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAc                                                               <  2:1280494/150‑1 (MQ=255)                                                                                gCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGAcc                                                              <  1:27115/150‑1 (MQ=255)                                                                                         acaacGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAggg                                                     >  1:1363035/1‑150 (MQ=255)                                                                                             cgcgACCTCTCATCAGTACTAATACATCATTTGTATTACAGAAACTGGTCGCAATGCCTG‑TTTTTTTTCGGGAGAATAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtg                                                 <  2:701214/150‑1 (MQ=255)                                                                                                   cGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGtata                                           <  2:229214/150‑1 (MQ=255)                                                                                                       catcaGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCat                                       >  1:1111269/1‑150 (MQ=255)                                                                                                          caGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCattat                                    >  2:1216946/1‑150 (MQ=255)                                                                                                            gTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCA‑TATTcc                                 <  1:1398066/150‑1 (MQ=255)                                                                                                             tACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTcc                                 <  2:1037154/150‑1 (MQ=255)                                                                                                             tACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTcc                                 <  2:539011/150‑1 (MQ=255)                                                                                                               cTAATACATAATTTGTATTACATAAACAGGGCTCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGTCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTg                               <  2:651997/150‑1 (MQ=255)                                                                                                                tAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑‑TTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGaa                             <  1:159238/150‑1 (MQ=255)                                                                                                                  aTACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAg                            <  1:46658/150‑1 (MQ=255)                                                                                                                     catcatTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCa                         >  1:1308466/1‑150 (MQ=255)                                                                                                                         atTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCActct                     >  1:984324/1‑150 (MQ=255)                                                                                                                            tGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTAc                  <  1:276818/150‑1 (MQ=255)                                                                                                                              tATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTc                <  1:1220256/150‑1 (MQ=255)                                                                                                                               aTTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCa               >  1:1023926/1‑150 (MQ=255)                                                                                                                                ttACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCaa              >  1:921801/1‑150 (MQ=255)                                                                                                                                  aaaaaaaCAGGGCTTGCGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc            >  2:207899/5‑150 (MQ=255)                                                                                                                                  aCAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATc            >  2:868079/1‑150 (MQ=255)                                                                                                                                  aCAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTACTGAAGTCACTCTTACTCAATc            >  2:862841/1‑150 (MQ=255)                                                                                                                                   cAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCt           >  1:997228/1‑150 (MQ=255)                                                                                                                                     gAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                               <  2:525758/112‑1 (MQ=255)                                                                                                                                     gAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGAtgtg                                               >  1:525758/1‑112 (MQ=255)                                                                                                                                     gAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTg         >  2:244180/1‑150 (MQ=255)                                                                                                                                            ggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc  <  2:1111269/150‑1 (MQ=255)                                                                                                                                            ggCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGc  <  2:921801/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104789‑1105077

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