BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	419,385:1	+T	39.0%	intergenic (+201/‑206)	phoR → / → brnQ	sensory histidine kinase in two‑component regulatory system with PhoB/branched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	419,383	1	.	T	39.0%	78.9 / 235.5	123	intergenic (+199/‑208)	phoR/brnQ	sensory histidine kinase in two‑component regulatory system with PhoB/branched‑chain amino acid transport system 2 carrier protein; LIV‑II transport system for Ile, Leu, and Val
Reads supporting (aligned to +/- strand): ref base . (34/41); new base T (24/24); total (58/65)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

Only 100 of 135 total aligned reads displayed.

TTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTG  >  NC_000913/419284‑419482

ttttGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCcctcc                                                                                                       >  2:922234/1‑100 (MQ=255)

ttttGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCcctcc                                                                                                       <  2:3382031/100‑1 (MQ=255)

 tttGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCcctcc‑t                                                                                                     <  1:1418623/100‑1 (MQ=255)

    gTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGc                                                                                                  >  2:870147/1‑100 (MQ=255)

      tCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑tttgct                                                                                                >  1:397623/1‑96 (MQ=255)

       cGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑tttgctt                                                                                               >  2:2483332/1‑95 (MQ=255)

       cGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑tttgctt                                                                                               <  2:1408938/100‑6 (MQ=255)

       cGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑tttgctt                                                                                               >  2:3436799/1‑95 (MQ=255)

       cGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCttt                                                                                               >  2:1895264/1‑100 (MQ=255)

       cGCATGATTAGCCATGTCTTTTTCACGAAAATAGTGTTTTATACTGGATGGTGATTTCTTATAGTTATATACCTCTGGTTTGTAGATCCCTCC‑TTGCttt                                                                                               >  1:2126437/1‑100 (MQ=255)

         cATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTa                                                                                              <  2:2510080/100‑1 (MQ=255)

         cATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTaa                                                                                             >  2:523484/1‑100 (MQ=255)

           tGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTaaa                                                                                            <  2:2233003/100‑1 (MQ=255)

           tGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTaaaa                                                                                           >  2:2933689/1‑100 (MQ=255)

           tGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTaaaa                                                                                           <  1:2657587/100‑1 (MQ=255)

            gATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAAc                                                                                          <  1:1586358/100‑1 (MQ=255)

             aTTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAAc                                                                                          >  1:1007899/1‑100 (MQ=255)

             aTTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACg                                                                                         >  2:49412/1‑100 (MQ=255)

              ttAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑t                                                                                       <  2:2558564/100‑1 (MQ=255)

               tAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑tt                                                                                      <  1:3451409/100‑1 (MQ=255)

               tAGCCATGTCCTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTGTGCTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑tt                                                                                      <  1:2275939/100‑1 (MQ=255)

                aGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑Tta                                                                                     <  1:1895248/100‑1 (MQ=255)

                 gCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑Tta                                                                                     <  1:1264664/100‑1 (MQ=255)

                  ccATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑Ttat                                                                                    <  1:414823/100‑1 (MQ=255)

                  ccATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑Ttata                                                                                   <  2:3174552/100‑1 (MQ=255)

                  ccATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑Ttata                                                                                   <  2:1587960/100‑1 (MQ=255)

                   cATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATaa                                                                                  >  2:1915966/1‑100 (MQ=255)

                    aTGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAg                                                                                 >  2:2457598/1‑100 (MQ=255)

                    aTGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAg                                                                                 <  1:598617/100‑1 (MQ=255)

                     tGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTCTAAAACG‑TTATAAGc                                                                                >  2:3402768/1‑100 (MQ=255)

                      gTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCg                                                                               <  1:193804/100‑1 (MQ=255)

                       tCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCg                                                                               <  1:675457/100‑1 (MQ=255)

                        cTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGt                                                                              <  2:2253707/100‑1 (MQ=255)

                        cTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGt                                                                              >  1:2910016/1‑100 (MQ=255)

                           tttCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTaa                                                                          <  1:3442498/100‑1 (MQ=255)

                            atctcGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTaaa                                                                         <  1:3364770/96‑1 (MQ=255)

                             tCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCTCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAAt                                                                        >  1:1969649/1‑100 (MQ=255)

                              cACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTCTACGCGTTTAAAt                                                                        >  2:2886282/1‑100 (MQ=255)

                              cACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAAt                                                                        >  1:1037055/1‑100 (MQ=255)

                              cACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACGCTTATAA‑CGTTTTAAtg                                                                       <  2:2333213/100‑2 (MQ=255)

                               aCGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATACGCGTTTAACTTg                                                                      >  1:1624290/1‑100 (MQ=255)

                               aCGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTg                                                                      <  2:1354675/100‑1 (MQ=255)

                                cGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTgc                                                                     <  2:2201496/100‑1 (MQ=255)

                                cGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTgc                                                                     <  1:1911430/100‑1 (MQ=255)

                                cGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTg                                                                      >  1:2862619/1‑99 (MQ=255)

                                 ggAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTgc                                                                     <  2:397626/100‑1 (MQ=255)

                                 ggAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    >  1:3255163/1‑100 (MQ=255)

                                  gAAATAGTGTTTTATACTGGTTTGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    <  1:3179876/100‑1 (MQ=255)

                                  gAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    <  2:3271092/100‑1 (MQ=255)

                                  gAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    <  1:2398213/100‑1 (MQ=255)

                                  gAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    <  1:2401842/100‑1 (MQ=255)

                                  gAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    >  2:710498/1‑100 (MQ=255)

                                  gAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTgcg                                                                    <  1:683665/99‑1 (MQ=255)

                                     aTAGTGTTTTATACTTGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTc                                                                <  2:3400352/100‑1 (MQ=255)

                                       aGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAg                                                              <  1:2658209/100‑1 (MQ=255)

                                        gtgtTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAgg                                                             <  1:2228531/100‑1 (MQ=255)

                                          gtTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGt                                                            >  2:2874459/1‑100 (MQ=255)

                                          gtTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGt                                                            >  2:3011135/1‑100 (MQ=255)

                                          gtTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGt                                                            >  1:3374500/1‑100 (MQ=255)

                                          gtTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAggag                                                           <  2:2345105/100‑3 (MQ=255)

                                            tttATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGc                                                          <  2:1007907/100‑1 (MQ=255)

                                             ttATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTg                                                        <  2:2135776/100‑1 (MQ=255)

                                               ataCTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGt                                                       <  2:26818/100‑1 (MQ=255)

                                               agaCTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTGAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGt                                                       >  2:2342027/3‑100 (MQ=255)

                                                taCTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTc                                                      >  2:1577310/1‑100 (MQ=255)

                                                taCTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCa                                                     <  2:3425429/100‑1 (MQ=255)

                                                 aCTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCAt                                                    >  1:405503/1‑100 (MQ=255)

                                                       tggtGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactga                                               >  1:3299112/1‑100 (MQ=255)

                                                       tggtGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgac                                              <  1:914097/100‑1 (MQ=255)

                                                        ggtGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgac                                              <  2:1506710/100‑1 (MQ=255)

                                                        ggtGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgact                                             >  1:1282393/1‑100 (MQ=255)

                                                         gtGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgact                                             >  2:776982/1‑100 (MQ=255)

                                                         gtGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgactg                                            <  2:331203/100‑1 (MQ=255)

                                                          tGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgactg                                            >  2:2761697/1‑100 (MQ=255)

                                                             tttCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCAtt                                        >  1:2667114/1‑100 (MQ=255)

                                                                   aTCGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAAcgc                                   >  1:2036147/1‑100 (MQ=255)

                                                                    tCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCgg                                 <  2:3432892/100‑1 (MQ=255)

                                                                    tCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCgg                                 <  1:1210259/100‑1 (MQ=255)

                                                                     cGCTATATACCTCTGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCgg                                 <  1:2373107/100‑1 (MQ=255)

                                                                         atataCCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAAt                            >  2:2775356/1‑100 (MQ=255)

                                                                          tataCCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATc                           <  2:3091736/100‑1 (MQ=255)

                                                                               ctctGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGaaa                       >  1:1672553/1‑100 (MQ=255)

                                                                               ctctGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATACGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAAACGcaa                       >  2:969776/1‑99 (MQ=255)

                                                                               ctctGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGaaaa                      >  2:1799561/1‑100 (MQ=255)

                                                                               ctctGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGaaaa                      >  2:2588414/1‑100 (MQ=255)

                                                                                 ctGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGaaaaa                     >  2:1054401/1‑100 (MQ=255)

                                                                                 ctGGTTTTTAGATCCCTCC‑TTGCTTTACAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAAc                    >  1:3293474/1‑100 (MQ=255)

                                                                                 ctGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAAc                    >  2:1345999/1‑100 (MQ=255)

                                                                                 ctGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAAc                    <  1:1530089/100‑1 (MQ=255)

                                                                                  tGGTTTTTAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAAc                    <  2:2591837/100‑1 (MQ=255)

                                                                                     tttttAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTAt                 >  1:1421863/1‑100 (MQ=255)

                                                                                       tttAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTAttct              <  1:1416708/100‑1 (MQ=255)

                                                                                        ttAGATCCCTCCTTTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTAttct              >  1:3021571/1‑100 (MQ=255)

                                                                                         tAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTAttcttc            <  1:49412/100‑1 (MQ=255)

                                                                                         tAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCAGTAAATCGAAAAACTAttcttc            <  1:1253005/100‑1 (MQ=255)

                                                                                           gATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTcgc          <  1:2235018/100‑1 (MQ=255)

                                                                                           gATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTcgc          >  2:1847576/1‑100 (MQ=255)

                                                                                            aTCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTcgcc         >  1:1943549/1‑100 (MQ=255)

                                                                                                 tcc‑tTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGAAAAAAAAAAAAAAACTATTCTTCGCCGCGcc    <  1:2945374/100‑1 (MQ=255)

                                                                                                   c‑tTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTg  <  1:3287163/100‑1 (MQ=255)

TTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCC‑TTGCTTTAAAACG‑TTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTG  >  NC_000913/419284‑419482

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