BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,127,469	G→A	G254G (GGC→GGT)	b2981 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,127,469	0	G	A	99.3%	1409.9 / NA	409	G254G (GGC→GGT)	b2981	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); major base A (216/189); minor base C (0/3); total (216/193)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 422 total aligned reads displayed.

GTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAT  >  NC_000913/3127379‑3127558

 ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                        <  2:6545701/93‑1 (MQ=255)

 ttCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCAccc                                                                                        <  2:1346244/93‑1 (MQ=255)

   cGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCAt                                                                                      <  1:5584025/93‑1 (MQ=255)

   cGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCAt                                                                                      <  2:8440387/93‑1 (MQ=255)

    ggCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATg                                                                                     <  1:3069389/93‑1 (MQ=255)

     gCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATgg                                                                                    <  2:6764756/93‑1 (MQ=255)

     gCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATgg                                                                                    <  1:2265417/93‑1 (MQ=255)

      cAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGt                                                                                   <  2:9109700/93‑1 (MQ=255)

        aGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcg                                                                                 >  1:7160170/1‑93 (MQ=255)

          ccGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTcgcg                                                                               >  1:9451008/1‑93 (MQ=255)

           cGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGa                                                                              >  1:2908893/1‑93 (MQ=255)

            gACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGaa                                                                             <  2:4460232/93‑1 (MQ=255)

            gACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGaa                                                                             >  2:3156711/1‑93 (MQ=255)

             aCAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAt                                                                            >  1:801270/1‑93 (MQ=255)

              cAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAAtt                                                                           >  2:4288403/1‑93 (MQ=255)

               aaTAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTc                                                                          >  2:7196046/1‑93 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                        <  2:10026926/93‑1 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                        >  1:5971001/1‑93 (MQ=255)

                 tAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCAt                                                                        <  2:9451803/93‑1 (MQ=255)

                  aaaCGCGCAGGGTGCAGAATCAATCTCAGCGGTATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATg                                                                       >  2:9752106/1‑93 (MQ=255)

                  aaaCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATg                                                                       >  1:7342155/1‑93 (MQ=255)

                   aaCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGt                                                                      <  1:1577309/93‑1 (MQ=255)

                    aCGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGtt                                                                     <  1:9358202/93‑1 (MQ=255)

                      gcgcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCg                                                                   <  2:7048460/93‑1 (MQ=255)

                        gcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGt                                                                 <  1:2027205/93‑1 (MQ=255)

                        gcAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGt                                                                 <  1:4174946/93‑1 (MQ=255)

                          aGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaa                                                               >  2:3681525/1‑93 (MQ=255)

                          aGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaa                                                               <  1:6655863/93‑1 (MQ=255)

                           gggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTaaa                                                              <  1:1690099/93‑1 (MQ=255)

                            ggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAAc                                                             >  1:9144193/1‑93 (MQ=255)

                            ggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAAc                                                             <  2:6362618/93‑1 (MQ=255)

                            ggTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAAc                                                             <  2:8255998/93‑1 (MQ=255)

                            ggTGCAGAATCAATCACAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAAc                                                             >  1:6975904/1‑93 (MQ=255)

                             gTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACg                                                            >  1:3971708/1‑93 (MQ=255)

                               gCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGc                                                          >  2:7990618/1‑93 (MQ=255)

                                 aGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGcac                                                        >  2:5252328/1‑93 (MQ=255)

                                   aatcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACaa                                                      >  1:5296550/1‑93 (MQ=255)

                                     tcaatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAAt                                                    >  2:5628688/1‑93 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                  >  2:4725172/1‑93 (MQ=255)

                                       aatcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATcc                                                  <  1:1425978/93‑1 (MQ=255)

                                        atcTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATccc                                                 <  1:8335162/93‑1 (MQ=255)

                                            cAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGt                                             <  2:1740491/93‑1 (MQ=255)

                                            cAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGt                                             <  1:9665721/93‑1 (MQ=255)

                                             aGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGtt                                            >  1:5140913/1‑93 (MQ=255)

                                               cGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgt                                          <  1:10007508/93‑1 (MQ=255)

                                                gATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgtg                                         >  2:78932/1‑93 (MQ=255)

                                                gATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTtgtg                                         >  1:6285498/1‑93 (MQ=255)

                                                  tataGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgct                                       >  2:4030102/1‑93 (MQ=255)

                                                    taGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgctgc                                     <  1:7852931/93‑1 (MQ=255)

                                                    taGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGtgctgc                                     >  2:2934065/1‑93 (MQ=255)

                                                     aGATCATATCTTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                    <  2:9208119/93‑1 (MQ=255)

                                                     aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                    <  1:8445071/93‑1 (MQ=255)

                                                     aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                    >  2:7395601/1‑93 (MQ=255)

                                                     aGATCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgcc                                    <  1:5875996/93‑1 (MQ=255)

                                                       aTCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccag                                  >  2:9067999/1‑93 (MQ=255)

                                                        tCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCCGCTGTGCTgccagc                                 >  2:5712392/1‑93 (MQ=255)

                                                        tCATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGGGCTgccagc                                 >  2:1145561/1‑93 (MQ=255)

                                                         cATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcc                                <  1:4775841/93‑1 (MQ=255)

                                                         cATATCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTgccagcc                                <  1:9246769/93‑1 (MQ=255)

                                                           tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGccagccag                              >  1:3884194/1‑93 (MQ=255)

                                                           tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGccagccag                              >  2:3073062/1‑93 (MQ=255)

                                                           tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGccagccag                              <  2:7628918/93‑1 (MQ=255)

                                                           tatCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTCCTGccagccag                              >  1:3673699/1‑93 (MQ=255)

                                                             tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            <  2:3358434/93‑1 (MQ=255)

                                                             tCGTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCagag                            <  1:8258629/93‑1 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                          <  2:3490064/93‑1 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                          >  1:2746602/1‑93 (MQ=255)

                                                               gTCCGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGtt                          >  1:4669061/1‑93 (MQ=255)

                                                                tCCGGCACCAGATATTTAAGTTGATCACCCATGGTGGCGAATTCCTGTTCGGGAAACGGCACAAATCCCCAGTTGTGCTCCCACCCATAGTTg                         >  1:3564260/1‑93 (MQ=255)

                                                                 ccGGCACCGGATATTTAAGTTGATCACCCTTGGTCGCGAGTTCATGTTCGGTAACCGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGa                        <  2:8960607/93‑1 (MQ=255)

                                                                 ccGGCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCACCCAGAGTTGa                        >  1:7477592/1‑93 (MQ=255)

                                                                    gCACCAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                     <  2:5061909/93‑1 (MQ=255)

                                                                    ccACCAGATATTTAAGTTGATCACCCATGGTCTCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGaaaa                     <  2:5469744/92‑1 (MQ=255)

                                                                      aCCAGATATTCAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtc                   <  1:6077474/93‑1 (MQ=255)

                                                                       ccAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtct                  >  2:9305218/1‑93 (MQ=255)

                                                                       ccAGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtct                  >  1:1246360/1‑93 (MQ=255)

                                                                         aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAAtctc                 <  1:3140180/92‑1 (MQ=255)

                                                                         aGATATTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCa                >  2:15951/1‑93 (MQ=255)

                                                                           atatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACg              >  1:6869068/1‑93 (MQ=255)

                                                                           atatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACg              <  2:2713343/93‑1 (MQ=255)

                                                                            tatTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGc             >  1:5519646/1‑93 (MQ=255)

                                                                             atTTAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCa            >  2:3869215/1‑93 (MQ=255)

                                                                              tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           <  1:538621/93‑1 (MQ=255)

                                                                              tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           <  2:6975846/93‑1 (MQ=255)

                                                                              tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           >  1:7039873/1‑93 (MQ=255)

                                                                              tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           >  2:9598298/1‑93 (MQ=255)

                                                                              tttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGGGCTGCCAGCCAGAGTTGAAAATCTCACGCAg           >  1:1617749/1‑93 (MQ=255)

                                                                               ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          <  2:2136359/93‑1 (MQ=255)

                                                                               ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          >  2:3570457/1‑93 (MQ=255)

                                                                               ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          >  1:4470260/1‑93 (MQ=255)

                                                                               ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          >  2:6060787/1‑93 (MQ=255)

                                                                               ttAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAgg          >  2:4938161/1‑93 (MQ=255)

                                                                                tAAGTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGa         <  1:8884076/93‑1 (MQ=255)

                                                                                   gTTGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCt      >  1:7643219/1‑93 (MQ=255)

                                                                                    ttGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTg     >  1:2312811/1‑93 (MQ=255)

                                                                                     tGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGc    >  2:2422292/1‑93 (MQ=255)

                                                                                     tGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGc    <  2:8816031/93‑1 (MQ=255)

                                                                                     tGATCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGc    <  1:3380664/93‑1 (MQ=255)

                                                                                       aTCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAt  <  2:1991669/93‑1 (MQ=255)

                                                                                       aTCACCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAt  >  1:3738679/1‑93 (MQ=255)

GTTCGGCAAGCCGACAATAAACGCGCAGGGTGCAGAATCAATCTCAGCGATATAGATCATATCGTCCGGCACCAGATATTTAAGTTGATCGCCCATGGTCGCGAATTCATGTTCGGTAAACGGCACAAATCCCCAGTTGTGCTGCCAGCCAGAGTTGAAAATCTCACGCAGGATCTGCAT  >  NC_000913/3127379‑3127558

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