Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,104,931	G→T	11.1%	intergenic (+29/‑30)	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,931	0	G	T	11.1%	215.7 / 10.0	68	intergenic (+29/‑30)	csgA/csgC	curlin, major subunit/inhibitor of CsgA amyloid formation
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base T (2/6);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104783‑1105077                                                                                                                                                     |                                                                                                                                                    aTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcgc                                                                                                                                                    >  1:651997/1‑150 (MQ=255) aTGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGgcgc                                                                                                                                                    <  1:146202/150‑1 (MQ=255)       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:219955/150‑7 (MQ=255)                  cGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCTTGCGCCCTG‑tttctgtaa                                                                                                                                  <  2:186076/150‑7 (MQ=255)                   ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                 >  2:123841/1‑150 (MQ=255)                   ggtggtGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCg                                                                                                                                 >  2:418808/1‑150 (MQ=255)                    gtggtgGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCgg                                                                                                                                >  2:1220256/1‑150 (MQ=255)                     tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                               <  2:891713/150‑1 (MQ=255)                     tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                               <  2:1076529/150‑1 (MQ=255)                     tggtggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCggg                                                                                                                               <  2:202143/150‑1 (MQ=255)                         ggCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaa                                                                                                                           >  2:447103/1‑150 (MQ=255)                          gCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaag                                                                                                                          <  2:561131/150‑1 (MQ=255)                          gCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGagaag                                                                                                                          >  2:126083/1‑150 (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                                                                       >  1:1037154/1‑150 (MQ=255)                                                                             cTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCttt                                                                       <  2:981886/150‑1 (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)                                                                                                                                                     |                                                                                                                                                    ATGACGGTTAAACAGTTCGGTGGTGGCAACGGTGCTGCAGTTGACCAGACTGCATCTAACTCCTCCGTCAACGTGACTCAGGTTGGCTTTGGTAACAACGCGACCGCTCATCAGTACTAATACATCATTTGTATTACAGAAACAGGGCGCAAGCCCTG‑TTTTTTTTCGGGAGAAGAATATGAATACGTTATTACTCCTTGCGGCACTTTCCAGTCAGATAACCTTTAATACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGC  >  NC_000913/1104783‑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: ‑