BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,113,378	A→T	6.3%	L34I (TTA→ATA)	yghG ←	lipoprotein YghG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,113,378	0	A	T	6.3%	435.3 / 12.0	137	L34I (TTA→ATA)	yghG	lipoprotein YghG
Reads supporting (aligned to +/- strand): ref base A (64/64); new base T (7/2); total (71/66)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 140 total aligned reads displayed.

CAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATA  >  NC_000913/3113231‑3113527

  ggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                                                                                   <  1:1049602/150‑1 (MQ=255)

  ggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                                                                                   >  1:1450732/1‑150 (MQ=255)

         gCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTTATCCTGTTTTTTCGCCAGTAAACTTTCAtttt                                                                                                                                            <  2:1396138/150‑1 (MQ=255)

         gCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAtttt                                                                                                                                            >  1:550050/1‑150 (MQ=255)

         gCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAtttt                                                                                                                                            >  2:1664104/1‑150 (MQ=255)

          cGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTc                                                                                                                                           <  2:1515591/150‑1 (MQ=255)

           gTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCa                                                                                                                                          <  1:43885/150‑1 (MQ=255)

             ggTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCAtt                                                                                                                                        <  2:1841721/150‑1 (MQ=255)

                ttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtg                                                                                                                                     >  2:1587723/1‑150 (MQ=255)

                ttGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtg                                                                                                                                     <  2:1017528/150‑1 (MQ=255)

                 tgtgTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtgg                                                                                                                                    >  2:92731/1‑150 (MQ=255)

                     tACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                                                                                <  1:737906/150‑1 (MQ=255)

                      aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                                                                               <  1:624082/150‑1 (MQ=255)

                      aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                                                                               >  2:535053/1‑150 (MQ=255)

                      aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctggc                                                                                                                               <  2:74954/150‑1 (MQ=255)

                       cccGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGca                                                                                                                              >  2:291309/1‑150 (MQ=255)

                        ccGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGcac                                                                                                                             <  2:1171609/150‑1 (MQ=255)

                             tCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCa                                                                                                                        <  2:172645/150‑1 (MQ=255)

                               acTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGTgccagccata                                                                                                                      <  1:1561065/149‑3 (MQ=255)

                                cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACtt                                                                                                                     >  1:10173/1‑150 (MQ=255)

                                cTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTAGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACtt                                                                                                                     >  2:1339896/1‑150 (MQ=255)

                                 tGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTta                                                                                                                    <  1:1745129/150‑1 (MQ=255)

                                  gATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGTgccagccataatg                                                                                                                   <  1:1498446/150‑6 (MQ=255)

                                   aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaat                                                                                                                  <  1:214248/150‑1 (MQ=255)

                                    tAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaata                                                                                                                 >  2:1457682/1‑150 (MQ=255)

                                      aTGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTaataat                                                                                                               >  1:1435589/1‑150 (MQ=255)

                                        ggTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                                                                             <  2:1083995/150‑1 (MQ=255)

                                        ggTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                                                                             <  1:1885882/150‑1 (MQ=255)

                                        ggTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                                                                             >  2:547192/1‑150 (MQ=255)

                                         gTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCt                                                                                                            >  1:937597/1‑150 (MQ=255)

                                          tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTg                                                                                                           <  2:766614/150‑1 (MQ=255)

                                           cATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGt                                                                                                          >  2:1270053/1‑150 (MQ=255)

                                           cATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGt                                                                                                          >  1:1317540/1‑150 (MQ=255)

                                            aTTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTa                                                                                                         <  1:1680233/150‑1 (MQ=255)

                                                  aCCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTc                                                                                                   <  1:456129/150‑1 (MQ=255)

                                                          gCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAg                                                                                           >  2:744681/1‑150 (MQ=255)

                                                           ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGc                                                                                          >  1:1949122/1‑150 (MQ=255)

                                                               tACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGTCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGata                                                                                      <  2:243107/150‑1 (MQ=255)

                                                                aCTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGatat                                                                                     >  2:1720948/1‑150 (MQ=255)

                                                                      cgcCAGCACTAAGGTATATCCCGCAGAATTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCAACTTAATAATCCTGTAACGCTCAACAATAGCTATATTAATAc                                                                               >  2:1883634/1‑150 (MQ=255)

                                                                       gcCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAct                                                                              <  1:1587723/150‑2 (MQ=255)

                                                                       gcCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                                                                              >  1:539159/1‑150 (MQ=255)

                                                                       gcCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                                                                              <  2:418316/150‑1 (MQ=255)

                                                                          aGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctc                                                                           >  1:381314/1‑150 (MQ=255)

                                                                            cACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccc                                                                         >  1:122261/1‑150 (MQ=255)

                                                                              cTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGGCAGTAAACTGGCATTTTCATTATTGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg                                                                       <  2:1131007/150‑1 (MQ=255)

                                                                              cTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg                                                                       <  1:1445989/150‑1 (MQ=255)

                                                                              cTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg                                                                       <  1:1664104/150‑1 (MQ=255)

                                                                               tAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTgg                                                                      >  1:1749547/1‑150 (MQ=255)

                                                                                aaGTTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCTCCAGTAAACTGGCATTTTGATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGc                                                                     <  2:1288329/150‑1 (MQ=255)

                                                                                aaGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGc                                                                     >  2:834409/1‑150 (MQ=255)

                                                                                 aGGTATATCCCGAAGATTTAATCGGCAGGTTTTGGCTGATAATTTGCGACTGTTTTTTCGCCAGTAAACTGTCATTTTCATTATGGCTGGCACAGCAACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCa                                                                    >  2:1358979/1‑150 (MQ=255)

                                                                                  ggTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAt                                                                   >  1:541851/1‑150 (MQ=255)

                                                                                  ggTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAt                                                                   >  1:368867/1‑150 (MQ=255)

                                                                                     atatCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTg                                                                >  2:877229/1‑150 (MQ=255)

                                                                                       atCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtt                                                              <  1:877229/150‑1 (MQ=255)

                                                                                        tCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGttt                                                             <  2:368867/150‑1 (MQ=255)

                                                                                        tCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGttt                                                             <  2:937597/150‑1 (MQ=255)

                                                                                         cccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt                                                                  <  2:262576/144‑1 (MQ=255)

                                                                                         cccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAtt                                                                  >  1:262576/1‑144 (MQ=255)

                                                                                         cccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt                                                            >  1:1547640/1‑150 (MQ=255)

                                                                                         cccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt                                                            >  2:969/1‑150 (MQ=255)

                                                                                          ccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTa                                                           >  2:21066/1‑150 (MQ=255)

                                                                                                    aaTCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGt                                                 >  2:292881/1‑150 (MQ=255)

                                                                                                     aTCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTGCCTGGCATTTGTTTTATCGACATGGTa                                                <  2:898955/150‑1 (MQ=255)

                                                                                                     aTCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTa                                                <  2:1450732/150‑1 (MQ=255)

                                                                                                      tCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTaa                                               >  1:1408013/1‑150 (MQ=255)

                                                                                                        ggCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAaat                                             <  2:510635/150‑1 (MQ=255)

                                                                                                         gCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAaata                                            >  2:1716205/1‑150 (MQ=255)

                                                                                                          cAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAaataa                                           >  1:101900/1‑150 (MQ=255)

                                                                                                            ggTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAAtc                                         >  2:47666/1‑150 (MQ=255)

                                                                                                              ttttGTCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAAtctc                                       <  1:254170/150‑1 (MQ=255)

                                                                                                               tttGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATctct                                      <  1:1952814/150‑1 (MQ=255)

                                                                                                                ttGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTa                                     <  2:1190540/150‑1 (MQ=255)

                                                                                                                    cTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGt                                 >  1:1697134/1‑150 (MQ=255)

                                                                                                                    cTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGt                                 >  1:179066/1‑150 (MQ=255)

                                                                                                                    cTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGt                                 >  1:1370804/1‑150 (MQ=255)

                                                                                                                        tatTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGttatt                             <  1:826964/150‑1 (MQ=255)

                                                                                                                         atTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGttatta                            >  1:451198/1‑150 (MQ=255)

                                                                                                                             tGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATtaataa                        >  2:1707318/1‑150 (MQ=255)

                                                                                                                               cgcCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGttat                              >  2:451117/1‑142 (MQ=255)

                                                                                                                               cgcCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAaga                      <  2:1228862/150‑1 (MQ=255)

                                                                                                                                gcCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAagag                     >  1:970281/1‑150 (MQ=255)

                                                                                                                                 ccTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGt                    >  2:1749335/1‑150 (MQ=255)

                                                                                                                                   tGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTa                  >  1:907162/1‑150 (MQ=255)

                                                                                                                                   tGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTa                  <  1:21066/150‑1 (MQ=255)

                                                                                                                                     tgaaaatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAagag                     <  2:320964/138‑1 (MQ=255)

                                                                                                                                     tgaaaatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAagag                     >  1:320964/8‑145 (MQ=255)

                                                                                                                                      gaaaatgCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTTGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTaaaa               >  1:1470115/7‑150 (MQ=255)

                                                                                                                                       ttttCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAAt              <  2:184102/150‑1 (MQ=255)

                                                                                                                                       aaaatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAAt              >  1:1014262/6‑150 (MQ=255)

                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGt            >  1:433963/4‑150 (MQ=255)

                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGt            >  1:290837/4‑150 (MQ=255)

                                                                                                                                           tgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCa          >  1:816301/2‑150 (MQ=255)

                                                                                                                                           cGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAATTTTTTTAATAAGAGTTAAAATGTCa          <  2:1530757/150‑1 (MQ=255)

                                                                                                                                           cGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAATATGTCa          >  1:844075/1‑150 (MQ=255)

                                                                                                                                            tccAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCAc         <  1:1339896/149‑1 (MQ=255)

                                                                                                                                              cAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACtt       >  1:1713778/1‑150 (MQ=255)

                                                                                                                                                  aaaCTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGAt   >  1:1446214/1‑150 (MQ=255)

                                                                                                                                                   aaCTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATa  <  2:1547640/150‑1 (MQ=255)

CAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTTGATA  >  NC_000913/3113231‑3113527

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: `‑`