BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,118,628	A→G	5.6%	F140S (TTT→TCT)	yceJ ←	putative cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,118,628	0	A	G	5.6%	325.9 / 6.8	108	F140S (TTT→TCT)	yceJ	putative cytochrome b561
Reads supporting (aligned to +/- strand): ref base A (41/60); new base G (3/3); total (45/63)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 108 total aligned reads displayed.

TATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTC  >  NC_000913/1118481‑1118775

tatcCTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAg                                                                                                                                                   <  2:675455/150‑1 (MQ=255)

 ataCTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGt                                                                                                                                                  >  1:626184/1‑150 (MQ=255)

 ataCTCCATAGGCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTAGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGt                                                                                                                                                  <  2:1417788/150‑1 (MQ=255)

    cTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTg                                                                                                                                               <  2:1066267/150‑1 (MQ=255)

    cTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTg                                                                                                                                               >  1:580182/1‑150 (MQ=255)

    cTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTg                                                                                                                                               >  1:304298/1‑150 (MQ=255)

      ccATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                                                                                             <  1:1600241/150‑1 (MQ=255)

      ccATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                                                                                             <  2:117094/150‑1 (MQ=255)

       cATAGTCAGATGAGGACTTTCCCAGCATGCGCTTCAGAGTGTCGGCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCt                                                                                                                                            <  1:1209508/150‑1 (MQ=255)

       cATAGTAAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCt                                                                                                                                            <  2:1255041/150‑1 (MQ=255)

           ctcccacgacgacTTTCCCACCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCGTAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgc                                                                                                                                        <  2:1665673/143‑1 (MQ=255)

            tCAGATGACGACTTTCCCAGCTTGCGCTTCAGAGTGGCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTgcg                                                                                                                                       <  1:327213/150‑1 (MQ=255)

             cAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgccc                                                                                                                                      <  1:532774/150‑3 (MQ=255)

             cAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACCCTCCCCGCGAGCCAAAAATGCCGGGCACCGGCAAAGTCTGCCTGTgcgc                                                                                                                                      >  1:516018/1‑150 (MQ=255)

               gATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACAGGCAGACTTTGCCGGTgccctg                                                                                                                                    <  1:1339413/150‑5 (MQ=255)

                aTGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAg                                                                                                                                   >  2:802781/1‑150 (MQ=255)

                 tGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCCGGGCACCGGCAAAGTCTGCCTGTGCGCCAGc                                                                                                                                  >  1:1375482/1‑150 (MQ=255)

                       cTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTGCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGc                                                                                                                            <  1:595601/150‑1 (MQ=255)

                       atttCCGAGCATGCGCTTCAGAGTTTCGTCATTAGAGATGAAATGATGCTTCAGGGCCATAAATCCGTTCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAATATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGc                                                                                                                            <  1:1785887/149‑1 (MQ=255)

                          tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                         <  1:1372449/150‑1 (MQ=255)

                          tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                         <  1:70412/150‑1 (MQ=255)

                          tCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGAAGAAATTATGCTACAGGGCCATAAATACGTGCATAACGGACAGTACTACGACCCTCCACGCGAGCCAAAAATGCAGGTCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                         >  1:876044/1‑150 (MQ=255)

                          tCCCAGCATGCGCTTCAGAGTGTCGTCTTAATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATACCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAg                                                                                                                         <  1:1167620/150‑1 (MQ=255)

                           cccAGCATGCGCTTCAGAGTGGCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAgg                                                                                                                        <  1:1025256/150‑1 (MQ=255)

                            ccAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGAACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAggg                                                                                                                       <  1:868167/150‑1 (MQ=255)

                              aGCATGAGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTc                                                                                                                     <  1:1737998/150‑1 (MQ=255)

                                  tGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCgg                                                                                                                 <  1:406260/150‑1 (MQ=255)

                                   gcgcTTCAGAGTGTCGTCTTTATCGTTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCggg                                                                                                                <  1:1107031/150‑1 (MQ=255)

                                      cTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGAc                                                                                                             >  2:737128/1‑150 (MQ=255)

                                         cAGAGTGTCGTCTTTATCGATGAAATGATGCTGCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTc                                                                                                          <  1:537241/150‑1 (MQ=255)

                                          agagTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATACAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCa                                                                                                         <  2:176113/150‑1 (MQ=255)

                                          agagTGTCGTCTTTATCGATGAAATGATGATTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCa                                                                                                         >  1:1289983/1‑150 (MQ=255)

                                              tgtCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAAcc                                                                                                     <  2:1825223/150‑1 (MQ=255)

                                               gtcgtcTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCCAACCa                                                                                                    >  1:8640/1‑150 (MQ=255)

                                                        aTCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                                                                                           <  2:899417/150‑1 (MQ=255)

                                                        aTCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGAACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                                                                                           <  1:504091/150‑1 (MQ=255)

                                                          cGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTATGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGc                                                                                         >  2:1078172/1‑150 (MQ=255)

                                                          cGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGGCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGc                                                                                         <  2:1174220/150‑1 (MQ=255)

                                                           gATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCt                                                                                        >  1:627934/1‑150 (MQ=255)

                                                            aTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                       <  2:227059/150‑1 (MQ=255)

                                                            aTGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGCCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                       <  2:896513/150‑1 (MQ=255)

                                                            aTGAAATGATGCTTCAGGGCCATAAATCCGAGCATACCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTg                                                                                       <  1:994925/150‑1 (MQ=255)

                                                             tGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGCGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGACAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGa                                                                                      >  2:1653294/1‑150 (MQ=255)

                                                             tGAAATGATGCTTCAGGGCCATAAATCAGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCCAAAACGCTGa                                                                                      >  1:892046/1‑150 (MQ=255)

                                                                  tgatgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGt                                                                                 >  2:790865/1‑150 (MQ=255)

                                                                    atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                               <  1:1613983/150‑1 (MQ=255)

                                                                    atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                               <  2:408283/150‑1 (MQ=255)

                                                                    atgCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGttt                                                                               >  2:499648/1‑150 (MQ=255)

                                                                       cTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTAcc                                                                            >  2:845728/1‑150 (MQ=255)

                                                                       cTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTATAcc                                                                            >  1:421293/1‑150 (MQ=255)

                                                                             ggCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGc                                                                      <  2:1294785/150‑1 (MQ=255)

                                                                              gCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGCCGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCa                                                                     >  1:1417209/1‑150 (MQ=255)

                                                                                 aTAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCCGCCAAAAACGCTGCTCGGTTTACCATCGGCAGtt                                                                  >  1:436274/1‑150 (MQ=255)

                                                                                  tAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTg                                                                 <  2:626184/150‑1 (MQ=255)

                                                                                  tAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACACTGATCGGTTTACCATCGGCAGTTg                                                                 <  2:1354108/150‑1 (MQ=255)

                                                                                   aaaTCCGTGCATACCGGACAGTAGTACGCCGCTCCACGAGAGCCAAAAATTCAGGGCACCTGCAAAGTCTGCCGGTGCGCCAGCGGCGGCAAGGGTCGCGGTGACGTCAAACCAGCAAAAAACGCTGATCGGTTTACCATCGGCAGTTga                                                                <  2:287301/150‑1 (MQ=255)

                                                                                     aTCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTgaga                                                              <  1:396689/150‑1 (MQ=255)

                                                                                       ccGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATc                                                            >  1:30143/1‑150 (MQ=255)

                                                                                        cGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCa                                                           >  2:1343045/1‑150 (MQ=255)

                                                                                        cGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCa                                                           <  2:1065733/150‑1 (MQ=255)

                                                                                         gTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCAGCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAg                                                          <  2:1226913/150‑1 (MQ=255)

                                                                                            cATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATa                                                       <  1:1517038/150‑1 (MQ=255)

                                                                                            cATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATa                                                       >  2:665181/1‑150 (MQ=255)

                                                                                              tAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGc                                                     >  1:686630/1‑150 (MQ=255)

                                                                                                aCGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGCCGTCACACCAGCCAAAACCGCTGATCGGTTTACCATCGGCAGTTGAGCTCAGATAGCCg                                                   >  1:1640517/1‑150 (MQ=255)

                                                                                                  ggACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCt                                                 <  1:1270143/150‑1 (MQ=255)

                                                                                                   gACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTg                                                <  2:16624/150‑1 (MQ=255)

                                                                                                     cAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGAACCGGCATAGTCTGCCTGTGCGCCAGCGTCGGGTAGGGTCGCGGGGACGTCAAACCAGCCAAATACGCTGATCGGTTTACCGTCGGCAGTTGAGATCAGAAAGCCTCTGAt                                              >  1:254720/1‑150 (MQ=255)

                                                                                                      aGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATg                                             <  1:802781/150‑1 (MQ=255)

                                                                                                       gTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTGTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGc                                            <  1:331189/150‑1 (MQ=255)

                                                                                                        tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                           <  2:304298/150‑1 (MQ=255)

                                                                                                        tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                           <  2:308071/150‑1 (MQ=255)

                                                                                                        tactacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGcc                                           >  1:1519372/1‑150 (MQ=255)

                                                                                                         actacGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCg                                          >  1:1200377/1‑150 (MQ=255)

                                                                                                              gacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGc                                     <  2:516018/150‑1 (MQ=255)

                                                                                                                cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGa                                   <  1:335998/150‑1 (MQ=255)

                                                                                                                     cACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGc                              >  1:675430/1‑150 (MQ=255)

                                                                                                                      aCGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGGGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCa                             >  2:378599/1‑150 (MQ=255)

                                                                                                                       cgcgAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCAGCGGCAGTTGCGATCAGATAGCCGCAGATGCCGATGGCGAAAAGaaa                            >  1:1618626/1‑150 (MQ=255)

                                                                                                                           aGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAgg                        <  2:876044/150‑1 (MQ=255)

                                                                                                                              cAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTAc                     >  1:1631811/1‑150 (MQ=255)

                                                                                                                              cAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCAACGGCAGTTGAGATCAGATATCCGCTGAAGCCGATGGCGAAAAGCAATACGTAc                     >  2:1106733/1‑150 (MQ=255)

                                                                                                                                 aaaTGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAAt                  >  1:464470/1‑150 (MQ=255)

                                                                                                                                    tGCAGGGCACCGGCAAAGTCTGCCTGTGCGCAAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCg               <  1:64686/150‑1 (MQ=255)

                                                                                                                                     gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga              <  2:130216/150‑1 (MQ=255)

                                                                                                                                     gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga              >  1:1618158/1‑150 (MQ=255)

                                                                                                                                     gCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCga              <  1:737128/150‑1 (MQ=255)

                                                                                                                                      ctgggCACCGGCCAAGTCTGCCTGGGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCgag             <  2:326131/148‑1 (MQ=255)

                                                                                                                                      cAGGGCACCGGCAAAGTCTGCTTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCgag             <  2:1507693/150‑1 (MQ=255)

                                                                                                                                      cAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCgag             >  1:1543640/1‑150 (MQ=255)

                                                                                                                                       ggcgCACAGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCATACCAGCCAAACACGATGATCGGTTTACCATCGGCAGGTTAGCTCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATTCgagg            >  1:244298/4‑149 (MQ=255)

                                                                                                                                       agtgCACAGGAAAAGTCTGCCTGTGAGTCAGAGTGGGCAATTGGCGCGTGGAAGTCAAACCATCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATTCCGATGGCGAAAAGCAATAGGTACAATGCgaga            <  1:1292344/147‑1 (MQ=255)

                                                                                                                                         ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATg          >  2:634778/1‑150 (MQ=255)

                                                                                                                                          gCACAGGCAGACTTTGCCGGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGt         >  1:783820/1‑150 (MQ=255)

                                                                                                                                            aCCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTcc       >  1:1006149/1‑150 (MQ=255)

                                                                                                                                              cGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAACCCAGCCAAAAACGCTGCTCGGTTTACCCTCGGCAGTTGAGATCAGATAGCCGCGGATGCCGATGGCGAAAAGCAAGCGGTCCCATGCGAGATGTcccg     >  1:319164/1‑148 (MQ=255)

                                                                                                                                              cGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAg     >  1:1028162/1‑150 (MQ=255)

                                                                                                                                               ggCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGc    <  2:1375482/150‑1 (MQ=255)

                                                                                                                                               ggCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGc    <  2:1052522/150‑1 (MQ=255)

                                                                                                                                                 cAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACAAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCAGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGctc  <  2:1289983/150‑1 (MQ=255)

TATACTCCATAGTCAGATGACGACTTTCCCAGCATGCGCTTCAGAGTGTCGTCTTTATCGATGAAATGATGCTTCAGGGCCATAAATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACCATCGGCAGTTGAGATCAGATAGCCGCTGATGCCGATGGCGAAAAGCAATAGGTACAATGCGAGATGTCCAGCTC  >  NC_000913/1118481‑1118775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`