BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	518,860	A→G	16.1%	F97S (TTC→TCC)	ybbO ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,860	0	A	G	16.1%	311.3 / 51.1	142	F97S (TTC→TCC)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (59/60); new base G (13/10); total (72/70)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 144 total aligned reads displayed.

ACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATC  >  NC_000913/518762‑518958

aCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGaa                                                                                                   <  2:909054/100‑1 (MQ=255)

aCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGaa                                                                                                   <  2:1162333/100‑1 (MQ=255)

 cAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGGGGAAAGGGGGCCATACATGCCGAAt                                                                                                  >  2:1893042/1‑100 (MQ=255)

  aGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATc                                                                                                 >  2:1632595/1‑100 (MQ=255)

     cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                              >  1:1938515/1‑100 (MQ=255)

     cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCgg                                                                                              <  1:2405905/100‑1 (MQ=255)

       cATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCa                                                                                            <  1:2881101/100‑1 (MQ=255)

         tGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtg                                                                                          <  1:1459935/100‑2 (MQ=255)

         tGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCAt                                                                                           >  2:1781774/1‑99 (MQ=255)

          ggTAAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtg                                                                                          <  1:2661881/99‑2 (MQ=255)

            tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtat                                                                                       <  2:213290/100‑5 (MQ=255)

            tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                       >  2:2450713/1‑100 (MQ=255)

            tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGtt                                                                                       >  1:2218693/1‑100 (MQ=255)

            tGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGt                                                                                        >  1:2218705/1‑99 (MQ=255)

             gAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTa                                                                                      >  2:1123762/1‑100 (MQ=255)

             gAGCTGGTGTGCCCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGGATTCGGCAtgtatg                                                                                      <  1:544644/100‑6 (MQ=255)

              tgCTGGTGTGCGCCGAATAAGTTGTCGGAAGACGGCGGTTCCGTCTGCGCACGGCTGAGGGTGGAAGGGGGGCCATACATGCCGGATTCGGTAtgtatgg                                                                                     <  1:1435035/99‑7 (MQ=255)

              aGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTaa                                                                                     >  1:367863/1‑100 (MQ=255)

                cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                   <  2:838859/100‑1 (MQ=255)

                cTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAg                                                                                   <  2:215791/100‑1 (MQ=255)

                 tGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGa                                                                                  <  1:497308/100‑1 (MQ=255)

                  ggTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGAt                                                                                 >  2:2891432/1‑100 (MQ=255)

                    tgtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATcc                                                                               <  1:658068/100‑1 (MQ=255)

                    tgtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATcc                                                                               <  1:251543/100‑1 (MQ=255)

                    tgtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATcc                                                                               <  1:237983/100‑1 (MQ=255)

                      tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                             >  2:1379105/1‑100 (MQ=255)

                      tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCa                                                                             >  2:2913797/1‑100 (MQ=255)

                        cgcCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGTCGAATCCGGCATTGTTAAAGATCCCATa                                                                           >  1:841606/1‑100 (MQ=255)

                         gacGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATa                                                                           <  1:1358358/97‑1 (MQ=255)

                          ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                         <  2:807735/100‑1 (MQ=255)

                          ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                         >  2:2693338/1‑100 (MQ=255)

                          ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACa                                                                         >  2:2693216/1‑100 (MQ=255)

                           aaaaaaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCAt                                                                            <  1:2792498/94‑1 (MQ=255)

                            gAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCAGAATCCGACATTGTGAAAGATCCCATACAGa                                                                       >  1:1701765/1‑100 (MQ=255)

                               aaaGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACaa                                                                    <  1:592675/100‑1 (MQ=255)

                                 aGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAAtt                                                                  >  1:2366715/1‑100 (MQ=255)

                                   ttGGCGGAAAACTGCTGTTCCATCTGCGCACGGCGGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTAt                                                                <  2:1290816/100‑1 (MQ=255)

                                    tGGCGGAACACTGCTGTTCCATCTGCGCACGGCTGAGGGTGGAAAGAGGGCCATACATGCCGAATCCGGCATTGTTAAAGCTCCCATACAGACAATTATc                                                               <  1:2058280/100‑1 (MQ=255)

                                      gCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCgg                                                             <  2:1671578/100‑1 (MQ=255)

                                      gCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCgg                                                             >  1:1681554/1‑100 (MQ=255)

                                       cGGAAAGCTCCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCg                                                              <  1:2564341/98‑1 (MQ=255)

                                       cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCgg                                                             >  1:107341/1‑99 (MQ=255)

                                       cGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAACCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGt                                                            >  1:1584898/1‑100 (MQ=255)

                                       cGGAAAACTGCTGTTCCATCTGCGCACGGCTGAGGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATGGTTAAAGATCCCATACAGACAATTATCGGt                                                            >  2:1937092/1‑100 (MQ=255)

                                        ggAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTc                                                           <  1:1278655/100‑1 (MQ=255)

                                           aaaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAgg                                                        <  2:1240621/100‑1 (MQ=255)

                                            aaCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAggg                                                       <  2:1748410/100‑1 (MQ=255)

                                             aCTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTCTCGGTCAGGGc                                                      >  1:2745335/1‑100 (MQ=255)

                                                gctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGAt                                                   >  1:358519/1‑100 (MQ=255)

                                                gctgTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGAt                                                   >  1:131405/1‑100 (MQ=255)

                                                 ctgTTCCAGCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCACACAGACAATTATCGGTCAGGGCGAtn                                                  <  1:1072818/100‑2 (MQ=255)

                                                   gTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                >  1:414611/1‑100 (MQ=255)

                                                   gTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAc                                                >  2:66084/1‑100 (MQ=255)

                                                    ttCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCAcc                                               <  2:597595/100‑1 (MQ=255)

                                                     tCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATAGATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCt                                              >  1:1800009/1‑100 (MQ=255)

                                                       cATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcg                                            <  2:2218718/100‑1 (MQ=255)

                                                       cATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcg                                            <  2:2218706/100‑1 (MQ=255)

                                                          cTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCtcgtcg                                         >  1:1279125/1‑100 (MQ=255)

                                                           tGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCgg                                        <  2:2362530/100‑1 (MQ=255)

                                                           tGCGCACGGCTGATGGTGGAAAGGGGGCCATAAATGCCGAATCCGGAATTGTTAAAGATCACATAAAGACAATTATCGGTAAGGGCGATAACCACGTCgg                                        >  1:116582/1‑100 (MQ=255)

                                                               cACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgc                                    >  2:2270131/1‑100 (MQ=255)

                                                                    cTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTgcgc                                  <  1:1901154/97‑1 (MQ=255)

                                                                     tGATGGTGGAAAGGGAGGCATACATGCAGAATCCGGCATTGTTAAAGATCCCATACCGACAATTATCGGTCAGGGTGATCACCTCGTCGTCTGCCCGATc                              <  1:2430737/100‑1 (MQ=255)

                                                                         ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAaca                          <  2:1567450/100‑1 (MQ=255)

                                                                         ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAaca                          >  2:389534/1‑100 (MQ=255)

                                                                         ggtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAaca                          >  2:311716/1‑100 (MQ=255)

                                                                          gtggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAacac                         >  2:1493899/1‑100 (MQ=255)

                                                                           tggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATTCAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACt                        >  2:1531942/1‑100 (MQ=255)

                                                                           tggAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACt                        >  2:2332338/1‑100 (MQ=255)

                                                                              aaaGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTc                     >  2:101723/1‑100 (MQ=255)

                                                                               aaGGGGGCCATACATGCCGATTCCGGCATTGTTAAAGATCCCATACAGTCAATTATCGGTCATTGCGATGATACCGCCGGCTCCGCGATCAACACTTTCt                    <  1:755416/100‑1 (MQ=255)

                                                                                aGGGGGCCATCCATGCCGAATCCGGCATTGTTAAAGAACCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTg                   <  2:2734019/100‑1 (MQ=255)

                                                                                 gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTgg                  <  2:1723144/100‑1 (MQ=255)

                                                                                 gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTgg                  <  1:867321/100‑1 (MQ=255)

                                                                                 gggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCAATCAACACTTTCTgg                  <  1:1379100/100‑1 (MQ=255)

                                                                                  ggggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                 <  1:2915491/100‑1 (MQ=255)

                                                                                  atctttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGt                 >  2:2441323/11‑100 (MQ=38)

                                                                                   gggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTg                <  2:583024/100‑1 (MQ=255)

                                                                                    ggCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGa               >  2:1056901/1‑100 (MQ=255)

                                                                                     tttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGaa              >  2:1197104/8‑100 (MQ=255)

                                                                                      ttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt             >  1:766783/7‑100 (MQ=255)

                                                                                      ttaacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt             >  2:39165/7‑100 (MQ=255)

                                                                                       taacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAt             >  1:2205026/6‑99 (MQ=255)

                                                                                        aacaaTGCCGGATTCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATcc           >  2:85204/5‑100 (MQ=255)

                                                                                         tACATGCCGAATCCGGCATTTTTAAAGATCCCATACAGACAATTATCTTTTTGTTTAATGACATCGTCAGCTGCGAGATCTACACTTTCTGGTGAATCCa          <  1:1923748/100‑1 (MQ=255)

                                                                                         tACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCa          >  2:521522/1‑100 (MQ=255)

                                                                                           cATGCCGAATTCTGCATTGTTAAAGATCCCATCCAGACAATGATCTTTCAGGTTGATGATCCGGCGGGCTCCGCGATCTACACTTTCTGGTGAATcc           <  1:1829962/97‑1 (MQ=255)

                                                                                           cATGCCGAATGCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTAAGCTGAGCGACCAACACTTTCTGGTGAATCCaaa        >  1:2095152/1‑98 (MQ=255)

                                                                                           cATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAga        >  2:2561165/1‑100 (MQ=255)

                                                                                           cATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAga        >  1:96067/1‑100 (MQ=255)

                                                                                            aTGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgat       >  1:2116175/1‑100 (MQ=255)

                                                                                             tGCCGAATCCGGTATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  1:970494/100‑1 (MQ=255)

                                                                                             tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATGACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  1:110940/100‑1 (MQ=255)

                                                                                             tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  2:2563876/100‑1 (MQ=255)

                                                                                             tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  2:1073123/100‑1 (MQ=255)

                                                                                             tGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatc      <  1:1139215/100‑1 (MQ=255)

                                                                                              gCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcg     <  2:320506/100‑1 (MQ=255)

                                                                                               ccGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcga    >  1:1871526/1‑100 (MQ=255)

                                                                                                 gAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc  <  2:686817/100‑1 (MQ=255)

                                                                                                 gAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAgatcgatc  <  1:1623926/100‑1 (MQ=255)

ACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATC  >  NC_000913/518762‑518958

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