BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	591,967	C→T	M404I (ATG→ATA)	nfrB ←	bacteriophage N4 receptor, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	591,967	0	C	T	100.0%	1017.8 / NA	249	M404I (ATG→ATA)	nfrB	bacteriophage N4 receptor, inner membrane subunit
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (121/128); total (121/128)

Only 100 of 253 total aligned reads displayed.

GAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCA  >  NC_000913/591869‑592064

 aaaGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTAtt                                                                                                >  1:786465/1‑101 (MQ=255)

 aaaGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTAtt                                                                                                >  2:786466/1‑101 (MQ=255)

     cagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAAcc                                                                                            >  1:786469/1‑101 (MQ=255)

      agcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCa                                                                                           >  2:786472/1‑101 (MQ=255)

       gcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCAt                                                                                          >  1:786474/1‑101 (MQ=255)

          ccTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATgcg                                                                                       >  1:786477/1‑101 (MQ=255)

           cTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATgcgc                                                                                      <  1:786160/101‑1 (MQ=255)

           cTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATgcgc                                                                                      >  2:786478/1‑101 (MQ=255)

            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCAt                                                                                          <  2:786242/96‑1 (MQ=255)

            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCt                                                                                     <  1:786062/101‑1 (MQ=255)

               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATgc                                                                                        <  2:786308/95‑1 (MQ=255)

               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGcc                                                                                  >  2:786483/1‑101 (MQ=255)

                 cAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGc                                                                                <  1:786099/101‑1 (MQ=255)

                  aGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCt                                                                               <  1:786183/101‑1 (MQ=255)

                    gCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGa                                                                             <  1:786145/101‑1 (MQ=255)

                    gCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGa                                                                             >  2:786487/1‑101 (MQ=255)

                     ccGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGaa                                                                            <  2:785995/101‑1 (MQ=255)

                      cGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTg                                                                                    <  2:786115/92‑1 (MQ=255)

                      cGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGaaa                                                                           >  1:786490/1‑101 (MQ=255)

                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGaaaa                                                                          <  1:786177/101‑1 (MQ=255)

                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGaaaa                                                                          <  2:786029/101‑1 (MQ=255)

                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGaaaa                                                                          <  1:786090/101‑1 (MQ=255)

                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGaaaa                                                                          <  1:786016/101‑1 (MQ=255)

                        tagtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAAt                                                                         >  2:786492/1‑101 (MQ=255)

                        tagtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAAt                                                                         <  2:786018/101‑1 (MQ=255)

                         agtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATa                                                                        >  1:786495/1‑101 (MQ=255)

                          gtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTg                                                                              <  2:786068/94‑1 (MQ=255)

                           tagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGa                                                                      >  2:786497/1‑101 (MQ=255)

                             gCCAGTAACGAAAATCCCCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGaaa                                                                    <  2:786249/101‑1 (MQ=255)

                              ccAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAg                                                                   <  1:786198/101‑1 (MQ=255)

                               cAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAg                                                                   >  1:786502/1‑100 (MQ=255)

                                 gTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGaa                                                                <  2:786126/101‑1 (MQ=255)

                                  tAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGaaa                                                               <  1:785914/101‑1 (MQ=255)

                                   aaCGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAAt                                                              >  2:786503/1‑101 (MQ=255)

                                      gAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGcc                                                           <  1:785968/101‑1 (MQ=255)

                                          aTCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGc                                                       <  2:786195/101‑1 (MQ=255)

                                          aTCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGc                                                       <  1:785885/101‑1 (MQ=255)

                                           tCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCa                                                      >  2:786508/1‑101 (MQ=255)

                                           tCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCa                                                      >  1:786507/1‑101 (MQ=255)

                                            cccccGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGaaa                                                                    <  2:786154/84‑1 (MQ=255)

                                            cACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGc                                                                                   <  1:786188/71‑1 (MQ=255)

                                            cACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                                                     <  2:786072/101‑1 (MQ=255)

                                             aCCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATc                                                    >  2:786510/1‑101 (MQ=255)

                                              cccGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                                                     <  2:786217/99‑1 (MQ=255)

                                              cccGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATc                                                    <  2:786207/100‑1 (MQ=255)

                                                cGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg                                                 >  1:786516/1‑101 (MQ=255)

                                                cGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg                                                 >  2:786517/1‑101 (MQ=255)

                                                 gCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGc                                                >  2:786520/1‑101 (MQ=255)

                                                 gCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGc                                                <  1:786104/101‑1 (MQ=255)

                                                  cTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGcc                                               >  1:786521/1‑101 (MQ=255)

                                                   tGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCca                                              <  1:785837/101‑1 (MQ=255)

                                                    gCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcac                                             >  2:786525/1‑101 (MQ=255)

                                                    gCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcac                                             <  2:786092/101‑1 (MQ=255)

                                                     cACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcaca                                            >  2:786527/1‑101 (MQ=255)

                                                     cACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcaca                                            <  1:786311/101‑1 (MQ=255)

                                                        gATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAAc                                         >  1:786529/1‑101 (MQ=255)

                                                         aTGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGATAGGAAATGCCAGGCATCGGGCCACAAACt                                        >  1:786532/1‑101 (MQ=255)

                                                         aTGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACt                                        <  2:786136/101‑1 (MQ=255)

                                                           gCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACttt                                      <  2:786291/101‑1 (MQ=255)

                                                              gTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCAt                                   >  1:786536/1‑101 (MQ=255)

                                                               ttAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcaca                                            >  2:786540/1‑91 (MQ=255)

                                                               ttAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATa                                  <  1:786025/101‑1 (MQ=255)

                                                               ttAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATa                                  >  2:786538/1‑101 (MQ=255)

                                                                tAACCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATAc                                 >  1:786542/1‑101 (MQ=255)

                                                                 aaCCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACg                                <  1:786298/101‑1 (MQ=255)

                                                                 aaCCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACg                                >  2:786544/1‑101 (MQ=255)

                                                                  aCCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCt                                                                                     >  2:786547/1‑47 (MQ=255)

                                                                  aCCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGc                               <  2:786162/101‑1 (MQ=255)

                                                                  aCCATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGc                               <  1:786194/101‑1 (MQ=255)

                                                                    cATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCca                             >  1:786552/1‑101 (MQ=255)

                                                                    cATCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCca                             >  2:786551/1‑101 (MQ=255)

                                                                     aTCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcag                            <  1:786245/101‑1 (MQ=255)

                                                                      tCAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg                                                 <  1:786256/79‑1 (MQ=255)

                                                                       cAAACCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcagca                          >  1:786554/1‑101 (MQ=255)

                                                                        aaaCCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcagcaa                         <  1:786156/101‑1 (MQ=255)

                                                                          aCCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcagcaaca                       >  1:786560/1‑101 (MQ=255)

                                                                          aCCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcagcaaca                       >  1:786558/1‑101 (MQ=255)

                                                                           ccAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAAt                                                                         <  2:785862/50‑1 (MQ=255)

                                                                           ccAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcagcaacag                      <  1:786227/101‑1 (MQ=255)

                                                                             aaaGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCcagcaacagca                    >  2:786562/1‑101 (MQ=255)

                                                                               aGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCaaa                  >  1:786564/1‑101 (MQ=255)

                                                                                gTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCaaaa                 <  1:786280/101‑1 (MQ=255)

                                                                                 tttAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAg                <  1:786130/101‑1 (MQ=255)

                                                                                  ttAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGc               >  2:786568/1‑101 (MQ=255)

                                                                                  ttAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGc               >  1:786569/1‑101 (MQ=255)

                                                                                  ttAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGc               <  1:786209/101‑1 (MQ=255)

                                                                                   tAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCt              <  1:786497/101‑1 (MQ=255)

                                                                                   tAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCt              <  1:786325/101‑1 (MQ=255)

                                                                                    aGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTg             <  2:786003/101‑1 (MQ=255)

                                                                                      ccACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGa           >  2:786574/1‑101 (MQ=255)

                                                                                      ccACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGa           >  1:786573/1‑101 (MQ=255)

                                                                                      ccACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGa           <  2:786240/101‑1 (MQ=255)

                                                                                          agcagGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGatcat       >  1:786579/1‑101 (MQ=255)

                                                                                           gcagGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGatcatc      >  1:786581/1‑101 (MQ=255)

                                                                                            cagGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGAtcatca     >  2:786584/1‑101 (MQ=255)

                                                                                             agGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCAc    >  1:786586/1‑101 (MQ=255)

                                                                                              gggTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGatca        <  2:786327/96‑1 (MQ=255)

                                                                                              gggTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCAcc   <  2:786388/101‑1 (MQ=255)

                                                                                               ggTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCa  <  1:786075/101‑1 (MQ=255)

                                                                                               ggTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCa  <  1:786053/101‑1 (MQ=255)

GAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATACGCCAGCAACAGCAAAAGCTGGATCATCACCA  >  NC_000913/591869‑592064

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