BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,524	C→A	intergenic (+391/+166)	yjhV → / ← fecE	pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,524	0	C	A	95.7%	609.7 / 8.2	187	intergenic (+391/+166)	yjhV/fecE	pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand): ref base C (2/6); new base A (81/98); total (83/104)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

TGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCA  >  NC_000913/4510428‑4510619

tGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGatc                                                                                               <  1:1639863/99‑1 (MQ=255)

  tcGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTg                                                                                             <  1:407671/98‑1 (MQ=255)

  ccGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGACCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTg                                                                                             <  1:357182/99‑1 (MQ=255)

   cGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCGAGACTGTTTTGAAGGCCGCGCCCAAGCCTGTACAGATCTgg                                                                                            >  2:4080623/1‑99 (MQ=255)

    gTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGt                                                                                           <  1:2092397/99‑1 (MQ=255)

     tAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGtt                                                                                          <  1:302085/99‑1 (MQ=255)

      aaGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGttt                                                                                         <  2:3713856/99‑1 (MQ=255)

      aaGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGttt                                                                                         >  2:1401645/1‑99 (MQ=255)

       aGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGACGCGCCAAAGCCTGTACAGATCTGGTTTg                                                                                        <  2:1925726/99‑1 (MQ=255)

           ttATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGtt                                                                                    <  2:2664764/99‑1 (MQ=255)

             ataCAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGa                                                                                  <  1:3756226/99‑1 (MQ=255)

             ataCAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGa                                                                                  >  2:2731755/1‑99 (MQ=255)

              taCAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGAt                                                                                 <  1:3103766/99‑1 (MQ=255)

              taCAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGAt                                                                                 >  2:3683647/1‑99 (MQ=255)

              taCAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGAt                                                                                 <  2:2149365/99‑1 (MQ=255)

               aCAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGAtt                                                                                >  2:4007840/1‑99 (MQ=255)

                cAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTa                                                                               <  2:1517289/99‑1 (MQ=255)

                 aGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTAt                                                                              >  2:641684/1‑99 (MQ=255)

                   aaaTGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGa                                                                            >  2:2599091/1‑99 (MQ=255)

                     ttGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGTTCACGACTGGTTTGAAAGCCGCGCCCATGCCTGTACAGATCTGGTTTGCGTTGATTATGAAc                                                                          <  2:3609276/98‑1 (MQ=255)

                     aTGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAAc                                                                          <  2:3900720/99‑1 (MQ=255)

                     aTGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAAc                                                                          <  1:3523058/99‑1 (MQ=255)

                       gTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCACGCCTGTACAGATCTGGTTGGCGTTGATTATGAACCt                                                                        >  2:3015508/1‑99 (MQ=255)

                       gTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCt                                                                        >  2:2875491/1‑99 (MQ=255)

                          cGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTc                                                                     >  2:105336/1‑99 (MQ=255)

                           gTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCa                                                                    <  1:2366655/99‑1 (MQ=255)

                           gTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCa                                                                    <  1:1845341/99‑1 (MQ=255)

                            ttAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAg                                                                   >  2:1169715/1‑99 (MQ=255)

                             tAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGc                                                                  >  2:2944139/1‑99 (MQ=255)

                             aaaGCAGAGTGCAAAATTGACTACCGTGATTACGAAACGGAAGCAAGCCGAGCCCAAGCCTGTACAGATCTGGTTTGCATTGATTAAGAACCTGTCAGc                                                                  <  2:817946/98‑1 (MQ=255)

                               aGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCt                                                                >  1:294705/1‑99 (MQ=255)

                                gCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTa                                                               <  2:2703211/99‑1 (MQ=255)

                                gCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTa                                                               >  1:1391936/1‑99 (MQ=255)

                                 cAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTaa                                                              >  2:2182552/1‑99 (MQ=255)

                                  agagTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTaaa                                                             <  1:1895067/99‑1 (MQ=255)

                                   gagTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaag                                                            <  1:996404/99‑1 (MQ=255)

                                   gagTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaag                                                            >  1:106114/1‑99 (MQ=255)

                                      ttCAAAATTGATTGCAGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCCGTACAGCTCTGGTTTGCGTTGATTATGAACCAGTCAGCCTAaagcaa                                                         >  2:56433/1‑99 (MQ=255)

                                        cAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaagc                                                       <  1:1319741/99‑1 (MQ=255)

                                        cAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaagc                                                       <  1:138513/99‑1 (MQ=255)

                                        cAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaagc                                                       <  1:3054116/99‑1 (MQ=255)

                                        cAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaagc                                                       >  1:283122/1‑99 (MQ=255)

                                        cAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaagc                                                       >  2:1924901/1‑99 (MQ=255)

                                        cAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaagc                                                       <  2:1465214/99‑1 (MQ=255)

                                          aaaTTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCgg                                                     <  1:3875835/99‑1 (MQ=255)

                                          aaaTTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCgg                                                     >  2:17783/1‑99 (MQ=255)

                                             ttgattgCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATg                                                  <  2:2214801/99‑1 (MQ=255)

                                               gattgCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGa                                                <  1:2550311/99‑1 (MQ=255)

                                                attgCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGAc                                               <  1:1533256/99‑1 (MQ=255)

                                                  tgCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGa                                             <  1:3177464/99‑1 (MQ=255)

                                                   gCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGAt                                            <  2:2095460/99‑1 (MQ=255)

                                                   gCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGAt                                            <  1:605320/99‑1 (MQ=255)

                                                   gCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGAt                                            <  2:2250830/99‑1 (MQ=255)

                                                   gCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGAt                                            >  2:3626536/1‑99 (MQ=255)

                                                   gCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGAt                                            <  2:1169429/99‑1 (MQ=255)

                                                    ccGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATg                                           >  1:1991483/1‑99 (MQ=255)

                                                    ccGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATg                                           >  2:3355978/1‑99 (MQ=255)

                                                    ccGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATg                                           >  1:2297906/1‑99 (MQ=255)

                                                     cGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGa                                          <  1:1776450/99‑1 (MQ=255)

                                                      gTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAg                                         >  1:58799/1‑99 (MQ=255)

                                                       tGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGt                                        >  1:43741/1‑99 (MQ=255)

                                                       tGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGt                                        >  2:3516029/1‑99 (MQ=255)

                                                       tGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGt                                        >  2:2435516/1‑99 (MQ=255)

                                                        gATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTa                                       >  1:1059621/1‑99 (MQ=255)

                                                        gATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCGAAAGCAAGCGGATGGACGATGAGTa                                       >  1:275652/1‑99 (MQ=255)

                                                           cACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTg                                    <  1:2627216/99‑1 (MQ=255)

                                                            aCGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTgg                                   <  1:1920085/99‑1 (MQ=255)

                                                             cGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGt                                  <  2:1804261/99‑1 (MQ=255)

                                                              gACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTa                                 >  1:2700062/1‑99 (MQ=255)

                                                               aCTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTaa                                <  1:1423639/99‑1 (MQ=255)

                                                               aCTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTaa                                <  2:3496983/99‑1 (MQ=255)

                                                                 tGGTTTGAAAGTCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATc                              >  1:3203209/1‑99 (MQ=255)

                                                                   gTTTGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCtt                            >  1:2204181/1‑99 (MQ=255)

                                                                    tttGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCttt                           >  1:759865/1‑99 (MQ=255)

                                                                    tttGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCttt                           <  2:3447298/99‑1 (MQ=255)

                                                                    tttGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCttt                           >  1:3666854/1‑99 (MQ=255)

                                                                    tgtGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCttt                           <  2:1611855/97‑1 (MQ=255)

                                                                     ttGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTc                          >  1:4166071/1‑99 (MQ=255)

                                                                      tGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCa                         <  1:258890/99‑1 (MQ=255)

                                                                      tGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCa                         <  2:3748743/99‑1 (MQ=255)

                                                                      tGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCa                         <  1:2436249/99‑1 (MQ=255)

                                                                      tGAAAGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCa                         <  1:949622/99‑1 (MQ=255)

                                                                          aGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGt                     >  2:2763851/1‑99 (MQ=255)

                                                                          aGCCGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGt                     >  1:1705821/1‑99 (MQ=255)

                                                                            ccGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTcc                   >  1:4193958/1‑99 (MQ=255)

                                                                            ccGCGCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTcc                   >  1:3716637/1‑99 (MQ=255)

                                                                             cgcgCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCg                  >  1:4228476/1‑99 (MQ=255)

                                                                             cgcgCCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCg                  >  1:1707098/1‑99 (MQ=255)

                                                                              gcgcCCAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCgg                 >  1:35893/1‑99 (MQ=255)

                                                                                  ccAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGaaaa             <  2:2288788/99‑1 (MQ=255)

                                                                                  ccAAGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGaaaa             <  2:321643/99‑1 (MQ=255)

                                                                                     aGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGtt          <  2:3250919/99‑1 (MQ=255)

                                                                                     aGCCTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGtt          <  2:1850107/99‑1 (MQ=255)

                                                                                      gcatGTACAGATCTGGTTTGGGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGACAAGTTc         >  2:3586904/4‑99 (MQ=255)

                                                                                        cTGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAg       <  2:972685/99‑1 (MQ=255)

                                                                                         tGTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGc      <  2:2846764/99‑1 (MQ=255)

                                                                                          gTACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGcc     >  1:111122/1‑99 (MQ=255)

                                                                                           tACAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGccc    >  1:4012236/1‑99 (MQ=255)

                                                                                             cAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCa  >  1:289297/1‑99 (MQ=255)

                                                                                             cAGATCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCa  >  1:2310370/1‑99 (MQ=255)

TGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCA  >  NC_000913/4510428‑4510619

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