BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,542,313	T→A	9.5%	W93R (TGG→AGG)	fimE →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,313	0	T	A	9.5%	627.9 / 44.6	167	W93R (TGG→AGG)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand): ref base T (74/77); new base A (13/3); total (87/80)
Fisher's exact test for biased strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

TAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCGGTATT  >  NC_000913/4542215‑4542411

taaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACt                                                                                                    <  1:5708/99‑1 (MQ=255)

taaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACt                                                                                                    <  2:2026363/99‑1 (MQ=255)

taaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACt                                                                                                    <  2:2307774/99‑1 (MQ=255)

 aaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTg                                                                                                   <  1:2184271/99‑1 (MQ=255)

 aaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTg                                                                                                   <  1:1620653/99‑1 (MQ=255)

  aataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTgg                                                                                                  >  2:894740/1‑99 (MQ=255)

   atatTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGa                                                                                                 >  2:1660414/1‑99 (MQ=255)

    tatTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGaa                                                                                                >  2:1213666/1‑99 (MQ=255)

    tatTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACAGGaa                                                                                                >  1:4031697/1‑99 (MQ=255)

    tatTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACAGGaa                                                                                                <  2:516329/99‑1 (MQ=255)

       tCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAgg                                                                                             >  2:3425562/1‑99 (MQ=255)

       tCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAgg                                                                                             <  1:181435/99‑1 (MQ=255)

       tCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAgg                                                                                             >  1:676644/1‑99 (MQ=255)

          ccGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGgcgc                                                                                          <  1:18635/99‑1 (MQ=255)

               tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGAcc                                                                                     >  1:3798978/1‑99 (MQ=255)

               tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGAcc                                                                                     <  2:3566760/99‑1 (MQ=255)

                gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCg                                                                                    <  1:2596820/99‑1 (MQ=255)

                  agaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGa                                                                                  <  1:1951629/99‑1 (MQ=255)

                  agaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGa                                                                                  >  1:763217/1‑99 (MQ=255)

                  agaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGa                                                                                  >  2:1414594/1‑99 (MQ=255)

                  agaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGa                                                                                  <  1:2407817/99‑1 (MQ=255)

                   gaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGAc                                                                                 <  1:1510264/99‑1 (MQ=255)

                   gaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGAc                                                                                 <  1:11101/99‑1 (MQ=255)

                   gaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGAc                                                                                 <  1:2168679/99‑1 (MQ=255)

                   gaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGAc                                                                                 >  2:3863257/1‑99 (MQ=255)

                    aaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACt                                                                                <  1:3138213/99‑1 (MQ=255)

                      cGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGa                                                                              >  2:749564/1‑99 (MQ=255)

                       ggATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGAc                                                                             <  2:2046017/99‑1 (MQ=255)

                       ggATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGAc                                                                             >  2:3284462/1‑99 (MQ=255)

                        gATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGACCGTGCTAACTGGAAAGGCGCTGACCGGACTGACg                                                                            <  2:1831290/99‑1 (MQ=255)

                          ttttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCt                                                                          >  2:1689433/1‑99 (MQ=255)

                           tttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCta                                                                         >  1:3145855/1‑99 (MQ=255)

                           tttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCta                                                                         <  1:3633767/99‑1 (MQ=255)

                            ttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtat                                                                        >  1:4213632/1‑99 (MQ=255)

                            ttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtat                                                                        >  1:1051568/1‑99 (MQ=255)

                             tCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACAGGAAAGGCGCTGACCGGACTGACGCtata                                                                       >  1:227719/1‑99 (MQ=255)

                              cTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTatat                                                                      >  2:18170/1‑99 (MQ=255)

                                aCCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtattt                                                                    >  2:3595606/1‑99 (MQ=255)

                                 ccGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtattta                                                                   >  2:749155/1‑99 (MQ=255)

                                  cGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatttat                                                                  <  1:1071082/99‑1 (MQ=255)

                                   gTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatttatt                                                                 <  2:3065897/99‑1 (MQ=255)

                                    ttCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATatttattt                                                                <  2:1946992/99‑1 (MQ=255)

                                    ttCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATatttattt                                                                <  1:719398/99‑1 (MQ=255)

                                     tCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtc                                                               <  2:1400216/99‑1 (MQ=255)

                                      cACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtct                                                              <  1:2301531/99‑1 (MQ=255)

                                       aCCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc                                                             >  2:4052060/1‑99 (MQ=255)

                                       aCCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc                                                             >  2:622203/1‑99 (MQ=255)

                                         ccGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgc                                                           >  1:1894767/1‑99 (MQ=255)

                                         ccGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgc                                                           <  1:3159170/99‑1 (MQ=255)

                                           gTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgccg                                                         >  2:489495/1‑99 (MQ=255)

                                            ttACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgccgc                                                        >  1:4069563/1‑99 (MQ=255)

                                            ttACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcgccgc                                                        <  1:4030332/99‑1 (MQ=255)

                                              aCGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCgg                                                      <  2:3899836/99‑1 (MQ=255)

                                                gTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGa                                                    >  1:904608/1‑99 (MQ=255)

                                                 ttttGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAg                                                   >  2:1786207/1‑99 (MQ=255)

                                                    tgatgaGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCg                                                >  2:3141385/1‑99 (MQ=255)

                                                     gatgaGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCgg                                               <  2:2295487/99‑1 (MQ=255)

                                                       tgaGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCt                                             >  1:3593267/1‑99 (MQ=255)

                                                         aGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACAGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCttt                                           >  1:2003996/1‑99 (MQ=255)

                                                          gCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTtc                                          >  1:1252651/1‑99 (MQ=255)

                                                          gCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTtc                                          <  2:1027628/99‑1 (MQ=255)

                                                            gTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTtctc                                        <  1:2886046/99‑1 (MQ=255)

                                                            gTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTtctc                                        >  2:786791/1‑99 (MQ=255)

                                                              gAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGc                                      <  1:900369/99‑1 (MQ=255)

                                                               aaGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGcc                                     >  2:4167158/1‑99 (MQ=255)

                                                               aaGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGcc                                     >  2:263138/1‑99 (MQ=255)

                                                               aaGCCGTGGAACGCTGGACCCAGGAACGTGCTAACAGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGcc                                     >  1:167235/1‑99 (MQ=255)

                                                                   cGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCcagca                                 <  2:1300793/99‑1 (MQ=255)

                                                                   cGTGGAACGCTGGACCCAGGAACGTGCTAACAGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCcagca                                 >  1:790744/1‑99 (MQ=255)

                                                                     tGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAgg                               >  1:3460328/1‑99 (MQ=255)

                                                                     tGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAgg                               >  1:3545512/1‑99 (MQ=255)

                                                                      ggAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGc                              <  2:2544821/99‑1 (MQ=255)

                                                                       gAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGcc                             >  1:52118/1‑99 (MQ=255)

                                                                         aCGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTa                           <  1:489757/99‑1 (MQ=255)

                                                                          cGCTGGACCCAGGAACGTGCTAACAGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTAt                          >  1:138688/1‑99 (MQ=255)

                                                                             tGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGc                       >  2:2324726/1‑99 (MQ=255)

                                                                               gACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCat                     >  2:2400936/1‑99 (MQ=255)

                                                                               gACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCat                     >  1:3855162/1‑99 (MQ=255)

                                                                                 cccAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCatta                   >  1:3205451/1‑99 (MQ=255)

                                                                                 cccAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCatta                   <  1:2398639/99‑1 (MQ=255)

                                                                                   cAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCattatt                 <  1:1488913/99‑1 (MQ=255)

                                                                                   cAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCattatt                 <  2:1553237/99‑1 (MQ=255)

                                                                                    aGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTc                <  1:3467137/99‑1 (MQ=255)

                                                                                     ggAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTcg               >  2:1021088/1‑99 (MQ=255)

                                                                                      gAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTcgc              <  1:1033916/99‑1 (MQ=255)

                                                                                        aCGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            >  1:3953135/1‑99 (MQ=255)

                                                                                        aCGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            <  2:949976/99‑1 (MQ=255)

                                                                                        aCGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            >  1:3075634/1‑99 (MQ=255)

                                                                                        aCGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            <  2:3665722/99‑1 (MQ=255)

                                                                                        aCGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            >  2:2970191/1‑99 (MQ=255)

                                                                                        aCGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            <  2:165025/99‑1 (MQ=255)

                                                                                        aCGTGCTAACAGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGa            >  1:3010157/1‑99 (MQ=255)

                                                                                          gTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATg          >  2:2700308/1‑99 (MQ=255)

                                                                                           tGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGc         >  1:2989867/1‑99 (MQ=255)

                                                                                             cTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCg       <  1:2111169/99‑1 (MQ=255)

                                                                                             cTAACAGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCg       >  2:3906045/1‑99 (MQ=255)

                                                                                                 cTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCGGTAt   <  1:3181491/99‑1 (MQ=255)

                                                                                                  tGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCGGTAtt  >  2:2823086/1‑99 (MQ=255)

                                                                                                  tGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCGGTAtt  <  2:546223/99‑1 (MQ=255)

                                                                                                  tGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCGGTAtt  >  1:2953084/1‑99 (MQ=255)

TAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGCTTTCTCGCCAGCAGGCCTATCGCATTATTCGCGATGCCGGTATT  >  NC_000913/4542215‑4542411

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