BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,857,430	A→G	I105T (ATC→ACC)	yidK ←	putative cotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,857,430	0	A	G	85.4%	709.0 / 108.2	302	I105T (ATC→ACC)	yidK	putative cotransporter
Reads supporting (aligned to +/- strand): ref base A (25/19); new base G (134/124); total (159/143)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 315 total aligned reads displayed.

TCGGCAGAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCAC  >  NC_000913/3857338‑3857521

tCGGCAGAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGa                                                                                             >  2:5991361/1‑93 (MQ=255)

  ggCAGAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGtg                                                                                           >  2:694291/1‑93 (MQ=255)

     aGAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGgtggtg                                                                                        >  2:3148084/1‑93 (MQ=255)

     aGAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGAtggtg                                                                                        <  2:7494792/93‑1 (MQ=255)

      gAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGtggtgg                                                                                       >  1:1547991/1‑93 (MQ=255)

      gAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGtggtgg                                                                                       >  2:6472452/1‑93 (MQ=255)

      gAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGtggtgg                                                                                       >  1:2695366/1‑93 (MQ=255)

      gAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGTAtggtgg                                                                                       >  1:2048063/1‑93 (MQ=255)

       aaaGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGc                                                                                      <  1:7469003/93‑1 (MQ=255)

        aaGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCa                                                                                     >  2:4141999/1‑93 (MQ=255)

         aGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCaa                                                                                    <  1:1654438/93‑1 (MQ=255)

            aGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGc                                                                                 >  2:2736077/1‑93 (MQ=255)

             gACGCCGGTGGCAATTAGGAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGcc                                                                                >  1:2208210/1‑93 (MQ=255)

                 ccGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCgcgc                                                                            <  1:941606/93‑1 (MQ=255)

                  cGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCt                                                                           >  1:6389855/1‑93 (MQ=255)

                   ggTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCtt                                                                          <  1:3190692/93‑1 (MQ=255)

                   ggTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCtt                                                                          >  2:1116692/1‑93 (MQ=255)

                    gTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTc                                                                         >  2:549175/1‑93 (MQ=255)

                     tGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCa                                                                        >  1:5895404/1‑93 (MQ=255)

                     tGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCa                                                                        <  2:200988/93‑1 (MQ=255)

                     tGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTCCCACCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCa                                                                        <  1:7404671/93‑1 (MQ=255)

                      ggCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAg                                                                       >  2:3819412/1‑93 (MQ=255)

                      ggCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAg                                                                       >  1:6907778/1‑93 (MQ=255)

                       gCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGa                                                                      >  1:2437006/1‑93 (MQ=255)

                       gCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGa                                                                      <  1:5471312/93‑1 (MQ=255)

                         aaTTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATa                                                                    >  2:320475/1‑93 (MQ=255)

                         aaTTAGGAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATa                                                                    >  1:1715347/1‑93 (MQ=255)

                         aaTTAGGAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATa                                                                    >  2:1775878/1‑93 (MQ=255)

                               gAAGCAGAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCgg                                                              <  1:3683633/93‑1 (MQ=255)

                                 aGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTa                                                            >  1:5099136/1‑93 (MQ=255)

                                  gCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTaa                                                           >  2:2452022/1‑93 (MQ=255)

                                  gCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTaa                                                           >  2:1361094/1‑93 (MQ=255)

                                     gAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTaaaaa                                                        >  2:4448989/1‑93 (MQ=255)

                                     gAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTaaaaa                                                        <  2:7111237/93‑1 (MQ=255)

                                     gAAGTCAATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTaaaaa                                                        >  2:370066/1‑93 (MQ=255)

                                      aaGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAg                                                       >  2:7315409/1‑93 (MQ=255)

                                       aGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGa                                                      <  2:4395580/93‑1 (MQ=255)

                                       aGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGa                                                      >  2:3542126/1‑93 (MQ=255)

                                        gTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGAt                                                     <  2:757817/93‑1 (MQ=255)

                                         tCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGAtt                                                    <  1:6821479/93‑1 (MQ=255)

                                         tCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGAtt                                                    <  2:5735989/93‑1 (MQ=255)

                                          caatgatAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTa                                                   <  1:3115437/91‑1 (MQ=255)

                                          cGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTa                                                   <  1:4351766/93‑1 (MQ=255)

                                          cGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTa                                                   <  1:469881/93‑1 (MQ=255)

                                           gatgatAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAg                                                  <  2:6357547/93‑1 (MQ=255)

                                            atgatAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAgc                                                 >  2:1640603/1‑93 (MQ=255)

                                            atgatAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAgc                                                 <  2:2682779/93‑1 (MQ=255)

                                              gatAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAgcgc                                               >  1:1239715/1‑93 (MQ=255)

                                                taataCGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGa                                             <  2:308812/93‑1 (MQ=255)

                                                taataCGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGa                                             >  2:6220618/1‑93 (MQ=255)

                                                  ataCGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAgg                                           >  1:5555425/1‑93 (MQ=255)

                                                  ataCGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAgg                                           >  2:2216274/1‑93 (MQ=255)

                                                   taCGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGa                                          >  1:3976824/1‑93 (MQ=255)

                                                    aCGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGaa                                         <  1:7362776/93‑1 (MQ=255)

                                                     cgcgTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAg                                        <  1:2634364/93‑1 (MQ=255)

                                                     cgcgTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAg                                        >  1:1035836/1‑93 (MQ=255)

                                                      gcgTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGa                                       <  2:3346082/93‑1 (MQ=255)

                                                      gcgTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGa                                       >  1:6748013/1‑93 (MQ=255)

                                                      gcgTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGa                                       >  2:4672947/1‑93 (MQ=255)

                                                      gcgTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGa                                       >  2:4034338/1‑93 (MQ=255)

                                                       cgtcgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGAt                                      <  1:6314015/93‑1 (MQ=255)

                                                        gtcgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGCGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATc                                     >  1:1859295/1‑93 (MQ=255)

                                                         tcgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCa                                    >  2:2099563/1‑93 (MQ=255)

                                                         tcgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCa                                    <  1:2886365/93‑1 (MQ=255)

                                                         tcgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCa                                    >  1:2982851/1‑93 (MQ=255)

                                                          cgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAg                                   <  1:4903027/93‑1 (MQ=255)

                                                          cgtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAg                                   <  2:5361472/93‑1 (MQ=255)

                                                           gtTTTATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGc                                  >  2:3899696/1‑93 (MQ=255)

                                                           gtTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGc                                  <  1:7134807/93‑1 (MQ=255)

                                                            ttttATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCg                                 >  1:5960094/1‑93 (MQ=255)

                                                            ttttATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCg                                 >  1:5687003/1‑93 (MQ=255)

                                                             tttATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCGATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGt                                >  1:6211501/1‑93 (MQ=255)

                                                              ttATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTc                               >  2:5210904/1‑93 (MQ=255)

                                                              ttATCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTc                               <  2:2884166/93‑1 (MQ=255)

                                                                aTCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCAc                             <  2:6718783/93‑1 (MQ=255)

                                                                 tCATAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCAcc                            <  1:5355795/93‑1 (MQ=255)

                                                                   aTAACGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGc                          >  1:282924/1‑93 (MQ=255)

                                                                   aTAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGc                          <  1:3301616/93‑1 (MQ=255)

                                                                       cGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGaa                      <  2:121169/93‑1 (MQ=255)

                                                                       cGTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGaa                      <  1:1469548/93‑1 (MQ=255)

                                                                        gTTCCTCCAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAg                     >  2:366074/1‑93 (MQ=255)

                                                                        gTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAg                     >  1:6470140/1‑93 (MQ=255)

                                                                          tcctccAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTc                   <  1:6596001/93‑1 (MQ=255)

                                                                           cctccAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCa                  >  2:4840075/1‑93 (MQ=255)

                                                                           cctccAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCa                  >  1:7261384/1‑93 (MQ=255)

                                                                           cctccAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCa                  <  2:2169373/93‑1 (MQ=255)

                                                                            ctccAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCAc                 >  2:6264808/1‑93 (MQ=255)

                                                                               cAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTc              <  1:3418709/93‑1 (MQ=255)

                                                                               cAGAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTc              <  1:2301184/93‑1 (MQ=255)

                                                                                 gAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTccc            >  2:5557927/1‑93 (MQ=255)

                                                                                 gAAAATCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTccc            >  2:1568529/1‑93 (MQ=255)

                                                                                  aaaaTCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCa           <  1:2540849/93‑1 (MQ=255)

                                                                                   aaaTCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAg          >  1:4265050/1‑93 (MQ=255)

                                                                                   aaaTCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAg          <  2:4300409/93‑1 (MQ=255)

                                                                                     aTCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGcc        <  2:4125902/93‑1 (MQ=255)

                                                                                     aTCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAATAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGcc        >  2:7071125/1‑93 (MQ=255)

                                                                                      tCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGccc       >  1:1160775/1‑93 (MQ=255)

                                                                                      tCGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGccc       >  2:4999286/1‑93 (MQ=255)

                                                                                       cGGGGGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCa      <  1:3573808/93‑1 (MQ=255)

                                                                                        ggggGTGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCAt     >  1:4997408/1‑93 (MQ=255)

TCGGCAGAAAGCAGACGCCGGTGGCAATTAGGAAGCAGAAGTCGATGATAATACGCGTCGTTTTATCATAACGTTCCTCCAGAAAATCGGGGATGGTGGCAATGCCGCGCTTCAGATAGCGCGGTAAAAAGATTAGCGCGAGGAAGATCAGCGTCACCGCTGAAGTCACTTCCCAGCCCATCAC  >  NC_000913/3857338‑3857521

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